1
|
Gardasevic T, Noy A. The impact of sequence periodicity on DNA mechanics: investigating the origin of A-tract's curvature. NANOSCALE 2024. [PMID: 39247956 DOI: 10.1039/d4nr02571g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/10/2024]
Abstract
Periodic sequences in phase with DNA helical shape are prevalent in genomes due to their capacity to modulate DNA elasticity on a global scale. However, how this occurs is not well understood. We use all-atom molecular dynamics simulations on 40 bp DNA fragments to assess the effect of periodicity on bending, twisting, and stretch elasticity. We observe that DNA static curvature is the mechanical parameter most influenced by periodicity, with A-tract sequences having the greatest effect. A-tracts generate global curvature by bending in distinct directions (minor groove and backbones) that complement the bending of the rest of DNA, which predominantly is towards the major groove. Even if A-tracts are rigid at the local scale, these small bends integrate with the greater bends from the sequences between, producing an amplifying effect. As a result, our findings support a 'delocalized bend' model in which the A-tract operates as an 'adaptable mechanical part'. By understanding how global curvature emerges from local fluctuations, we reconcile previous contradictory theories and open an avenue for manipulating DNA mechanics through sequence design.
Collapse
Affiliation(s)
- Tania Gardasevic
- School of Physics, Engineering and Technology, University of York, York, YO10 5DD, UK.
| | - Agnes Noy
- School of Physics, Engineering and Technology, University of York, York, YO10 5DD, UK.
| |
Collapse
|
2
|
Duenas-Herrera M, Bonthuis DJ, Loche P, Netz RR, Scalfi L. Force field for halide and alkali ions in water based on single-ion and ion-pair thermodynamic properties for a wide range of concentrations. J Chem Phys 2024; 161:074506. [PMID: 39158049 DOI: 10.1063/5.0217998] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2024] [Accepted: 07/20/2024] [Indexed: 08/20/2024] Open
Abstract
A classical non-polarizable force field for the common halide (F-, Cl-, Br-, and I-) and alkali (Li+, Na+, K+, and Cs+) ions in SPC/E water is presented. This is an extension of the force field developed by Loche et al. for Na+, K+, Cl-, and Br- (JPCB 125, 8581-8587, 2021): in the present work, we additionally optimize Lennard-Jones parameters for Li+, I-, Cs+, and F- ions. Li+ and F- are particularly challenging ions to model due to their small size. The force field is optimized with respect to experimental solvation free energies and activity coefficients, which are the necessary and sufficient quantities to accurately reproduce the electrolyte thermodynamics. Good agreement with experimental reference data is achieved for a wide range of concentrations (up to 4 mol/l). We find that standard Lorentz-Berthelot combination rules are sufficient for all ions except F-, for which modified combination rules are necessary. With the optimized parameters, we show that, although the force field is only optimized based on thermodynamic properties, structural properties are reproduced quantitatively, while ion diffusion coefficients are in qualitative agreement with experimental values.
Collapse
Affiliation(s)
- Maria Duenas-Herrera
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Douwe Jan Bonthuis
- Institute of Theoretical and Computational Physics, Graz University of Technology, Graz, Austria
| | - Philip Loche
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
- Laboratory of Computational Science and Modeling, IMX, Ecole Polytechnique Federale de Lausanne, 1015 Lausanne, Switzerland
| | - Roland R Netz
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Laura Scalfi
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| |
Collapse
|
3
|
Karatrantos AV, Middendorf M, Nosov DR, Cai Q, Westermann S, Hoffmann K, Nürnberg P, Shaplov AS, Schönhoff M. Diffusion and structure of propylene carbonate-metal salt electrolyte solutions for post-lithium-ion batteries: From experiment to simulation. J Chem Phys 2024; 161:054502. [PMID: 39087537 DOI: 10.1063/5.0216222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2024] [Accepted: 07/11/2024] [Indexed: 08/02/2024] Open
Abstract
The diffusion of cations in organic solvent solutions is important for the performance of metal-ion batteries. In this article, pulsed field gradient nuclear magnetic resonance experiments and fully atomistic molecular dynamic simulations were employed to study the temperature-dependent diffusive behavior of various liquid electrolytes representing 1M propylene carbonate solutions of metal salts with bis(trifluoromethylsulfonyl)imide (TFSI-) or hexafluorophosphate (PF6-) anions commonly used in lithium-ion batteries and beyond. The experimental studies revealed the temperature dependence of the diffusion coefficients for the propylene carbonate (PC) solvent and for the anions following an Arrhenius type of behavior. It was observed that the PC molecules are the faster species. For the monovalent cations (Li+, Na+, K+), the PC solvent diffusion was enhanced as the cation size increased, while for the divalent cations (Mg2+, Ca2+, Sr2+, Ba2+), the opposite trend was observed, i.e., the diffusion coefficients decreased as the cation size increased. The anion diffusion in LiTFSI and NaTFSI solutions was found to be similar, while in electrolytes with divalent cations, a decrease in anion diffusion with increasing cation size was observed. It was shown that non-polarizable charge-scaled force fields could correspond perfectly to the experimental values of the anion and PC solvent diffusion coefficients in salt solutions of both monovalent (Li+, Na+, K+) and divalent (Mg2+, Ca2+, Sr2+, Ba2+) cations at a range of operational temperatures. Finally, after calculating the radial distribution functions between cations, anions, and solvent molecules, the increase in the PC diffusion coefficient established with the increase in cation size for monovalent cations was clearly explained by the large hydration shell of small Li+ cations, due to their strong interaction with the PC solvent. In solutions with larger monovalent cations, such as Na+, and with a smaller solvation shell of PC, the PC diffusion is faster due to more liberated solvent molecules. In the salt solutions with divalent cations, both the anion and the PC diffusion coefficients decreased as the cation size increased due to an enhanced cation-anion coordination, which was accompanied by an increase in the amount of PC in the cation solvation shell due to the presence of anions.
Collapse
Affiliation(s)
- Argyrios V Karatrantos
- Luxembourg Institute of Science and Technology (LIST), 5 avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
- School of Chemistry and Chemical Engineering, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford GU2 7EX, United Kingdom
| | - Maleen Middendorf
- Institute of Physical Chemistry, University of Münster, Corrensstraße 28/30, 48149 Münster, Germany
- International Graduate School on Battery Chemistry, Characterization, Analysis, Recycling and Application (BACCARA), Münster, Germany
| | - Daniil R Nosov
- Luxembourg Institute of Science and Technology (LIST), 5 avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
- Department of Physics and Materials Science, University of Luxembourg, 2 Avenue de l'Université, L-4365 Esch-sur-Alzette, Luxembourg
| | - Qiong Cai
- School of Chemistry and Chemical Engineering, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford GU2 7EX, United Kingdom
| | - Stephan Westermann
- Luxembourg Institute of Science and Technology (LIST), 5 avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
| | - Katja Hoffmann
- Institute of Physical Chemistry, University of Münster, Corrensstraße 28/30, 48149 Münster, Germany
| | - Pinchas Nürnberg
- Institute of Physical Chemistry, University of Münster, Corrensstraße 28/30, 48149 Münster, Germany
| | - Alexander S Shaplov
- Luxembourg Institute of Science and Technology (LIST), 5 avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
| | - Monika Schönhoff
- Institute of Physical Chemistry, University of Münster, Corrensstraße 28/30, 48149 Münster, Germany
| |
Collapse
|
4
|
Vogt J. CO 2 ultrathin film growth on a monolayer of CO 2 adsorbed on the NaCl(100) surface: sticking coefficient and IR-optical signatures in the ν 3 region. Phys Chem Chem Phys 2024; 26:21019-21029. [PMID: 39051430 DOI: 10.1039/d4cp02311k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
CO2 ultrathin molecular films were grown onto a preadsorbed monolayer NaCl(100)/p(2 × 1)-CO2 at 40 K. Polarization infrared spectroscopy (PIRS) reveals that so-prepared films have better quality than directly grown films. A sticking probability of 0.74 ± 0.1 was deduced from the integrated IR absorption. The presence of the monolayer doublet in the film spectra suggests a Stranski-Krastanov film growth with locally varying film thicknesses on the surface. In the region of the ν3(12C16O2) band, fine structure was observed between the well-known transverse-optical (TO) and longitudinal optical (LO) bands. Two independent computational models were applied to analyze the nature of the observed fine structure. Both pair potential calculations in combination with a vibrational exciton model as well as plane-wave density functional theory (DFT) in combination with phonon calculations of IR intensities at the Γ-point reveal that a weak mode visible in s-polarization and p-polarization originates from a vibrational film excitation located near the substrate interface. A series of p-polarized weak bands appearing and partly disappearing upon film-growth is assigned to film stacks of unique local thickness.
Collapse
Affiliation(s)
- Jochen Vogt
- Chemisches Institut der Universität Magdeburg, Universitätsplatz 2, Magdeburg, Germany.
| |
Collapse
|
5
|
Coelho FM, Vinogradov J, Derksen JJ, Franco LFM. Electrokinetic properties of NaCl solution via molecular dynamics simulations with scaled-charge electrolytes. J Chem Phys 2024; 161:044508. [PMID: 39072421 DOI: 10.1063/5.0219098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2024] [Accepted: 06/30/2024] [Indexed: 07/30/2024] Open
Abstract
Scaling ionic charges has become an alternative to polarizable force fields for representing indirect charge transfer effects in molecular simulations. In our work, we apply molecular dynamics simulations to investigate the properties of NaCl aqueous solutions in homogeneous and confined media. We compare classical integer- and scaled-charge force fields for the ions. In the bulk, we validate the force fields by computing equilibrium and transport properties and comparing them with experimental data. Integer-charge ions overestimate dielectric saturation and ionic association. Both force fields present an excess in ion-ion correlation, which leads to a deviation in the ionic conductivity at higher ionic strengths. Negatively charged quartz is used to simulate the confinement effect. Electrostatic interactions dominate counter-ion adsorption. Full-charge ions have stronger and more defined adsorption planes. We obtain the electroosmotic mobility of the solution by combining the shear plane location from non-equilibrium simulations with the ionic distribution from equilibrium simulations. From the Helmholtz-Smoluchowski equation, the zeta potential and the streaming potential coupling coefficient are computed. From an atomic-scale perspective, our molecular dynamics simulations corroborate the hypothesis of maximum packing of the Stern layer, which results in a stable and non-zero zeta potential at high salinity. The scaled-charge model representation of both properties is in excellent qualitative and quantitative agreement with experimental data. With our work, we demonstrate how useful and precise simple scaled-charge models for electrolytes can be to represent complex systems, such as the electrical double layer.
Collapse
Affiliation(s)
- Felipe M Coelho
- Universidade Estadual de Campinas (UNICAMP), Faculdade de Engenharia Química, Campinas-SP 13083-852, Brazil
| | - Jan Vinogradov
- Department of Mechanical Engineering and Mechatronics, Ariel University, 40700 Ariel, Israel
| | - Jos J Derksen
- School of Engineering, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom
| | - Luís F M Franco
- Universidade Estadual de Campinas (UNICAMP), Faculdade de Engenharia Química, Campinas-SP 13083-852, Brazil
| |
Collapse
|
6
|
Cheng K, Xu X, Song J, Chen Y, Kan Z, Li C. Molecular dynamics simulations of choline chloride and ascorbic acid deep eutectic solvents: Investigation of structural and dynamics properties. J Mol Graph Model 2024; 130:108784. [PMID: 38692127 DOI: 10.1016/j.jmgm.2024.108784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2024] [Revised: 04/22/2024] [Accepted: 04/23/2024] [Indexed: 05/03/2024]
Abstract
Deep eutectic solvents (DESs) composed of choline chloride (ChCl) and ascorbic acid (AA) were investigated using the molecular dynamics (MD) simulations. The analyses of the configuration, radial distribution function (RDFs), coordination number, spatial distribution function (SDFs), interaction energies, hydrogen bond number, and self-diffusion coefficient of the ChCl/AA binary systems of different concentrations showed that the stability of the hydrogen bond network and the mutual attraction between systems were the strongest at the experimental eutectic concentration (molar ratio of 2:1). In our simulated temperature range from 303.15 to 353.15 K, the hydrogen bonding network of ChCl/AA DES does not undergo considerable alterations, indicating that its stability was insensitive to temperature. In addition, the influence of the water content on the ChCl/AA DES system was further investigated. The simulated results revealed that the water molecules could disrupt the formation of the hydrogen bonding network by occupyin positions that are essential for the formation of hydrogen bonds within the DES system.
Collapse
Affiliation(s)
- Kexin Cheng
- Department of Chemistry, School of Science, China Pharmaceutical University, Nanjing, 211198, Jiangsu, PR China
| | - Xuchen Xu
- Department of Chemistry, School of Science, China Pharmaceutical University, Nanjing, 211198, Jiangsu, PR China
| | - Jilun Song
- Department of Chemistry, School of Science, China Pharmaceutical University, Nanjing, 211198, Jiangsu, PR China
| | - Yu Chen
- Department of Chemistry, School of Science, China Pharmaceutical University, Nanjing, 211198, Jiangsu, PR China
| | - Zigui Kan
- Department of Chemistry, School of Science, China Pharmaceutical University, Nanjing, 211198, Jiangsu, PR China.
| | - Caolong Li
- Department of Chemistry, School of Science, China Pharmaceutical University, Nanjing, 211198, Jiangsu, PR China.
| |
Collapse
|
7
|
O’Neill N, Shi BX, Fong K, Michaelides A, Schran C. To Pair or not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water. J Phys Chem Lett 2024; 15:6081-6091. [PMID: 38820256 PMCID: PMC11181334 DOI: 10.1021/acs.jpclett.4c01030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 05/23/2024] [Accepted: 05/24/2024] [Indexed: 06/02/2024]
Abstract
The extent of ion pairing in solution is an important phenomenon to rationalize transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However, direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wave function based machine learning potentials with the random phase approximation (RPA) and second order Møller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.
Collapse
Affiliation(s)
- Niamh O’Neill
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Cavendish
Laboratory, Department of Physics, University
of Cambridge, Cambridge CB3 0HE, United
Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Benjamin X. Shi
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Kara Fong
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Christoph Schran
- Cavendish
Laboratory, Department of Physics, University
of Cambridge, Cambridge CB3 0HE, United
Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| |
Collapse
|
8
|
Ishraaq R, Das S. All-atom molecular dynamics simulations of polymer and polyelectrolyte brushes. Chem Commun (Camb) 2024; 60:6093-6129. [PMID: 38819435 DOI: 10.1039/d4cc01557f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Densely grafted polymer and polyelectrolyte (PE) brushes, owing to their significant abilities to functionalize surfaces for a plethora of applications in sensing, diagnostics, current rectification, surface wettability modification, drug delivery, and oil recovery, have attracted significant attention over the past several decades. Unfortunately, most of the attention has primarily focused on understanding the properties of the grafted polymer and the PE chains with little attention devoted to studying the behavior of the brush-supported ions (counterions needed to screen the PE chains) and water molecules. Over the past few years, our group has been at the forefront of addressing this gap: we have employed all-atom molecular dynamics (MD) simulations for studying a wide variety of polymer and PE brush systems with specific attention to unraveling the properties and behavior of the brush-supported water molecules and ions. Our findings have revealed some of the most fascinating properties of such brush-supported ions and water molecules, including the most remarkable control of nanofluidic transport afforded by the specific ion and water responses induced by the PE brushes grafted on the inner walls of the nanochannel. This feature article aims to summarize some of our key contributions associated with such atomistic simulations of polymer and PE brushes and brush-supported water molecules and counterions.
Collapse
Affiliation(s)
- Raashiq Ishraaq
- Department of Mechanical Engineering, University of Maryland, College Park, MD 20742, USA.
| | - Siddhartha Das
- Department of Mechanical Engineering, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
9
|
Córdoba A, Montes de Oca JM, Darling SB, de Pablo JJ. Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores. ACS NANO 2024; 18:12569-12579. [PMID: 38696274 DOI: 10.1021/acsnano.4c03546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2024]
Abstract
In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. These findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.
Collapse
Affiliation(s)
- Andrés Córdoba
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Joan Manuel Montes de Oca
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Seth B Darling
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Juan J de Pablo
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
10
|
Goswami A, Peña-Torres A, Jónsson EÖ, Egorov SA, Jónsson H. Evidence of Sharp Transitions between Octahedral and Capped Trigonal Prism States of the Solvation Shell of the Aqueous Fe 3+ Ion. J Phys Chem Lett 2024; 15:4523-4530. [PMID: 38634894 DOI: 10.1021/acs.jpclett.4c00756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2024]
Abstract
The structure of the solvation shell of the aqueous Fe3+ ion has been a subject of controversy due to discrepancies between experiments and different levels of theory. We address this issue by performing simulations for a wide range of ion concentrations, using various potential energy functions, supplemented by density functional theory calculations of selected configurations. The solvation shell undergoes abrupt transitions between two states: a hexacoordinated octahedral (OH) state and a capped trigonal prism (CTP) state with 7-fold coordination. The lifetime of these states is dependent on concentration. In dilute FeCl3 solutions, the lifetimes of both are similar (≈1 ns). However, the lifetime of the OH state increases with ion concentration, while that of the CTP state decreases slightly. When a uniform negative background charge is used instead of explicit counterions, the lifetime of the OH state is greatly overestimated. These findings underscore the need for further experimental measurements and higher-level simulations.
Collapse
Affiliation(s)
- Amrita Goswami
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Alejandro Peña-Torres
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Elvar Ö Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Sergei A Egorov
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
- Department of Chemistry, University of Virginia, Charlottesville, Virginia 22901, United States
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| |
Collapse
|
11
|
Ma L, Jiang J. Vehicular Motions Dominate the Ion Transport in Concentrated LiTFSI Aqueous Solutions? J Phys Chem Lett 2024; 15:4531-4537. [PMID: 38635898 DOI: 10.1021/acs.jpclett.4c00791] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/20/2024]
Abstract
Water-in-salt electrolytes (WiSEs) show great promise for applications in grid-scale energy storage. The design of high-performance WiSEs requires a comprehensive understanding of their microstructures and ion transport properties. In the present work, based on the CL&Pol force field, we have developed a polarizable force field (PFF) tailored for high-concentration LiTFSI aqueous solutions, which accurately reproduces the structural and dynamical properties. Unlike the literature, we do not observe the presence of bulk-like water in LiTFSI solutions exceeding 19 mol/kg. Furthermore, we find that the vast majority of Li(H20)n+ are short-lived, and thus, the structural motion rather than the vehicular motion is the main mode of ion transport. Our results have significant implications for understanding the ion dynamics in WiSEs. Additionally, further in-depth experimental analyses are imperative.
Collapse
Affiliation(s)
- Linbo Ma
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Jian Jiang
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| |
Collapse
|
12
|
Summa CM, Langford DP, Dinshaw SH, Webb J, Rick SW. Calculations of Absolute Free Energies, Enthalpies, and Entropies for Drug Binding. J Chem Theory Comput 2024; 20:2812-2819. [PMID: 38538531 DOI: 10.1021/acs.jctc.4c00057] [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: 04/10/2024]
Abstract
Computer simulation methods can aid in the rational design of drugs aimed at a specific target, typically a protein. The affinity of a drug for its target is given by the free energy of binding. Binding can be further characterized by the enthalpy and entropy changes in the process. Methods exist to determine exact free energies, enthalpies, and entropies that are dependent only on the quality of the potential model and adequate sampling of conformational degrees of freedom. Entropy and enthalpy are roughly an order of magnitude more difficult to calculate than the free energy. This project combines a replica exchange method for enhanced sampling, designed to be efficient for protein-sized systems, with free energy calculations. This approach, replica exchange with dynamical scaling (REDS), uses two conventional simulations at different temperatures so that the entropy can be found from the temperature dependence of the free energy. A third replica is placed between them, with a modified Hamiltonian that allows it to span the temperature range of the conventional replicas. REDS provides temperature-dependent data and aids in sampling. It is applied to the bromodomain-containing protein 4 (BRD4) system. We find that for the force fields used, the free energies are accurate but the entropies and enthalpies are not, with the entropic contribution being too positive. Reproducing the entropy and enthalpy of binding appears to be a more stringent test of the force fields than reproducing the free energy.
Collapse
Affiliation(s)
- Christopher M Summa
- Department of Computer Science, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Dillon P Langford
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Sam H Dinshaw
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Jennifer Webb
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Steven W Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| |
Collapse
|
13
|
Yang Y, Narayanan Nair AK, Lau D, Sun S. Interfacial properties of the brine + carbon dioxide + oil + silica system. J Chem Phys 2024; 160:114702. [PMID: 38497476 DOI: 10.1063/5.0197087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 03/01/2024] [Indexed: 03/19/2024] Open
Abstract
Molecular dynamics simulations of the H2O + CO2 + aromatic hydrocarbon and H2O + CO2 + benzene + silica (hydrophilic) systems are performed to gain insights into CO2-enhanced oil recovery (EOR) processes. For comparison purposes, an overview of the previous simulation studies of the interfacial properties of the brine + CO2 + alkane + silica system is also presented. In general, the water contact angle (CA) of the H2O + CO2 + silica (hydrophilic) system increased with pressure and decreased with temperature. The CAs of the H2O + hydrocarbon + silica (hydrophilic) system are not significantly affected by temperature and pressure. The simulated CAs were in the ranges of about 58°-77° and 81°-93° for the H2O + hexane + silica (hydrophilic) and the H2O + aromatic hydrocarbon + silica (hydrophilic) systems, respectively. In general, these CAs were not significantly influenced by the addition of CO2. The simulated CAs were in the ranges of about 51.4°-95.0°, 69.1°-86.0°, and 72.0°-87.9° for the brine + CO2 + silica (hydrophilic), brine + hexane + silica (hydrophilic), and brine + CO2 + hexane + silica (hydrophilic) systems, respectively. All these CAs increased with increasing NaCl concentration. The adhesion tension of the brine + silica (hydrophilic) system in the presence of CO2 and/or hexane decreased with increasing salt concentration. The simulated CAs were in the range of about 117°-139° for the H2O + alkane + silica (hydrophobic) system. These CAs are increased by the addition of CO2. At high pressures, the distributions of H2O normal to the silica (hydrophobic) surface in the droplet region of the H2O + silica system were found to be strongly affected by the presence of CO2. These insights might be key for optimizing the performance of the miscible CO2 water-alternating-gas injection schemes widely used for EOR.
Collapse
Affiliation(s)
- Yafan Yang
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou 221116, Jiangsu, China
| | - Arun Kumar Narayanan Nair
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Denvid Lau
- Department of Architecture and Civil Engineering, City University of Hong Kong, Hong Kong 999077, China
| | - Shuyu Sun
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
14
|
Mondal S, Bagchi B. Ion diffusion captures composition-dependent anomalies in water-DMSO binary mixtures. J Chem Phys 2024; 160:114505. [PMID: 38506289 DOI: 10.1063/5.0189259] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 03/01/2024] [Indexed: 03/21/2024] Open
Abstract
Aqueous dimethyl sulfoxide (Aq-DMSO) binary mixture exhibits many fascinating composition-dependent anomalies that are explained by using the peculiarities of the water-DMSO hydrogen bond. Ions can couple strongly to these composition-dependent anomalies to produce exotic dynamics of their own. We carry out theoretical studies using computer simulations to understand the structural and dynamical aspects of rigid monovalent cations (Li+, Na+, K+, Rb+, and Cs+) in aqueous DMSO solutions, with chloride as the counterion. We uncover a number of composition-dependent ion diffusion anomalies, which can be traced back to the interplay between the size-dependent charge density of the ion and the resulting difference in interactions of the ion with water and DMSO molecules. Size and composition dependence of the diffusion coefficients of the five ions exhibit fascinating variations that can be explained partially.
Collapse
Affiliation(s)
- Sangita Mondal
- SSCU, Indian Institute of Science, Bangalore 560012, India
| | - Biman Bagchi
- SSCU, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
15
|
Xie M, Zhang M, Jin Z. Machine Learning-Based Interfacial Tension Equations for (H 2 + CO 2)-Water/Brine Systems over a Wide Range of Temperature and Pressure. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:5369-5377. [PMID: 38417158 DOI: 10.1021/acs.langmuir.3c03831] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/01/2024]
Abstract
Large-scale underground hydrogen storage (UHS) plays a vital role in energy transition. H2-brine interfacial tension (IFT) is a crucial parameter in structural trapping in underground geological locations and gas-water two-phase flow in subsurface porous media. On the other hand, cushion gas, such as CO2, is often co-injected with H2 to retain reservoir pressure. Therefore, it is imperative to accurately predict the (H2 + CO2)-water/brine IFT under UHS conditions. While there have been a number of experimental measurements on H2-water/brine and (H2 + CO2)-water/brine IFT, an accurate and efficient (H2 + CO2)-water/brine IFT model under UHS conditions is still lacking. In this work, we use molecular dynamics (MD) simulations to generate an extensive (H2 + CO2)-water/brine IFT databank (840 data points) over a wide range of temperature (from 298 to 373 K), pressure (from 50 to 400 bar), gas composition, and brine salinity (up to 3.15 mol/kg) for typical UHS conditions, which is used to develop an accurate and efficient machine learning (ML)-based IFT equation. Our ML-based IFT equation is validated by comparing to available experimental data and other IFT equations for various systems (H2-brine/water, CO2-brine/water, and (H2 + CO2)-brine/water), rendering generally good performance (with R2 = 0.902 against 601 experimental data points). The developed ML-based IFT equation can be readily applied and implemented in reservoir simulations and other UHS applications.
Collapse
Affiliation(s)
- Minjunshi Xie
- School of Mining and Petroleum Engineering, Department of Civil and Environmental Engineering, University of Alberta, Edmonton, AB T6G 1H9, Canada
| | - Mingshan Zhang
- Key Laboratory of Ministry of Education on Safe Mining of Deep Metal Mines, School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China
- Key Laboratory of Liaoning Province on Deep Engineering and Intelligent Technology, Northeastern University, Shenyang 110819, China
| | - Zhehui Jin
- School of Mining and Petroleum Engineering, Department of Civil and Environmental Engineering, University of Alberta, Edmonton, AB T6G 1H9, Canada
| |
Collapse
|
16
|
Xiao T, Zhou Y, Li B. Energy-Scaled Debye-Hückel Theory for the Electrostatic Solvation Free Energy in Size-Asymmetric Electrolyte Solutions. J Phys Chem B 2024; 128:1029-1039. [PMID: 38235680 DOI: 10.1021/acs.jpcb.3c07233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
In this report, an energy-scaled Debye-Hückel theory is developed for fast and accurate evaluation of the electrostatic solvation free energy in size-asymmetric electrolyte solutions. A size-asymmetric electrolyte solution is mapped to a dielectric continuum medium with Debye-Hückel-like response. Based on the scaling relation of the electrostatic energy of a spherical ion in the small and large size limits, a Padé polynomial is used to interpolate the electrostatic energy at finite size. The Padé polynomial is further interpreted as the electrostatic energy of an effective Debye-Hückel mean field model, depicted by a modified Debye parameter and a surface charge density due to the size asymmetry of the solvent ions. This theory can distinguish the electrostatic energies and the electrostatic solvation free energies of solutes with the same size but opposite charges. Application to charged hard and charged soft spheres demonstrates the accuracy of our approach.
Collapse
Affiliation(s)
- Tiejun Xiao
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, Guiyang 550018, People's Republic of China
| | - Yun Zhou
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, Guiyang 550018, People's Republic of China
| | - Bo Li
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, Guiyang 550018, People's Republic of China
| |
Collapse
|
17
|
Blazquez S, Bourg IC, Vega C. Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2. J Chem Phys 2024; 160:046101. [PMID: 38270237 DOI: 10.1063/5.0186233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2023] [Accepted: 12/27/2023] [Indexed: 01/26/2024] Open
Abstract
In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Ian C Bourg
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, USA
- High Meadows Environmental Institute, Princeton University, Princeton, New Jersey 08544, USA
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
18
|
Li X, Bourg IC. Hygroscopic Growth of Adsorbed Water Films on Smectite Clay Particles. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2024; 58:1109-1118. [PMID: 38164899 PMCID: PMC10795194 DOI: 10.1021/acs.est.3c08253] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2023] [Revised: 12/07/2023] [Accepted: 12/12/2023] [Indexed: 01/03/2024]
Abstract
Hygroscopic growth of adsorbed water films on clay particles underlies a number of environmental science questions, from the air quality and climate impacts of mineral dust aerosols to the hydrology and mechanics of unsaturated soils and sedimentary rocks. Here, we use molecular dynamics (MD) simulations to establish the relation between adsorbed water film thickness (h) and relative humidity (RH) or disjoining pressure (Π), which has long been uncertain due to factors including sensitivity to particle shape, surface roughness, and aqueous chemistry. We present a new MD simulation approach that enables precise quantification of Π in films up to six water monolayers thick. We find that the hygroscopicity of phyllosilicate mineral surfaces increases in the order mica < K-smectite < Na-smectite. The relationship between Π and h on clay surfaces follows a double exponential decay with e-folding lengths of 2.3 and 7.5 Å. The two decay length scales are attributed to hydration repulsion and osmotic phenomena in the electrical double layer (EDL) at the clay-water interface.
Collapse
Affiliation(s)
- Xiaohan Li
- Department
of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ian C. Bourg
- Department
of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
- High
Meadows Environmental Institute, Princeton
University, Princeton, New Jersey 08544, United States
| |
Collapse
|
19
|
Li K, Chen B, Li M, Jiang L, Song Y, Yang M. Facilitation of Hydrate Dissociation and Structural Evolution by Major Marine Anions under Static Electric Fields. J Phys Chem B 2023; 127:10447-10457. [PMID: 37991934 DOI: 10.1021/acs.jpcb.3c06012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2023]
Abstract
Electric fields have been proven to be capable of significantly affecting the equilibrium state of hydrates. In this study, the thermodynamic properties and structural changes of methane hydrate (MH) in various anion solutions in an electric field at 0.7 V/nm were investigated by molecular dynamics simulations. The presence of anions significantly enhances the instability of methane hydrates under electric fields, leading to a staged dissociation process. First, the anions coexist with MH to form a temporary metastable structure under the action of an electric field. Then, the migration of anions causes the dissociation of nearby hydrates and the formation of flow channels in the hydrate layer, which leads to the complete dissociation of MH after a period. The promotive effects of F-, Br-, I-, and Cl- ions were close, while SO42- was relatively weak. The anions are still in hydration shells in the MH phase, but the structure of the hydration shells differs slightly from that in solution (the coordination numbers of I- and SO42- ions increased). The migration resistances of multiple anions to cross the surface of the hydrate layer are similar. However, inside the hydrate phase, the anions with a larger radius have a higher migration resistance. It is difficult for SO42- ions to migrate inside the hydrate phase, and they tend to form a metastable structure on the hydrate surface. Combining our previous studies, SrCl2 solution has the best hydrate promotion under an electric field environment.
Collapse
Affiliation(s)
- Kehan Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Bingbing Chen
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Mingjun Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Lanlan Jiang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Yongchen Song
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Mingjun Yang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| |
Collapse
|
20
|
Seltzer AM, Shackleton SA, Bourg IC. Solubility Equilibrium Isotope Effects of Noble Gases in Water: Theory and Observations. J Phys Chem B 2023; 127:9802-9812. [PMID: 37937341 DOI: 10.1021/acs.jpcb.3c05651] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2023]
Abstract
The abundance and isotopic composition of noble gases dissolved in water have many applications in the geosciences. In recent years, new analytical techniques have opened the door to the use of high-precision measurements of noble gas isotopes as tracers for groundwater hydrology, oceanography, mantle geochemistry, and paleoclimatology. These analytical advances have brought about new measurements of solubility equilibrium isotope effects (SEIEs) in water (i.e., the relative solubilities of noble gas isotopes) and their sensitivities to the temperature and salinity. Here, we carry out a suite of classical molecular dynamics (MD) simulations and employ the theoretical method of quantum correction to estimate SEIEs for comparison with experimental observations. We find that classical MD simulations can accurately predict SEIEs for the isotopes of Ar, Kr, and Xe to order 0.01‰, on the scale of analytical uncertainty. However, MD simulations consistently overpredict the SEIEs of Ne and He by up to 40% of observed values. We carry out sensitivity tests at different temperatures, salinities, and pressures and employ different sets of interatomic potential parameters and water models. For all noble gas isotopes, the TIP4P water model is found to reproduce observed SEIEs more accurately than the SPC/E and TIP4P/ice models. Classical MD simulations also accurately capture the sign and approximate magnitude of temperature and salinity sensitivities of SEIEs for heavy noble gases. We find that experimental and modeled SEIEs generally follow an inverse-square mass dependence, which implies that the mean-square force experienced by a noble gas atom within a solvation shell is similar for all noble gases. This inverse-square mass proportionality is nearly exact for Ar, Kr, and Xe isotopes, but He and Ne exhibit a slightly weaker mass dependence. We hypothesize that the apparent dichotomy between He-Ne and Ar-Kr-Xe SEIEs may result from atomic size differences, whereby the smaller noble gases are more likely to spontaneously fit within cavities of water without breaking water-water H-bonds, thereby experiencing softer collisions during translation within a solvation shell. We further speculate that the overprediction of simulated He and Ne SEIEs may result from the neglection of higher-order quantum corrections or the overly stiff representation of van der Waals repulsion by the widely used Lennard-Jones 6-12 potential model. We suggest that new measurements of SEIEs of heavy and light noble gases may represent a novel set of constraints with which to refine hydrophobic solvation theories and optimize the set of interatomic potential models used in MD simulations of water and noble gases.
Collapse
Affiliation(s)
- Alan M Seltzer
- Woods Hole Oceanographic Institution, Woods Hole, Massachusetts 02543-1050,United States
| | - Sarah A Shackleton
- Department of Geosciences, Princeton University, Princeton, New Jersey 08544,United States
| | - Ian C Bourg
- Department of Civil and Environmental Engineering (CEE) and High Meadows Environmental Institute (HMEI), Princeton University, Princeton, New Jersey 08544,United States
| |
Collapse
|
21
|
Villamil V, Vairoletti F, Tijman A, López G, Peixoto de Abreu Lima A, Saiz C, Iglesias C, Mahler G. Novel Kinetic Resolution of Thiazolo-Benzimidazolines Using MAO Enzymes. ACS OMEGA 2023; 8:42114-42125. [PMID: 38024698 PMCID: PMC10652373 DOI: 10.1021/acsomega.3c03223] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 10/06/2023] [Accepted: 10/10/2023] [Indexed: 12/01/2023]
Abstract
The kinetic resolution of racemic 1H,3H-thiazolo[3,4-a]benzimidazoline (TBIM) heterocycles was achieved using E. coli whole cells expressing the MAO-N D11 enzyme. Several cosolvents were screened using TBIM 2a as the substrate. DMF was the best cosolvent, affording the pure enantiomer (+)-2a in 44% yield, 94% ee. The stereochemistry of TBIM was predicted by means of ab initio calculations of optical rotation and circular dichroism spectra. The reaction scope was investigated for 11 substituted (±) TBIM using an optimized protocol. The best yield and % ee were obtained for the nonsubstituted 2a. Among the substituted compounds, the 5-substituted-TBIM showed better % ee than the 4-substituted one. The small electron donor group (Me) led to better % ee than the electron-withdrawing groups (-NO2 and -CO2Et), and the bulky naphthyl group was detrimental for the kinetic resolution. Docking experiments and molecular dynamics (MD) simulations were employed to further understand the interactions between MAO-N D11 and the thiazolo-benzimidazoline substrates. For 2a, the MD showed favorable positioning and binding energy for both enantiomers, thus suggesting that this kinetic resolution is influenced not only by the active site but also by the entry tunnel. This work constitutes the first report of the enzymatic kinetic resolution applied to TBIM heterocycles.
Collapse
Affiliation(s)
- Valentina Villamil
- Departamento
de Química Orgánica, Laboratorio de Quimica Farmaceutica,
Facultad de Quimica, Universidad de la República, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
| | - Franco Vairoletti
- Departamento
de Química Orgánica, Laboratorio de Quimica Farmaceutica,
Facultad de Quimica, Universidad de la República, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
- Programa
de Posgrado en Quimica, Universidad de la
República Uruguay, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
| | - Ariel Tijman
- Programa
de Posgrado en Quimica, Universidad de la
República Uruguay, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
- Departamento
de Biociencias, Laboratorio de Microbiología Molecular, Facultad
de Quimica, Universidad de la Republica, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
- Departamento
de Biociencias y Departamento de Quimica Organica, Laboratorio de
Biocatalisis y Biotransformaciones, Facultad de Quimica, Universidad de la Republica, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
| | - Gonzalo López
- Programa
de Posgrado en Quimica, Universidad de la
República Uruguay, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
- Departamento
de Biociencias, Laboratorio de Microbiología Molecular, Facultad
de Quimica, Universidad de la Republica, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
- Departamento
de Biociencias y Departamento de Quimica Organica, Laboratorio de
Biocatalisis y Biotransformaciones, Facultad de Quimica, Universidad de la Republica, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
| | - Alejandro Peixoto de Abreu Lima
- Departamento
de Química Orgánica, Laboratorio de Síntesis
Orgánica, Facultad de Quimica, Universidad
de la Republica, Gral
Flores 2124, Montevideo, Montevideo 11800, Uruguay
| | - Cecilia Saiz
- Departamento
de Química Orgánica, Laboratorio de Quimica Farmaceutica,
Facultad de Quimica, Universidad de la República, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
| | - César Iglesias
- Departamento
de Biociencias, Laboratorio de Microbiología Molecular, Facultad
de Quimica, Universidad de la Republica, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
- Departamento
de Biociencias y Departamento de Quimica Organica, Laboratorio de
Biocatalisis y Biotransformaciones, Facultad de Quimica, Universidad de la Republica, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
| | - Graciela Mahler
- Departamento
de Química Orgánica, Laboratorio de Quimica Farmaceutica,
Facultad de Quimica, Universidad de la República, Gral Flores 2124, Montevideo, Montevideo 11800, Uruguay
| |
Collapse
|
22
|
Sun Y, Zhang D, Bashir A, Li C, Fan Z. Scaling Solute-Solvent Distances to Improve Solubility and Ion Paring Predictions in Rigid Ion Models. J Phys Chem B 2023; 127:9575-9586. [PMID: 37906589 DOI: 10.1021/acs.jpcb.3c05993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
Force fields based on the rigid ion model (RIM) have been developed to accurately predict the various physical and chemical properties of salts and water. However, the combined use of these models often fails to accurately predict the solubility of salts in water. To address this issue, several approaches, such as charge scaling or reparameterization, have been proposed. Nevertheless, these methods require laborious reparameterization of nonbonded force field parameters. In this article, we propose a scaling solute-solvent distance (SSSD) method to improve force fields in predicting salt solubility without changing the solute-solute and solvent-solvent interactions in the original force fields. This method can also tune the ion pairing of salt in water. One main advantage of the SSSD method is that reparameterization of the crystal and water models is not needed. We use two RIMs for the NaCl-water system (JC-SPC/E and SD-SPC/E) and the CHARMM force field for the KCl-water system to demonstrate the improved accuracy in predicting solubility by the SSSD method. Furthermore, we use the RDG-SPC/Fw force field to show that the SSSD method can also be used to tune the ion pairing of CaCO3 in water. Limitations of this method are also discussed.
Collapse
Affiliation(s)
- Ying Sun
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
| | - Dan Zhang
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- Department of Chemistry, Shanghai University, Shanghai 200444, P. R. China
| | - Ayesha Bashir
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- School of Nano-Tech and Nano-Bionics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Chuncheng Li
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
| | - Zhaochuan Fan
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- School of Nano-Tech and Nano-Bionics, University of Science and Technology of China, Hefei 230026, P. R. China
| |
Collapse
|
23
|
Zhang J, Clennell MB, Dewhurst DN. Transport Properties of NaCl in Aqueous Solution and Hydrogen Solubility in Brine. J Phys Chem B 2023; 127:8900-8915. [PMID: 37794729 PMCID: PMC10591480 DOI: 10.1021/acs.jpcb.3c03863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Revised: 08/17/2023] [Indexed: 10/06/2023]
Abstract
Ion transport properties and hydrogen solubility in brine play pivotal roles in various engineering and scientific scopes including chemical, physical, geochemical, and geothermal domains. Molecular dynamics simulations were performed to obtain transport properties of NaCl in the binary H2O + NaCl system using different force fields. Brine density, ion diffusivity, molar conductivity, conductivity, and hydrogen solubilities were obtained as functions of temperature and salt concentration. We compared the performance of different force fields against the experimental correlation model and developed three mathematical models. The first was the modified brine density model based on the simulated brine density over a wide range of salinity levels, and the second and third analytical mathematical models were derived for the ion diffusivity and molar conductivity as a function of salinity and temperature. The results of this study illustrated that the modified brine density model not only produced the same results of the previous model for lower salinity levels but also applied well to predict the brine density for a higher salinity level. The derived mathematical models indicated that the ion diffusivity and molar conductivity decreased linearly with salinity, and the slope and y-intercept of the lines of diffusivity and molar conductivity versus temperature were third-order polynomials of temperature. The developed models provided the mechanism for the behavior of decreasing molar conductivity with increasing salinity and increasing conductivity with increasing salinity. The directions of the effect of salinity on the molar conductivity and conductivity were opposite. The molar conductivity increased with a decreasing salinity level. However, the conductivity increased with increasing salinity, as the contribution of the ion concentration or salinity level to conductivity dominated over that of the ion movement.
Collapse
Affiliation(s)
- Junfang Zhang
- CSIRO Energy, 26 Dick Perry Avenue, Kensington, Western Australia 6151, Australia
| | - Michael B. Clennell
- CSIRO Energy, 26 Dick Perry Avenue, Kensington, Western Australia 6151, Australia
| | - David N. Dewhurst
- CSIRO Energy, 26 Dick Perry Avenue, Kensington, Western Australia 6151, Australia
| |
Collapse
|
24
|
Zhang Q, Li J. Benchmark computational investigations for the basic model of the salt-water complex: NaCl(H 2O) and its anion NaCl(H 2O) . Phys Chem Chem Phys 2023; 25:27215-27229. [PMID: 37791409 DOI: 10.1039/d3cp03421f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The microsolvation of salts in water is a fundamental physicochemical process. In this work, the aqueous salt complex NaCl(H2O) and its anion NaCl(H2O)- were investigated using comprehensive calculations, including the costly and accurate CCSD(T)-F12a and focal point analysis (FPA) methods. For the neutral NaCl(H2O), three isomers exist, two of which are mirror-symmetric with almost identical structures and their corresponding anions are also mirror-symmetric. For the NaCl(H2O)- anion, there are four isomers. Several transition states are found for the first time. The structural rearrangements of neutral NaCl(H2O) and NaCl(H2O)- anions are mainly caused by breaking and forming of the hydrogen bonds and the enhancement and weakening of interactions between Na and O atoms. The distributions of the anion complexes from 15-300 K are computed and compared to recent experimental results. The analysis of the intermolecular weak interactions shows the weak van der Waals interactions between Na and O atoms, as well as hydrogen bonding between H and Cl. Moreover, the theoretically predicted anion photoelectron spectra are assigned and analyzed in detail, and they agree with experimental spectra satisfactorily. The Na-Cl stretching vibrational mode dominates the vibrational structure in both anion spectra with some minor contributions from the intermolecular motions between H2O and NaCl.
Collapse
Affiliation(s)
- Qi Zhang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| |
Collapse
|
25
|
Zheng X, Bourg IC. Nanoscale Prediction of the Thermal, Mechanical, and Transport Properties of Hydrated Clay on 10 6- and 10 15-Fold Larger Length and Time Scales. ACS NANO 2023; 17:19211-19223. [PMID: 37774155 PMCID: PMC10569101 DOI: 10.1021/acsnano.3c05751] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/25/2023] [Accepted: 09/25/2023] [Indexed: 10/01/2023]
Abstract
Coupled thermal, hydraulic, mechanical, and chemical (THMC) processes, such as desiccation-driven cracking or chemically driven fluid flow, significantly impact the performance of composite materials formed by fluid-mediated nanoparticle assembly, including energy storage materials, ordinary Portland cement, bioinorganic nanocomposites, liquid crystals, and engineered clay barriers used in the isolation of hazardous wastes. These couplings are particularly important in the isolation of high-level radioactive waste (HLRW), where heat generated by radioactive decay can drive the temperature up to at least 373 K in the engineered barrier. Here, we use large-scale all-atom molecular dynamics simulations of hydrated smectite clay nanoparticle assemblages to predict the fundamental THMC properties of hydrated compacted clay over a wide range of temperatures (up to 373 K) and dry densities relevant to HLRW management. Equilibrium simulations of clay-water mixtures at different hydration levels are analyzed to quantify material properties, including thermal conductivity, heat capacity, thermal expansion, suction, water and ion self-diffusivity, and hydraulic conductivity. Predictions are validated against experimental results for the properties of compacted bentonite clay. Our results demonstrate the feasibility of using atomistic-level simulations of assemblages of clay nanoparticles on scales of tens of nanometers and nanoseconds to infer the properties of compacted bentonite on scales of centimeters and days, a direct upscaling over 6 orders of magnitude in space and 15 orders of magnitude in time.
Collapse
Affiliation(s)
- Xiaojin Zheng
- Department
of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ian C. Bourg
- Department
of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
- High
Meadows Environmental Institute, Princeton
University, Princeton, New Jersey 08544, United States
| |
Collapse
|
26
|
Omrani S, Ghasemi M, Singh M, Mahmoodpour S, Zhou T, Babaei M, Niasar V. Interfacial Tension-Temperature-Pressure-Salinity Relationship for the Hydrogen-Brine System under Reservoir Conditions: Integration of Molecular Dynamics and Machine Learning. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:12680-12691. [PMID: 37650690 PMCID: PMC10501201 DOI: 10.1021/acs.langmuir.3c01424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 08/14/2023] [Indexed: 09/01/2023]
Abstract
Hydrogen (H2) underground storage has attracted considerable attention as a potentially efficient strategy for the large-scale storage of H2. Nevertheless, successful execution and long-term storage and withdrawal of H2 necessitate a thorough understanding of the physical and chemical properties of H2 in contact with the resident fluids. As capillary forces control H2 migration and trapping in a subsurface environment, quantifying the interfacial tension (IFT) between H2 and the resident fluids in the subsurface is important. In this study, molecular dynamics (MD) simulation was employed to develop a data set for the IFT of H2-brine systems under a wide range of thermodynamic conditions (298-373 K temperatures and 1-30 MPa pressures) and NaCl salinities (0-5.02 mol·kg-1). For the first time to our knowledge, a comprehensive assessment was carried out to introduce the most accurate force field combination for H2-brine systems in predicting interfacial properties with an absolute relative deviation (ARD) of less than 3% compared with the experimental data. In addition, the effect of the cation type was investigated for brines containing NaCl, KCl, CaCl2, and MgCl2. Our results show that H2-brine IFT decreases with increasing temperature under any pressure condition, while higher NaCl salinity increases the IFT. A slight decrease in IFT occurs when the pressure increases. Under the impact of cation type, Ca2+ can increase IFT values more than others, i.e., up to 12% with respect to KCl. In the last step, the predicted IFT data set was used to provide a reliable correlation using machine learning (ML). Three white-box ML approaches of the group method of data handling (GMDH), gene expression programming (GEP), and genetic programming (GP) were applied. GP demonstrates the most accurate correlation with a coefficient of determination (R2) and absolute average relative deviation (AARD) of 0.9783 and 0.9767%, respectively.
Collapse
Affiliation(s)
- Sina Omrani
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| | - Mehdi Ghasemi
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| | - Mrityunjay Singh
- Institute
of Applied Geosciences, Geothermal Science and Technology, Technische Universität Darmstadt, 64289 Darmstadt, Germany
| | - Saeed Mahmoodpour
- Group
of Geothermal Technologies, Technische Universität
Munchen, 80333 Munich, Germany
| | - Tianhang Zhou
- College
of Carbon Neutrality Future Technology, China University of Petroleum (Beijing), 102249 Beijing, China
| | - Masoud Babaei
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| | - Vahid Niasar
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| |
Collapse
|
27
|
Zeng L, Chen M, Wang Z, Qiao R, Feng G. Structural Evolution Governs Reversible Heat Generation in Electrical Double Layers. PHYSICAL REVIEW LETTERS 2023; 131:096201. [PMID: 37721826 DOI: 10.1103/physrevlett.131.096201] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Revised: 05/01/2023] [Accepted: 07/07/2023] [Indexed: 09/20/2023]
Abstract
Electrical double layer (EDL) formation determines the reversible heat generation of supercapacitors. While classical theories suggest an exothermic nature, experiments revealed that it can be endothermic, depending on the polarization and electrolyte. Here, we perform constant-potential molecular dynamics simulations and develop a lattice gas model to explore the reversible heat of EDL formation in aqueous and ionic liquid (IL) electrolytes. Our Letter reveals that EDL formation in aqueous electrolytes exhibits endothermicity under negative polarization; it shows new complexity of endothermicity followed by exothermicity in ILs, regardless of electrode polarity. These thermal behaviors are determined by the structural evolution during EDL formation, dominated by adsorbed solvent molecules rather than ions in aqueous electrolytes but governed by "demixing" and "vacancy occupation" phenomena in ILs. This Letter provides new insights into the reversible heat of supercapacitors and presents a theoretical approach to investigating thermal behaviors involving the dynamics of EDLs.
Collapse
Affiliation(s)
- Liang Zeng
- State Key Laboratory of Coal Combustion, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Ming Chen
- State Key Laboratory of Coal Combustion, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Zhenxiang Wang
- State Key Laboratory of Coal Combustion, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Rui Qiao
- Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Guang Feng
- State Key Laboratory of Coal Combustion, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| |
Collapse
|
28
|
Blazquez S, Abascal JLF, Lagerweij J, Habibi P, Dey P, Vlugt TJH, Moultos OA, Vega C. Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property. J Chem Theory Comput 2023; 19:5380-5393. [PMID: 37506381 PMCID: PMC10448725 DOI: 10.1021/acs.jctc.3c00562] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Indexed: 07/30/2023]
Abstract
In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.
Collapse
Affiliation(s)
- Samuel Blazquez
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L. F. Abascal
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jelle Lagerweij
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Parsa Habibi
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Poulumi Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Carlos Vega
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
29
|
Parisis N, Dans PD, Jbara M, Singh B, Schausi-Tiffoche D, Molina-Serrano D, Brun-Heath I, Hendrychová D, Maity SK, Buitrago D, Lema R, Nait Achour T, Giunta S, Girardot M, Talarek N, Rofidal V, Danezi K, Coudreuse D, Prioleau MN, Feil R, Orozco M, Brik A, Wu PYJ, Krasinska L, Fisher D. Histone H3 serine-57 is a CHK1 substrate whose phosphorylation affects DNA repair. Nat Commun 2023; 14:5104. [PMID: 37607906 PMCID: PMC10444856 DOI: 10.1038/s41467-023-40843-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2018] [Accepted: 08/12/2023] [Indexed: 08/24/2023] Open
Abstract
Histone post-translational modifications promote a chromatin environment that controls transcription, DNA replication and repair, but surprisingly few phosphorylations have been documented. We report the discovery of histone H3 serine-57 phosphorylation (H3S57ph) and show that it is implicated in different DNA repair pathways from fungi to vertebrates. We identified CHK1 as a major human H3S57 kinase, and disrupting or constitutively mimicking H3S57ph had opposing effects on rate of recovery from replication stress, 53BP1 chromatin binding, and dependency on RAD52. In fission yeast, mutation of all H3 alleles to S57A abrogated DNA repair by both non-homologous end-joining and homologous recombination, while cells with phospho-mimicking S57D alleles were partly compromised for both repair pathways, presented aberrant Rad52 foci and were strongly sensitised to replication stress. Mechanistically, H3S57ph loosens DNA-histone contacts, increasing nucleosome mobility, and interacts with H3K56. Our results suggest that dynamic phosphorylation of H3S57 is required for DNA repair and recovery from replication stress, opening avenues for investigating the role of this modification in other DNA-related processes.
Collapse
Affiliation(s)
- Nikolaos Parisis
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
- BPMP, CNRS, INRA, Montpellier SupAgro, University of Montpellier, Montpellier, France
- Institut Jacques Monod, CNRS, University Paris Diderot, Paris, France
| | - Pablo D Dans
- IRB Barcelona, BIST, Barcelona, Spain
- Bioinformatics Unit, Institute Pasteur of Montevideo, Montevideo, Uruguay
- Department of Biological Sciences, CENUR North Riverside, University of the Republic (UdelaR), Salto, Uruguay
| | - Muhammad Jbara
- Schulich Faculty of Chemistry, Technion Israel Institute of Technology, Haifa, Israel
- School of Chemistry, Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
| | | | | | | | | | - Denisa Hendrychová
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
- Department of Experimental Biology, Faculty of Science, Palacký University Olomouc, Olomouc, Czech Republic
| | - Suman Kumar Maity
- Schulich Faculty of Chemistry, Technion Israel Institute of Technology, Haifa, Israel
| | | | | | - Thiziri Nait Achour
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
| | - Simona Giunta
- The Rockefeller University, New York, NY, USA
- Laboratory of Genome Evolution, Department of Biology and Biotechnology "Charles Darwin", University of Rome Sapienza, Rome, Italy
| | - Michael Girardot
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
| | - Nicolas Talarek
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
| | - Valérie Rofidal
- BPMP, CNRS, INRA, Montpellier SupAgro, University of Montpellier, Montpellier, France
| | - Katerina Danezi
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
| | - Damien Coudreuse
- IGDR, CNRS, University of Rennes, Rennes, France
- IBGC, CNRS, University of Bordeaux, Bordeaux, France
| | | | - Robert Feil
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
| | | | - Ashraf Brik
- Schulich Faculty of Chemistry, Technion Israel Institute of Technology, Haifa, Israel
| | - Pei-Yun Jenny Wu
- IGDR, CNRS, University of Rennes, Rennes, France
- IBGC, CNRS, University of Bordeaux, Bordeaux, France
| | - Liliana Krasinska
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France.
- Equipe labellisée Ligue contre le Cancer, Paris, France.
| | - Daniel Fisher
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France.
- Equipe labellisée Ligue contre le Cancer, Paris, France.
| |
Collapse
|
30
|
Alexiou TS, Likos CN. Effective Interactions between Double-Stranded DNA Molecules in Aqueous Electrolyte Solutions: Effects of Molecular Architecture and Counterion Valency. J Phys Chem B 2023; 127:6969-6981. [PMID: 37493448 PMCID: PMC10424236 DOI: 10.1021/acs.jpcb.3c02216] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Revised: 06/12/2023] [Indexed: 07/27/2023]
Abstract
A computational investigation of the effects of molecular topology, namely, linear and circular, as well as counterion valency, on the ensuing pairwise effective interactions between DNA molecules in an unlinked state is presented. Umbrella sampling simulations have been performed through the introduction of bias potential along a reaction coordinate defined as the distance between the centers-of-mass of pairs of DNA molecules, and effective pair interaction potentials have been computed by employing the weighted histogram analysis method. An interesting comparison can be drawn between the different DNA topologies studied here, especially with regard to the contrasting effects of divalent counterions on the effective pair potentials: while DNA-DNA repulsion in short center-of-mass distances decreases significantly in the presence of divalent counterion-ions (as compared to monovalent ions) for linear DNA, the opposite effect occurs for the DNA minicircles. This can be attributed to the fact that linear DNA fragments can easily adopt relative orientations that minimize electrostatic and steric repulsions by rotating relative to one another and by exhibiting more pronounced bending due to the presence of free ends.
Collapse
Affiliation(s)
| | - Christos N Likos
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| |
Collapse
|
31
|
Saurabh S, Li Z, Hollowell P, Waigh T, Li P, Webster J, Seddon JM, Kalonia C, Lu JR, Bresme F. Structure and interaction of therapeutic proteins in solution: a combined simulation and experimental study. Mol Phys 2023; 121:e2236248. [PMID: 38107421 PMCID: PMC10721229 DOI: 10.1080/00268976.2023.2236248] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Accepted: 06/30/2023] [Indexed: 12/19/2023]
Abstract
The aggregation of therapeutic proteins in solution has attracted significant interest, driving efforts to understand the relationship between microscopic structural changes and protein-protein interactions determining aggregation processes in solution. Additionally, there is substantial interest in being able to predict aggregation based on protein structure as part of molecular developability assessments. Molecular Dynamics provides theoretical tools to complement experimental studies and to interrogate and identify the microscopic mechanisms determining aggregation. Here we perform all-atom MD simulations to study the structure and inter-protein interaction of the Fab and Fc fragments of the monoclonal antibody (mAb) COE3. We unravel the role of ion-protein interactions in building the ionic double layer and determining effective inter-protein interaction. Further, we demonstrate, using various state-of-the-art force fields (charmm, gromos, amber, opls/aa), that the protein solvation, ionic structure and protein-protein interaction depend significantly on the force field parameters. We perform SANS and Static Light Scattering experiments to assess the accuracy of the different forcefields. Comparison of the simulated and experimental results reveal significant differences in the forcefields' performance, particularly in their ability to predict the protein size in solution and inter-protein interactions quantified through the second virial coefficients. In addition, the performance of the forcefields is correlated with the protein hydration structure.
Collapse
Affiliation(s)
- Suman Saurabh
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London, United Kingdom
| | - Zongyi Li
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Peter Hollowell
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Thomas Waigh
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
- Photon Science Institute, The University of Manchester, Manchester, UK
| | - Peixun Li
- STFC ISIS Facility, Rutherford Appleton Laboratory, Didcot, UK
| | - John Webster
- STFC ISIS Facility, Rutherford Appleton Laboratory, Didcot, UK
| | - John M. Seddon
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London, United Kingdom
| | - Cavan Kalonia
- Dosage Form Design and Development, BioPharmaceutical Development, BioPharmaceuticals R&D, AstraZeneca, Gaithersburg, MD, USA
| | - Jian R. Lu
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London, United Kingdom
| |
Collapse
|
32
|
Chen WQ, Jivkov AP, Sedighi M. Thermo-Osmosis in Charged Nanochannels: Effects of Surface Charge and Ionic Strength. ACS APPLIED MATERIALS & INTERFACES 2023. [PMID: 37428544 PMCID: PMC10360061 DOI: 10.1021/acsami.3c02559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/11/2023]
Abstract
Thermo-osmosis refers to fluid migration due to the temperature gradient. The mechanistic understanding of thermo-osmosis in charged nano-porous media is still incomplete, while it is important for several environmental and energy applications, such as low-grade waste heat recovery, wastewater recovery, fuel cells, and nuclear waste storage. This paper presents results from a series of molecular dynamics simulations of thermo-osmosis in charged silica nanochannels that advance the understanding of the phenomenon. Simulations with pure water and water with dissolved NaCl are considered. First, the effect of surface charge on the sign and magnitude of the thermo-osmotic coefficient is quantified. This effect was found to be mainly linked to the structural modifications of an aqueous electrical double layer (EDL) caused by the nanoconfinement and surface charges. In addition, the results illustrate that the surface charges reduce the self-diffusivity and thermo-osmosis of interfacial liquid. The thermo-osmosis was found to change direction when the surface charge density exceeds -0.03C · m-2. It was found that the thermo-osmotic flow and self-diffusivity increase with the concentration of NaCl. The fluxes of solvent and solute are decoupled by considering the Ludwig-Soret effect of NaCl ions to identify the main mechanisms controlling the behavior. In addition to the advance in microscopic quantification and mechanistic understanding of thermo-osmosis, the work provides approaches to investigate a broader category of coupled heat and mass transfer problems in nanoscale space.
Collapse
Affiliation(s)
- Wei Qiang Chen
- School of Engineering, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - Andrey P Jivkov
- School of Engineering, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - Majid Sedighi
- School of Engineering, The University of Manchester, Manchester M13 9PL, United Kingdom
| |
Collapse
|
33
|
Tarrés-Solé A, Battistini F, Gerhold JM, Piétrement O, Martínez-García B, Ruiz-López E, Lyonnais S, Bernadó P, Roca J, Orozco M, Le Cam E, Sedman J, Solà M. Structural analysis of the Candida albicans mitochondrial DNA maintenance factor Gcf1p reveals a dynamic DNA-bridging mechanism. Nucleic Acids Res 2023; 51:5864-5882. [PMID: 37207342 PMCID: PMC10287934 DOI: 10.1093/nar/gkad397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Revised: 04/01/2023] [Accepted: 05/03/2023] [Indexed: 05/21/2023] Open
Abstract
The compaction of mitochondrial DNA (mtDNA) is regulated by architectural HMG-box proteins whose limited cross-species similarity suggests diverse underlying mechanisms. Viability of Candida albicans, a human antibiotic-resistant mucosal pathogen, is compromised by altering mtDNA regulators. Among them, there is the mtDNA maintenance factor Gcf1p, which differs in sequence and structure from its human and Saccharomyces cerevisiae counterparts, TFAM and Abf2p. Our crystallographic, biophysical, biochemical and computational analysis showed that Gcf1p forms dynamic protein/DNA multimers by a combined action of an N-terminal unstructured tail and a long helix. Furthermore, an HMG-box domain canonically binds the minor groove and dramatically bends the DNA while, unprecedentedly, a second HMG-box binds the major groove without imposing distortions. This architectural protein thus uses its multiple domains to bridge co-aligned DNA segments without altering the DNA topology, revealing a new mechanism of mtDNA condensation.
Collapse
Affiliation(s)
- Aleix Tarrés-Solé
- Structural MitoLab, Molecular Biology Institute Barcelona (IBMB-CSIC), Barcelona Science Park, Barcelona 08028, Spain
| | - Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona 08028, Spain
- Department of Biochemistry, University of Barcelona, Barcelona 08028, Spain
| | - Joachim M Gerhold
- Institute of Molecular and Cell Biology, University of Tartu, Tartu 51010, Estonia
| | - Olivier Piétrement
- Genome Integrity and Cancer UMR 9019 CNRS, Université Paris Saclay, Gustave Roussy Campus, 114 rue Edouard Vaillant 94805VillejuifCedex, France
| | | | - Elena Ruiz-López
- Structural MitoLab, Molecular Biology Institute Barcelona (IBMB-CSIC), Barcelona Science Park, Barcelona 08028, Spain
| | - Sébastien Lyonnais
- Structural MitoLab, Molecular Biology Institute Barcelona (IBMB-CSIC), Barcelona Science Park, Barcelona 08028, Spain
| | - Pau Bernadó
- Centre de Biologie Structurale (CBS), Inserm, CNRS and Université de Montpellier, France, Sébastien Lyonnais, UAR 3725 CNRS, Université de Montpellier, 34000 Montpellier, France
| | - Joaquim Roca
- Molecular Biology Institute Barcelona (IBMB-CSIC), Barcelona 08028, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona 08028, Spain
- Department of Biochemistry, University of Barcelona, Barcelona 08028, Spain
| | - Eric Le Cam
- Genome Integrity and Cancer UMR 9019 CNRS, Université Paris Saclay, Gustave Roussy Campus, 114 rue Edouard Vaillant 94805VillejuifCedex, France
| | - Juhan Sedman
- Institute of Molecular and Cell Biology, University of Tartu, Tartu 51010, Estonia
| | - Maria Solà
- Structural MitoLab, Molecular Biology Institute Barcelona (IBMB-CSIC), Barcelona Science Park, Barcelona 08028, Spain
| |
Collapse
|
34
|
Ma L, Zhong Z, Hu J, Qing L, Jiang J. Long-Lived Weak Ion Pairs in Ionic Liquids: An Insight from All-Atom Molecular Dynamics Simulations. J Phys Chem B 2023. [PMID: 37262343 DOI: 10.1021/acs.jpcb.3c01559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
The microstructure and local dynamics of ions in room-temperature ionic liquids (RTILs) have drawn a lot of attention due to their extensive potential applications in numerous fields. It is well-known that the widely used definitions of ion pairs (IPs) cannot reflect the full picture of RTILs. In this study, we find a universal residence time (τMR), which is regardless of the number of counterions in the first solvation shell in RTILs. Inspired by this, we propose a weak IP (WIP) model from a spatiotemporal perspective and demonstrate that the WIPs are long-lived and that their lifetimes obey a log-normal distribution, which is different from the literature. In addition, the electrostatic interactions are the main factors in the formation of WIPs, and the reorientations of ions are vital to the ruptures of WIPs. This research provides a new perspective for understanding the microstructural and dynamical properties of RTILs.
Collapse
Affiliation(s)
- Linbo Ma
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Zhixuan Zhong
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Junbao Hu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Leying Qing
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
| | - Jian Jiang
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| |
Collapse
|
35
|
Wei Z, Xu H, Xu X, Feng G, Zheng W, Li T. Solvation of magnesium chloride dimer in water: The case of anionic and neutral clusters. J Chem Phys 2023; 158:2888211. [PMID: 37140000 DOI: 10.1063/5.0146319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2023] [Accepted: 04/24/2023] [Indexed: 05/05/2023] Open
Abstract
The structures of magnesium chloride dimer-water clusters, (MgCl2)2(H2O)n-/0, were investigated with size-selected anion photoelectron spectroscopy and theoretical calculations to understand the dissolution of magnesium chloride in water. The most stable structures were confirmed by comparing vertical detachment energies (VDEs) with the experimental measurements. A dramatic drop of VDE at n = 3 has been observed in the experiment, which is in accordance with the structural change of (MgCl2)2(H2O)n-. Compared to the neutral clusters, the excess electron induces two significant phenomena in (MgCl2)2(H2O)n-. First, the planar D2h geometry can be converted into a C3v structure at n = 0, making the Mg-Cl bonds easier to be broken by water molecules. More importantly, a negative charge-transfer-to-solvent process occurs after adding three water molecules (i.e., at n = 3), which leads to an obvious deviation in the evolution of the clusters. Such electron transfer behavior was noticed at n = 1 in monomer MgCl2(H2O)n-, indicating that the dimerization between two MgCl2 molecules can make the cluster more capable of binding electron. In neutral (MgCl2)2(H2O)n, this dimerization provides more sites for the added water molecules, which can stabilize the entire cluster and maintain its initial structure. Specifically, filling the coordination number to be 6 for Mg atoms can be seen as a link between structural preferences in the dissolution of the monomers, dimers, and extended bulk-state of MgCl2. This work represents an important step forward into fully understanding the solvation of MgCl2 crystals and other multivalent salt oligomers.
Collapse
Affiliation(s)
- Zhiyou Wei
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325001, China
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hongguang Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xiling Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Gang Feng
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China
| | - Weijun Zheng
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Tao Li
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325001, China
- Oujiang Laboratory (Zhejiang Lab for Regenerative Medicine, Vision and Brain Health), Wenzhou, Zhejiang 325001, China
| |
Collapse
|
36
|
Abbaspour M. Polyoxometalate ionic liquid between graphene oxide surfaces as a new membrane in the desalination process: a molecular dynamics study. Phys Chem Chem Phys 2023; 25:13654-13664. [PMID: 37145119 DOI: 10.1039/d2cp05486h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
In this study, the performance of the positioning of polyoxometalate ionic liquid ([Keggin][emim]3 IL) between graphene oxide (GO) plates with different concentrations (nIL-GO (n = 1-4)) were examined in the desalination process at different external pressures using molecular dynamics (MD) simulations. The use of Keggin anions with charged GO layers was also investigated in the desalination process. The potential of the mean force, average number of hydrogen bonds, self-diffusion coefficient, and angle distribution function were calculated and discussed. The results showed that although the presence of polyoxometalate ILs between the GO plates decreases water flux, they efficiently increase salt rejection. The positioning of one IL increases salt rejection to two times at lower pressure and increases it up to four times at higher pressure. Moreover, the positioning of four ILs results in almost complete salt rejection at all pressures. The use of only Keggin anions between the charged GO plates (n[Keggin]-GO+3n) presents more water flux and a smaller salt rejection rate than the nIL-GO systems. However, the n[Keggin]-GO+3n systems show a nearly complete salt rejection at high concentrations of Keggin anions. These systems also have a smaller risk of the contamination of the desalinated water by the probable escape of cations from the nanostructure to the desalinated water at very high pressures.
Collapse
Affiliation(s)
- Mohsen Abbaspour
- Department of Chemistry, Hakim Sabzevari University, Sabzevar, Iran.
| |
Collapse
|
37
|
Lemay AC, Sontarp EJ, Martinez D, Maruri P, Mohammed R, Neapole R, Wiese M, Willemsen JAR, Bourg IC. Molecular Dynamics Simulation Prediction of the Partitioning Constants ( KH, Kiw, Kia) of 82 Legacy and Emerging Organic Contaminants at the Water-Air Interface. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2023; 57:6296-6308. [PMID: 37014786 DOI: 10.1021/acs.est.3c00267] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
The tendency of organic contaminants (OCs) to partition between different phases is a key set of properties that underlie their human and ecological health impacts and the success of remediation efforts. A significant challenge associated with these efforts is the need for accurate partitioning data for an ever-expanding list of OCs and breakdown products. All-atom molecular dynamics (MD) simulations have the potential to help generate these data, but existing studies have applied these techniques only to a limited variety of OCs. Here, we use established MD simulation approaches to examine the partitioning of 82 OCs, including many compounds of critical concern, at the water-air interface. Our predictions of the Henry's law constant (KH) and interfacial adsorption coefficients (Kiw, Kia) correlate strongly with experimental results, indicating that MD simulations can be used to predict KH, Kiw, and Kia values with mean absolute deviations of 1.1, 0.3, and 0.3 logarithmic units after correcting for systematic bias, respectively. A library of MD simulation input files for the examined OCs is provided to facilitate future investigations of the partitioning of these compounds in the presence of other phases.
Collapse
Affiliation(s)
- Amélie C Lemay
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ethan J Sontarp
- Department of Geosciences, Princeton University, Princeton, New Jersey 08544, United States
| | - Daniela Martinez
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Philip Maruri
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Raneem Mohammed
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ryan Neapole
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Morgan Wiese
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Jennifer A R Willemsen
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ian C Bourg
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
- High Meadows Environmental Institute, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
38
|
Abbaspour M, Fotourechi F, Akbarzadeh H, Salemi S. Investigation of small inhibitor effects on methane hydrate formation in a carbon nanotube using molecular dynamics simulation. RSC Adv 2023; 13:6800-6807. [PMID: 36865572 PMCID: PMC9971844 DOI: 10.1039/d2ra06518e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2022] [Accepted: 02/16/2023] [Indexed: 03/04/2023] Open
Abstract
In this work, we simulated water molecules in fixed and rigid (15,0) CNTs and the confined water molecules formed a hexagonal ice nanotube in the CNT. After the addition of methane molecules in the nanotube, the hexagonal structure of confined water molecules disappeared and were replaced by almost all the guest methane molecules. The replaced molecules formed a row of water molecules in the middle of the hollow space of the CNT. We also added five small inhibitors with different concentrations (0.8 mol% and 3.8 mol%) to methane clathrates in CNT: benzene, 1-ethyl-3-methylimidazolium chloride ionic liquid ([emim+][Cl-] IL), methanol, NaCl, and tetrahydrofuran (THF). We investigated the thermodynamic and kinetic inhibition behaviors of the different inhibitors on the methane clathrate formation in the CNT using the radial distribution function (RDF), hydrogen bonding (HB), and angle distribution function (ADF). Our results showed that the [emim+][Cl-] IL is the best inhibitor from both aspects. It was also shown that the effect of THF and benzene is better than that of NaCl and methanol. Furthermore, our results showed that the THF inhibitors tended to aggregate in the CNT, but the benzene and IL molecules were distributed along the CNT and can affect the inhibitor behavior of THF in the CNT. We have also examined the effect of CNT chirality using the armchair (9,9) CNT, the effect of CNT size using the (17,0) CNT, and the effect of CNT flexibility using the (15,0) CNT by the DREIDING force field. Our results showed that the IL has stronger thermodynamic and kinetic inhibition effects in the armchair (9,9) and the flexible (15,0) CNT than the other systems, respectively.
Collapse
Affiliation(s)
| | | | - Hamed Akbarzadeh
- Dep. of Physical Chemistry, Faculty of Chemistry, Kharazmi UniversityTehranIran
| | - Sirous Salemi
- Dep. of Chemistry, Hakim Sabzevari University Sabzevar Iran
| |
Collapse
|
39
|
Matthes D, de Groot BL. Molecular dynamics simulations reveal the importance of amyloid-beta oligomer β-sheet edge conformations in membrane permeabilization. J Biol Chem 2023; 299:103034. [PMID: 36806684 PMCID: PMC10033322 DOI: 10.1016/j.jbc.2023.103034] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 02/03/2023] [Accepted: 02/13/2023] [Indexed: 02/18/2023] Open
Abstract
Oligomeric aggregates of the amyloid-beta peptide(1-42) (Aβ42) are regarded as a primary cause of cytotoxicity related to membrane damage in Alzheimer's disease. However, a dynamical and structural characterization of pore-forming Aβ42 oligomers at atomic detail has not been feasible. Here, we used Aβ42 oligomer structures previously determined in a membrane-mimicking environment as putative model systems to study the pore formation process in phospholipid bilayers with all-atom molecular dynamics simulations. Multiple Aβ42 oligomer sizes, conformations, and N-terminally truncated isoforms were investigated on the multi-μs time scale. We found that pore formation and ion permeation occur via edge conductivity and exclusively for β-sandwich structures that feature exposed side-by-side β-strand pairs formed by residues 9 to 21 of Aβ42. The extent of pore formation and ion permeation depends on the insertion depth of hydrophilic residues 13 to 16 (HHQK domain) and thus on subtle differences in the overall stability, orientation, and conformation of the aggregates in the membrane. Additionally, we determined that backbone carbonyl and polar side-chain atoms from the edge strands directly contribute to the coordination sphere of the permeating ions. Furthermore, point mutations that alter the number of favorable side-chain contacts correlate with the ability of the Aβ42 oligomer models to facilitate ion permeation in the bilayer center. Our findings suggest that membrane-inserted, layered β-sheet edges are a key structural motif in pore-forming Aβ42 oligomers independent of their size and play a pivotal role in aggregate-induced membrane permeabilization.
Collapse
Affiliation(s)
- Dirk Matthes
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany.
| | - Bert L de Groot
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany
| |
Collapse
|
40
|
A Fijivirus Major Viroplasm Protein Shows RNA-Stimulated ATPase Activity by Adopting Pentameric and Hexameric Assemblies of Dimers. mBio 2023; 14:e0002323. [PMID: 36786587 PMCID: PMC10128069 DOI: 10.1128/mbio.00023-23] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023] Open
Abstract
Fijiviruses replicate and package their genomes within viroplasms in a process involving RNA-RNA and RNA-protein interactions. Here, we demonstrate that the 24 C-terminal residues (C-arm) of the P9-1 major viroplasm protein of the mal de Río Cuarto virus (MRCV) are required for its multimerization and the formation of viroplasm-like structures. Using an integrative structural approach, the C-arm was found to be dispensable for P9-1 dimer assembly but essential for the formation of pentamers and hexamers of dimers (decamers and dodecamers), which favored RNA binding. Although both P9-1 and P9-1ΔC-arm catalyzed ATP with similar activities, an RNA-stimulated ATPase activity was only detected in the full-length protein, indicating a C-arm-mediated interaction between the ATP catalytic site and the allosteric RNA binding sites in the (do)decameric assemblies. A stronger preference to bind phosphate moieties in the decamer was predicted, suggesting that the allosteric modulation of ATPase activity by RNA is favored in this structural conformation. Our work reveals the structural versatility of a fijivirus major viroplasm protein and provides clues to its mechanism of action. IMPORTANCE The mal de Río Cuarto virus (MRCV) causes an important maize disease in Argentina. MRCV replicates in several species of Gramineae plants and planthopper vectors. The viral factories, also called viroplasms, have been studied in detail in animal reovirids. This work reveals that a major viroplasm protein of MRCV forms previously unidentified structural arrangements and provides evidence that it may simultaneously adopt two distinct quaternary assemblies. Furthermore, our work uncovers an allosteric communication between the ATP and RNA binding sites that is favored in the multimeric arrangements. Our results contribute to the understanding of plant reovirids viroplasm structure and function and pave the way for the design of antiviral strategies for disease control.
Collapse
|
41
|
Blazquez S, Conde MM, Vega C. Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water. J Chem Phys 2023; 158:054505. [PMID: 36754806 DOI: 10.1063/5.0136498] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
42
|
Building a Hofmeister-like series for the maximum in density temperature of aqueous electrolyte solutions. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
|
43
|
Tong J, Peng B, Kontogeorgis GM, Liang X. Behavior of the aqueous sodium chloride solutions from molecular simulations and theories. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2022.121086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
44
|
Yao S, Van R, Pan X, Park JH, Mao Y, Pu J, Mei Y, Shao Y. Machine learning based implicit solvent model for aqueous-solution alanine dipeptide molecular dynamics simulations. RSC Adv 2023; 13:4565-4577. [PMID: 36760282 PMCID: PMC9900604 DOI: 10.1039/d2ra08180f] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 01/20/2023] [Indexed: 02/05/2023] Open
Abstract
Inspired by the recent work from Noé and coworkers on the development of machine learning based implicit solvent model for the simulation of solvated peptides [Chen et al., J. Chem. Phys., 2021, 155, 084101], here we report another investigation of the possibility of using machine learning (ML) techniques to "derive" an implicit solvent model directly from explicit solvent molecular dynamics (MD) simulations. For alanine dipeptide, a machine learning potential (MLP) based on the DeepPot-SE representation of the molecule was trained to capture its interactions with its average solvent environment configuration (ASEC). The predicted forces on the solute deviated only by an RMSD of 0.4 kcal mol-1 Å-1 from the reference values, and the MLP-based free energy surface differed from that obtained from explicit solvent MD simulations by an RMSD of less than 0.9 kcal mol-1. Our MLP training protocol could also accurately reproduce combined quantum mechanical molecular mechanical (QM/MM) forces on the quantum mechanical (QM) solute in ASEC environment, thus enabling the development of accurate ML-based implicit solvent models for ab initio-QM MD simulations. Such ML-based implicit solvent models for QM calculations are cost-effective in both the training stage, where the use of ASEC reduces the number of data points to be labelled, and the inference stage, where the MLP can be evaluated at a relatively small additional cost on top of the QM calculation of the solute.
Collapse
Affiliation(s)
- Songyuan Yao
- Department of Chemistry and Biochemistry, University of Oklahoma Norman OK 73019 USA
| | - Richard Van
- Department of Chemistry and Biochemistry, University of Oklahoma Norman OK 73019 USA
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma Norman OK 73019 USA
| | - Ji Hwan Park
- School of Computer Science, University of Oklahoma Norman OK 73019 USA
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University San Diego CA 92182 USA
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis Indianapolis IN 46202 USA
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200062 China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai Shanghai 200062 China
- Collaborative Innovation Center of Extreme Optics, Shanxi University Taiyuan Shanxi 030006 China
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma Norman OK 73019 USA
| |
Collapse
|
45
|
Chen B, Li K, Sun H, Jiang L, Yang M, Song Y. Promoting Effect of Common Marine Cations on Hydrate Dissociation and Structural Evolution under a Static Electric Field. J Phys Chem B 2023; 127:698-709. [PMID: 36629388 DOI: 10.1021/acs.jpcb.2c05382] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Natural gas hydrate, a potential energy resource, is attracting worldwide attention. In this study, we propose a new method of hydrate dissociation which uses seawater and electrostatic fields (SE method) cooperatively. The hydrate molecular dissociation mechanism of gas hydrate is a key issue in studying the kinetic properties of gas hydrate using the SE method. Therefore, molecular dynamics simulations were used to investigate the thermodynamic properties and structural changes of methane hydrate (MH) in multiple kinds of salt solutions under an electrostatic field. The results show that the electric field can drive cations into the MH phase to form a series of random semiopen cages, which are essentially temporary and metastable. The variation in free energy indicates that it is more difficult for divalent cations to enter the hydrate phase than monovalent cations, meaning that the hydrate structures formed with divalent cations are more unstable. Then, the ion current occurred in the hydrate phase (called ion migration in this study), which greatly accelerated hydrate dissociation. In contrast, the promotion effect of cations with the same charge on MH dissociation is as follows: Sr2+ > K+ ≈ Na+ > Ca2+ ≈ Mg2+. In general, the presence of common marine cations enhanced the promotion effect of the electric field on gas hydrate dissociation.
Collapse
Affiliation(s)
- Bingbing Chen
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian116024, China
| | - Kehan Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian116024, China
| | - Huiru Sun
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian116024, China
| | - Lanlan Jiang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian116024, China
| | - Mingjun Yang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian116024, China
| | - Yongchen Song
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian116024, China
| |
Collapse
|
46
|
Park SJ, Narvaez WA, Schwartz BJ. Ab Initio Studies of Hydrated Electron/Cation Contact Pairs: Hydrated Electrons Simulated with Density Functional Theory Are Too Kosmotropic. J Phys Chem Lett 2023; 14:559-566. [PMID: 36630724 DOI: 10.1021/acs.jpclett.2c03705] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
We have performed the first DFT-based ab initio MD simulations of a hydrated electron (eaq-) in the presence of Na+, a system chosen because ion-pairing behavior in water depends sensitively on the local hydration structure. Experiments show that eaq-'s interact weakly with Na+; the eaq-'s spectrum blue shifts by only a few tens of meV upon ion pairing without changing shape. We find that the spectrum of the DFT-simulated eaq- red shifts and changes shape upon interaction with Na+, in contrast with experiment. We show that this is because the hydration structure of the DFT-simulated eaq- is too ordered or kosmotropic. Conversely, simulations that produce eaq-'s with a less ordered or chaotropic hydration structure form weaker ion pairs with Na+, yielding predicted spectral blue shifts in better agreement with experiment. Thus, ab initio simulations based on hybrid GGA DFT functionals fail to produce the correct solvation structure for the hydrated electron.
Collapse
Affiliation(s)
- Sanghyun J Park
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Wilberth A Narvaez
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Benjamin J Schwartz
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| |
Collapse
|
47
|
Chen WQ, Sedighi M, Jivkov AP. Thermal diffusion of ionic species in charged nanochannels. NANOSCALE 2022; 15:215-229. [PMID: 36468769 DOI: 10.1039/d2nr05504j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Diffusion of ions due to temperature gradients (known as thermal diffusion) in charged nanochannels is of interest in several engineering fields, including energy recovery and environmental protection. This paper presents a fundamental investigation of the thermal diffusion of sodium chloride in charged silica nanochannels performed by molecular dynamics (MD). The results reveal the effects of nanoconfinement and surface charges on the sign and magnitude of the Soret coefficient. It is shown that the sign and magnitude of the Soret coefficient are controlled by the structural modifications of the interfacial solutions. These modifications include the ionic solvation and hydrogen bond structure induced by the nanoconfinement and surface charges. The results show that both nanoconfinement and surface charges can make the solutions more thermophilic. Furthermore, the thermal diffusion of solutions in boundary layers is significantly different from that of solutions in bulk fluid, contributing to the overall difference between the thermal diffusivity of pore fluid and that associated with bulk fluid. The findings provide further understanding of thermal diffusion in nano-porous systems. The proposed MD simulation methodology is applicable to a wider category of coupled heat and mass transfer problems in nanoscale spaces.
Collapse
Affiliation(s)
- Wei Qiang Chen
- Department of Mechanical, Aerospace and Civil Engineering, School of Engineering, The University of Manchester, Manchester, M13 9PL, UK.
| | - Majid Sedighi
- Department of Mechanical, Aerospace and Civil Engineering, School of Engineering, The University of Manchester, Manchester, M13 9PL, UK.
| | - Andrey P Jivkov
- Department of Mechanical, Aerospace and Civil Engineering, School of Engineering, The University of Manchester, Manchester, M13 9PL, UK.
| |
Collapse
|
48
|
Shen X, Bourg IC. Interaction between Hydrated Smectite Clay Particles as a Function of Salinity (0-1 M) and Counterion Type (Na, K, Ca). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:20990-20997. [PMID: 37881773 PMCID: PMC10595998 DOI: 10.1021/acs.jpcc.2c04636] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 11/11/2022] [Indexed: 10/27/2023]
Abstract
Swelling clay minerals control the hydrologic and mechanical properties of many soils, sediments, and sedimentary rocks. This important and well-known phenomenon remains challenging to predict because it emerges from complex multiscale couplings between aqueous chemistry and colloidal interaction mechanics in nanoporous clay assemblages, for which predictive models remain elusive. In particular, the predominant theory of colloidal interactions across fluid films, the widely used Derjaguin-Landau-Verwey-Overbeek model, fails to predict the ubiquitous existence of stable swelling states at interparticle distances below 3 nm that are stabilized by specific inter-atomic interactions in overlapping electrical double layers between the charged clay surfaces. Atomistic simulations have the potential to generate detailed insights into the mechanisms of these interactions. Recently, we developed a metadynamics-based molecular dynamics simulation methodology that can predict the free energy of interaction between parallel smectite clay particles in a wide range of interparticle distances (from 0.3 to 3 nm) and salinities (from 0.0 to 1.0 M NaCl). Here, we extend this work by characterizing the sensitivity of interparticle interactions to counterion type (Na, K, Ca). We establish a detailed picture of the free energy of interaction of parallel clay particles across water films as the sum of five interaction mechanisms with different sensitivities to salinity, counterion type, and interparticle distance.
Collapse
Affiliation(s)
- Xinyi Shen
- Department of Civil and Environmental
Engineering and High Meadows Environmental Institute, Princeton University, Princeton, New Jersey08544, United States
| | - Ian C. Bourg
- Department of Civil and Environmental
Engineering and High Meadows Environmental Institute, Princeton University, Princeton, New Jersey08544, United States
| |
Collapse
|
49
|
Song Y, Li K, Sun H, Chen B, Yang M. New Sights on derived behaviors of methane hydrate molecular structure in Na+/Cl- ions invading process. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
|
50
|
Chen YJ, Yu HY. Enthalpic Interactions and Solution Behaviors of Solvent-Free Polymer Brushes. Polymers (Basel) 2022; 14:polym14235237. [PMID: 36501627 PMCID: PMC9740690 DOI: 10.3390/polym14235237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 11/25/2022] [Accepted: 11/28/2022] [Indexed: 12/04/2022] Open
Abstract
We performed molecular dynamics simulations to characterize the role of enthalpic interaction in impacting the static and dynamic properties of solvent-free polymer brushes. The intrinsic enthalpic interaction in the simulation was introduced using different attraction strengths between distinct species. Two model systems were considered: one consisting of binary brushes of two different polymer types and the other containing a mixture of homopolymer brushes and free molecules. In the first system, we observed that, when two originally incompatible polymers were grafted to opposing surfaces, the miscibility between them was significantly enhanced. A less favorable intrinsic enthalpic interaction in the brushes resulted in a more stretched chain configuration, a lower degree of inter-brush penetration, and faster segmental relaxation. In the second system, we characterized the solvent capacity of the homopolymer brushes from variations in the energy components of the system as a function of the number of free molecules. We determined that molecular absorption was driven by the release of the entropic frustration for the grafted chains in conjunction with the chemical affinity between the solutes and polymers. The solute distribution function within the inter-wall space showed that solute-polymer mixing in the middle of the gap occurred preferentially when the enthalpic interaction was more favorable. When this was not the case, absorption was predominantly localized near the grafting surface. From the mean square displacement of the solute, we found that the brush profiles restrained the molecular diffusion perpendicular to the grafting wall; the weaker the attraction from the brush, the higher the solute mobility.
Collapse
|