Computer simulations of NaCl association in polarizable water

The impact of sequence periodicity on DNA mechanics: investigating the origin of A-tract's curvature

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.

Force field for halide and alkali ions in water based on single-ion and ion-pair thermodynamic properties for a wide range of concentrations

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.

Diffusion and structure of propylene carbonate-metal salt electrolyte solutions for post-lithium-ion batteries: From experiment to simulation

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.

CO2 ultrathin film growth on a monolayer of CO2 adsorbed on the NaCl(100) surface: sticking coefficient and IR-optical signatures in the ν3 region

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.

Electrokinetic properties of NaCl solution via molecular dynamics simulations with scaled-charge electrolytes

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.

Molecular dynamics simulations of choline chloride and ascorbic acid deep eutectic solvents: Investigation of structural and dynamics properties

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.

To Pair or not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water

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.

All-atom molecular dynamics simulations of polymer and polyelectrolyte brushes

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.

Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores

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.

Evidence of Sharp Transitions between Octahedral and Capped Trigonal Prism States of the Solvation Shell of the Aqueous Fe3+ Ion

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.

Vehicular Motions Dominate the Ion Transport in Concentrated LiTFSI Aqueous Solutions?

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.

Calculations of Absolute Free Energies, Enthalpies, and Entropies for Drug Binding

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.

Interfacial properties of the brine + carbon dioxide + oil + silica system

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.

Ion diffusion captures composition-dependent anomalies in water-DMSO binary mixtures

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.

Machine Learning-Based Interfacial Tension Equations for (H2 + CO2)-Water/Brine Systems over a Wide Range of Temperature and Pressure

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.

Energy-Scaled Debye-Hückel Theory for the Electrostatic Solvation Free Energy in Size-Asymmetric Electrolyte Solutions

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.

Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2

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.

Hygroscopic Growth of Adsorbed Water Films on Smectite Clay Particles

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.

Facilitation of Hydrate Dissociation and Structural Evolution by Major Marine Anions under Static Electric Fields

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.

Solubility Equilibrium Isotope Effects of Noble Gases in Water: Theory and Observations

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.

Novel Kinetic Resolution of Thiazolo-Benzimidazolines Using MAO Enzymes

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.

Scaling Solute-Solvent Distances to Improve Solubility and Ion Paring Predictions in Rigid Ion Models

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.

Transport Properties of NaCl in Aqueous Solution and Hydrogen Solubility in Brine

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.

Benchmark computational investigations for the basic model of the salt-water complex: NaCl(H2O) and its anion NaCl(H2O). Phys Chem Chem Phys 2023;25:27215-27229. [PMID: 37791409 DOI: 10.1039/d3cp03421f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

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.

Nanoscale Prediction of the Thermal, Mechanical, and Transport Properties of Hydrated Clay on 106- and 1015-Fold Larger Length and Time Scales

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.

Interfacial Tension-Temperature-Pressure-Salinity Relationship for the Hydrogen-Brine System under Reservoir Conditions: Integration of Molecular Dynamics and Machine Learning

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.

Structural Evolution Governs Reversible Heat Generation in Electrical Double Layers

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.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

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.

Histone H3 serine-57 is a CHK1 substrate whose phosphorylation affects DNA repair

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.

Effective Interactions between Double-Stranded DNA Molecules in Aqueous Electrolyte Solutions: Effects of Molecular Architecture and Counterion Valency

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.

Structure and interaction of therapeutic proteins in solution: a combined simulation and experimental study

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.

Thermo-Osmosis in Charged Nanochannels: Effects of Surface Charge and Ionic Strength

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.

Structural analysis of the Candida albicans mitochondrial DNA maintenance factor Gcf1p reveals a dynamic DNA-bridging mechanism

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.

Long-Lived Weak Ion Pairs in Ionic Liquids: An Insight from All-Atom Molecular Dynamics Simulations

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.

Solvation of magnesium chloride dimer in water: The case of anionic and neutral clusters

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.

Polyoxometalate ionic liquid between graphene oxide surfaces as a new membrane in the desalination process: a molecular dynamics study

In this study, the performance of the positioning of polyoxometalate ionic liquid ([Keggin][emim]₃ 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⁺³ⁿ) presents more water flux and a smaller salt rejection rate than the nIL-GO systems. However, the n[Keggin]-GO⁺³ⁿ 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.

Molecular Dynamics Simulation Prediction of the Partitioning Constants (KH, Kiw, Kia) of 82 Legacy and Emerging Organic Contaminants at the Water-Air Interface

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 (K_H) and interfacial adsorption coefficients (K_iw, K_ia) correlate strongly with experimental results, indicating that MD simulations can be used to predict K_H, K_iw, and K_ia 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.

Investigation of small inhibitor effects on methane hydrate formation in a carbon nanotube using molecular dynamics simulation

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.

Molecular dynamics simulations reveal the importance of amyloid-beta oligomer β-sheet edge conformations in membrane permeabilization

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.

A Fijivirus Major Viroplasm Protein Shows RNA-Stimulated ATPase Activity by Adopting Pentameric and Hexameric Assemblies of Dimers

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.

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water

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.

Machine learning based implicit solvent model for aqueous-solution alanine dipeptide molecular dynamics simulations

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.

Promoting Effect of Common Marine Cations on Hydrate Dissociation and Structural Evolution under a Static Electric Field

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: Sr²⁺ > K⁺ ≈ Na⁺ > Ca²⁺ ≈ Mg²⁺. In general, the presence of common marine cations enhanced the promotion effect of the electric field on gas hydrate dissociation.

Ab Initio Studies of Hydrated Electron/Cation Contact Pairs: Hydrated Electrons Simulated with Density Functional Theory Are Too Kosmotropic

We have performed the first DFT-based ab initio MD simulations of a hydrated electron (e_aq^-) in the presence of Na⁺, a system chosen because ion-pairing behavior in water depends sensitively on the local hydration structure. Experiments show that e_aq^-'s interact weakly with Na⁺; the e_aq^-'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 e_aq^- 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 e_aq^- is too ordered or kosmotropic. Conversely, simulations that produce e_aq^-'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.

Thermal diffusion of ionic species in charged nanochannels

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.

Interaction between Hydrated Smectite Clay Particles as a Function of Salinity (0-1 M) and Counterion Type (Na, K, Ca)

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.

Enthalpic Interactions and Solution Behaviors of Solvent-Free Polymer Brushes

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.