The multi-scale polarizable pseudo-particle solvent coarse-grained approach: From NaCl salt solutions to polyelectrolyte hydration.

We discuss key parameters that affect the reliability of hybrid simulations in the aqueous phase based on an efficient multi-scale coarse-grained polarizable pseudo-particle approach, denoted as pppl, to model the solvent water, whereas solutes are modeled using an all atom polarizable force field. Among those parameters, the extension of the solvent domain (SD) at the solute vicinity (domain in which each solvent particle corresponds to a single water molecule) and the magnitude of solute/solvent short range polarization damping effects are shown to be pivotal to model NaCl salty aqueous solutions and the hydration of charged systems, such as the hydrophobic polyelectrolyte polymer that we have recently investigated [Masella et al., J. Chem. Phys. 155, 114903 (2021)]. Strong short range damping is pivotal to simulate aqueous salt NaCl solutions at moderate concentration (up to 1.0M). The SD extension (as well as short range damping) has a weak effect on the polymer conformation; however, it plays a pivotal role in computing accurate polymer/solvent interaction energies. As the pppl approach is up to two orders of magnitude computationally more efficient than all atom polarizable force field methods, our results show it to be an efficient alternative route to investigate the equilibrium properties of complex charged molecular systems in extended chemical environments.

Early Steps in the O2 Scavenger Process in the Aqueous Phase: Hydrazine vs DEHA.

We investigate the first direct proton abstraction reactions from reducing agents (RAHs) hydrazine and diethyl hydroxylamine (DEHA), toward dioxygen (O2) in the aqueous phase, spanning ambient to high-temperature conditions. Quantum chemistry methods and molecular dynamics simulations are employed in this study. Quantum chemistry methods are used to analyze the quasi-equilibrium between a reactive conformation and a transition state in the [RAH,O2] cluster. On the other hand, molecular dynamics simulations estimate the probability of observing a reactive conformation of the [RAH,O2] cluster in the solution. In this study, we assume that the energy barrier of the quasi-equilibrium is sufficiently high for the RAH/O2 association process to be at equilibrium. Our findings indicate that the first proton abstraction process from a reactive conformation cluster by DEHA is energetically favored compared to hydrazine. Conversely, the association process of hydrazine and O2 in solution is more favorable than that of DEHA. Consequently, the rate constant for the first proton abstraction process is similar for both hydrazine and DEHA, particularly at high temperatures, with activation energies of approximately 21.5 ± 1.5 kcal mol-1 for both compounds. These results align with recent experiments investigating the complete O2 scavenger process in liquid water with hydrazine and DEHA. Therefore, our findings support the assumption that first proton abstraction reactions are the rate-determining steps in O2 scavenger processes in the aqueous phase.

Chitosan Polysaccharides from a Polarizable Multiscale Approach.

We report simulations of chitosan polysaccharides in the aqueous phase, at infinite dilute conditions and zero ionic strength. Those simulations are performed by means of a polarizable multiscale modeling scheme that relies on a polarizable all atom force field to model solutes and on a polarizable solvent coarse grained approach. Force field parameters are assigned only from quantum chemistry ab initio data. We simulate chitosan monomer units, dimers and 50-long chains. Regarding the 50-long chains we simulate three sets of ten randomly built chain replica at three different pH conditions (corresponding to different chain protonation states, the chain degree of deacetylation is 85%). Our simulations show the persistence length of 50-long chitosan chains at strong acidic conditions (pH <5) to be 24 ± 2 nm (at weak/negligible ionic strength conditions), and to be 1 order of magnitude shorter at usual pH conditions. Our simulation data support the most recent simulation and experimental studies devoted to chitosan polysaccharides in the aqueous phase.

Protein-ligand interactions from a quantum fragmentation perspective: The case of the SARS-CoV-2 main protease interacting with α-ketoamide inhibitors.

We present a hybrid, multi-method, computational scheme for protein/ligand systems well suited to be used on modern and forthcoming massively parallel computing systems. The scheme relies on a multi-scale polarizable molecular modeling, approach to perform molecular dynamics simulations, and on an efficient Density Functional Theory (DFT) linear scaling method to post-process simulation snapshots. We use this scheme to investigate recent α-ketoamide inhibitors targeting the main protease of the SARS-CoV-2 virus. We assessed the reliability and the coherence of the hybrid scheme, in particular, by checking the ability of MM and DFT to reproduce results from high-end ab initio computations regarding such inhibitors. The DFT approach enables an a posteriori fragmentation of the system and an investigation into the strength of interaction among identified fragment pairs. We show the necessity of accounting for a large set of plausible protease/inhibitor conformations to generate reliable interaction data. Finally, we point out ways to further improve α-ketoamide inhibitors to more strongly interact with particular protease domains neighboring the active site.

Hybrid polarizable simulations of a conventional hydrophobic polyelectrolyte. Toward a theoretical tool for green science innovation.

Hybrid modeling approaches based on all-atom force fields to handle a solute and coarse-grained models to account for the solvent are promising numerical tools that can be used to understand the properties of large and multi-components solutions and thus to speed up the development of new industrial products that obey the standard of green and sustainable chemistry. Here, we discuss the ability of a full polarizable hybrid approach coupled to a standard molecular dynamics scheme to model the behavior in the aqueous phase and at infinite dilution conditions of a standard hydrophobic polyelectrolyte polymer whose charge is neutralized by explicit counterions. Beyond the standard picture of a polyelectrolyte behavior governed by an interplay between opposite intra-polyelectrolyte and inter-polyelectrolyte/counterion Coulombic effects, our simulations show the key role played by both intra-solute polarization effects and long range solute/solvent electrostatics to stabilize compact globular conformations of that polyelectrolyte. Our full polarizable hybrid modeling approach is thus a new theoretical tool well suited to be used in digital strategies for accelerating innovation for green science, for instance.

Properties of the tetravalent actinide series in aqueous phase from a microscopic simulation self-consistent engine.

In the context of nuclear fuel recycling and environmental issues, the understanding of the properties of radio-elements with various approaches remains a challenge regarding their dangerousness. Moreover, experimentally, some issues are also of importance; first, it is imperative to work at sufficiently high concentrations to reach the sensitivities of the analytical tools, however this condition often leads to precipitation for some of them; second, stabilizing specific oxidation states of some actinides remains a challenge, thus making it difficult to extract general trends across the actinide series. Complementary to experiments, modeling can be used to unbiasedly probe the actinide's properties in an aquatic environment and offers a predictive tool. We report the first molecular dynamics simulations based on homogeneously built force fields for the whole series of the tetravalent actinides in aqueous phase from ThIV to BkIV and including PuIV. The force fields used to model the interactions among the constituents include polarization and charge donation microscopic effects. They are built from a self-consistent iterative ab initio based engine that can be included in future developments as an element of a potential machine learning procedure devoted to generating accurate force fields. The comparison of our simulated hydrated actinide properties to available experimental data shows the model robustness and the relevance of our parameter assignment engine. Moreover, our simulated structural, dynamical and evolution of the hydration free energy data show that, apart from AmIV and CmIV, the actinide properties change progressively along the series.

Ion hydration free energies and water surface potential in water nano drops: The cluster pair approximation and the proton hydration Gibbs free energy in solution.

We estimate both single ion hydration Gibbs free energies in water droplets, comprising from 50 to 1000 molecules, and water/vacuum surface potentials in pure water droplets comprising up to 10 000 molecules. We consider four ions, namely, Li+, NH4 +, F-, and Cl-, and we model their hydration process and water/water interactions using polarizable force fields based on an induced point dipole approach. We show both ion hydration Gibbs free energies and water surface potentials to obey linear functions of the droplet radius as soon as droplets comprising a few hundred water molecules. Moreover, we also show that the differences in anion/cation hydration Gibbs free energies in droplets obey a different regime in large droplets than in small clusters comprising no more than six water molecules, in line with the earlier results computed from standard additive point charge force fields. Hence, both point charge and more sophisticated induced point dipole molecular modeling approaches suggest that methods considering only the thermodynamical properties of small ion/water clusters to estimate the absolute proton hydration Gibbs free energy in solution are questionable. In particular, taking into account the data of large ion/water droplets may yield a proton hydration Gibbs free energy in solution value to be shifted by several kBT units compared to small clusters-based approaches.

Improving the description of solvent pairwise interactions using local solute/solvent three-body functions. The case of halides and carboxylates in aqueous environment.

We propose a general strategy to remediate force-field artifacts in describing pairwise interactions among similar molecules M in the vicinity of another chemical species, C, like water molecules interacting at short distance from a monoatomic ion. This strategy is based on introducing a three-body potential energy term that alters the pairwise interactions among M-type molecules when they lie at short range from the species C. In other words the species C is the center of a space domain where the pairwise interactions among the molecules M is altered. Here, we apply it to improve the description of the water interactions provided by the polarizable water model TCPE/2013 in the vicinity of halides, from F- to At- , and of the prototypical carboxylate anion CH3 COO- . We show the accuracy and the transferability of such an approach to investigate not only the hydration process of single anions but also of a salt solution NH4+/Cl- in aqueous phase. This strategy can be used to remediate the drawbacks of any kind of force fields. © 2019 Wiley Periodicals, Inc.

Solvation of the Guanidinium Ion in Pure Aqueous Environments: A Theoretical Study from an "Ab Initio"-Based Polarizable Force Field.

We report simulation results regarding the hydration process of the guanidinium cation in water droplets and in bulk liquid water, at a low concentration of 0.03 M, performed using a polarizable approach to model both water/water and ion/water interactions. In line with earlier theoretical studies, our simulations show a preferential orientation of guanidinium at water-vacuum interfaces, i.e., a parallel orientation of the guanidinium plane to the aqueous surface. In an apparent contradiction with earlier simulation studies, we show also that guanidinium has a stronger propensity for the cores of aqueous systems than the ammonium cation. However, our bulk simulation conditions correspond to weaker cation concentrations than in earlier studies, by 2 orders of magnitude, and that the same simulations performed using a standard nonpolarizable force field leads to the same conclusion. From droplet data, we extrapolate the guanidinium single hydration enthalpy value to be -82.9 ± 2.2 kcal mol-1. That is about half as large as the sole experimental estimate reported to date, about -144 kcal mol-1. Our result yields a guanidinium absolute bulk hydration free energy at ambiant conditions to be -78.4 ± 2.6 kcal mol-1, a value smaller by 3 kcal mol-1 compared to ammonium. The relatively large magnitude of our guanidinium hydration free energy estimate suggests the Gdm+ protein denaturing properties to result from a competition between the cation hydration effects and the cation/protein interactions, a competition that can be modulated by weak differences in the protein or in the cation chemical environment.

Organic ion association in aqueous phase and ab initio-based force fields: The case of carboxylate/ammonium salts.

We performed molecular dynamics simulations of carboxylate/methylated ammonium ion pairs solvated in bulk water and of carboxylate/methylated ammonium salt solutions at ambient conditions using an ab initio-based polarizable force field whose parameters are assigned to reproduce only high end quantum computations, at the Møller-Plesset second-order perturbation theory/complete basis set limit level, regarding single ions and ion pairs as isolated and micro-hydrated in gas phase. Our results agree with the available experimental results regarding carboxylate/ammonium salt solutions. For instance, our force field approach predicts the percentage of acetate associated with ammonium ions in CH3COO-/CH3NH3+ solutions at the 0.2-0.8M concentration scale to range from 14% to 35%, in line with the estimates computed from the experimental ion association constant in liquid water. Moreover our simulations predict the number of water molecules released from the ion first hydration shell to the bulk upon ion association to be about 2.0 ± 0.6 molecules for acetate/protonated amine ion pairs, 3.1 ± 1.5 molecules for the HCOO-/NH4+ pair and 3.3 ± 1.2 molecules for the CH3COO-/(CH3)4N+ pair. For protonated amine-based ion pairs, these values are in line with experiment for alkali/halide pairs solvated in bulk water. All these results demonstrate the promising feature of ab initio-based force fields, i.e., their capacity in accurately modeling chemical systems that cannot be readily investigated using available experimental techniques.

Complexity Reduction in Large Quantum Systems: Fragment Identification and Population Analysis via a Local Optimized Minimal Basis.

We present, within Kohn-Sham density functional theory calculations, a quantitative method to identify and assess the partitioning of a large quantum-mechanical system into fragments. We then show how within this framework simple generalizations of other well-known population analyses can be used to extract, from first-principles, reliable electrostatic multipoles for the identified fragments. Our approach reduces arbitrariness in the fragmentation procedure and enables the possibility to assess quantitatively whether the corresponding fragment multipoles can be interpreted as observable quantities associated with a system moiety. By applying our formalism within the code BigDFT, we show that the use of a minimal set of in situ-optimized basis functions allows at the same time a proper fragment definition and an accurate description of the electronic structure.

Extrapolating Single Organic Ion Solvation Thermochemistry from Simulated Water Nanodroplets.

We compute the ion/water interaction energies of methylated ammonium cations and alkylated carboxylate anions solvated in large nanodroplets of 10 000 water molecules using 10 ns molecular dynamics simulations and an all-atom polarizable force-field approach. Together with our earlier results concerning the solvation of these organic ions in nanodroplets whose molecular sizes range from 50 to 1000, these new data allow us to discuss the reliability of extrapolating absolute single-ion bulk solvation energies from small ion/water droplets using common power-law functions of cluster size. We show that reliable estimates of these energies can be extrapolated from a small data set comprising the results of three droplets whose sizes are between 100 and 1000 using a basic power-law function of droplet size. This agrees with an earlier conclusion drawn from a model built within the mean spherical framework and paves the road toward a theoretical protocol to systematically compute the solvation energies of complex organic ions.

Structural, dynamical, and transport properties of the hydrated halides: How do At(-) bulk properties compare with those of the other halides, from F(-) to I(-)?

The properties of halides from the lightest, fluoride (F(-)), to the heaviest, astatide (At(-)), have been studied in water using a polarizable force-field approach based on molecular dynamics (MD) simulations at the 10 ns scale. The selected force-field explicitly treats the cooperativity within the halide-water hydrogen bond networks. The force-field parameters have been adjusted to ab initio data on anion/water clusters computed at the relativistic Möller-Plesset second-order perturbation theory level of theory. The anion static polarizabilities of the two heaviest halides, I(-) and At(-), were computed in the gas phase using large and diffuse atomic basis sets, and taking into account both electron correlation and spin-orbit coupling within a four-component framework. Our MD simulation results show the solvation properties of I(-) and At(-) in aqueous phase to be very close. For instance, their first hydration shells are structured and encompass 9.2 and 9.1 water molecules at about 3.70 ± 0.05 Å, respectively. These values have to be compared to the F(-), Cl(-), and Br(-) ones, i.e., 6.3, 8.4, and 9.0 water molecules at 2.74, 3.38, and 3.55 Å, respectively. Moreover our computations predict the solvation free energy of At(-) in liquid water at ambient conditions to be 68 kcal mol(-1), a value also close the I(-) one, about 70 kcal mol(-1). In all, our simulation results for I(-) are in excellent agreement with the latest neutron- and X-ray diffraction studies. Those for the At(-) ion are predictive, as no theoretical or experimental data are available to date.

Simulated Solvation of Organic Ions II: Study of Linear Alkylated Carboxylate Ions in Water Nanodrops and in Liquid Water. Propensity for Air/Water Interface and Convergence to Bulk Solvation Properties.

We investigated the solvation of carboxylate ions from formate to hexanoate, in droplets of 50 to 1000 water molecules and neat water, by computations using standard molecular dynamics and sophisticated polarizable models. The carboxylate ions from methanoate to hexanoate show strong propensity for the air/water interface in small droplets. Only the ions larger than propanoate retain propensity for the interface in larger droplets, where their enthalpic stabilization by ion/water dispersion is reduced there by 3 kcal mol(-1) per CH2 group. This is compensated by entropy effects over +3.3 cal mol(-1) K(-1) per CH2 group. On the surface, the anionic headgroups are strongly oriented toward the aqueous core, while the hydrophobic alkyl chains are repelled into air and lose their structure-making effects. These results reproduce the structure-making effects of alkyl groups in solution, and suggest that the hydrocarbon chains of ionic headgroups and alkyl substituents solvate independently. Extrapolation to bulk solution using standard extrapolation schemes yields absolute carboxylate solvation energies. The results for formate and acetate yield a proton solvation enthalpy of about 270 kcal mol(-1), close to the experiment-based value. The largest carboxylate ions yield a value smaller by about 10 kcal mol(-1), which requires studies in much larger droplets.

The fast multipole method and point dipole moment polarizable force fields.

We present an implementation of the fast multipole method for computing Coulombic electrostatic and polarization forces from polarizable force-fields based on induced point dipole moments. We demonstrate the expected O(N) scaling of that approach by performing single energy point calculations on hexamer protein subunits of the mature HIV-1 capsid. We also show the long time energy conservation in molecular dynamics at the nanosecond scale by performing simulations of a protein complex embedded in a coarse-grained solvent using a standard integrator and a multiple time step integrator. Our tests show the applicability of fast multipole method combined with state-of-the-art chemical models in molecular dynamical systems.

Simulated solvation of organic ions: protonated methylamines in water nanodroplets. Convergence toward bulk properties and the absolute proton solvation enthalpy.

We applied an alternative, purely theoretical route to estimate thermodynamical properties of organic ions in bulk solution. The method performs a large ensemble of simulations of ions solvated in water nanodroplets of different sizes, using a polarizable molecular dynamics approach. We consider protonated ammonia and methylamines, and K(+) for comparison, solvated in droplets of 50-1000 water molecules. The parameters of the model are assigned from high level quantum computations of small clusters. All the bulk phase results extrapolated from droplet simulations match, and confirm independently, the relative and absolute experiment-based ion solvation energies. Without using experiment-based parameters or assumptions, the results confirm independently the solvation enthalpy of the proton, as -270.3 ± 1.1 kcal mol(-1). The calculated relative solvation enthalpies of these ions are constant from small water clusters, where only the ionic headgroups are solvated, up to bulk solution. This agrees with experimental thermochemistry, that the relative solvation energies of alkylammonium ions by only four H2O molecules reproduce the relative bulk solvation energies, although the small clusters lack major bulk solvation factors. The droplet results also show a slow convergence of ion solvation properties toward their bulk limit, and predict that the stepwise solvation enthalpies of ion/water droplets are very close to those of pure neutral water droplets already after 50 water molecules. Both the ionic and neutral clusters approach the bulk condensation energy very gradually up to 10,000 water molecules, consistent with the macroscopic liquid drop model for pure water droplets. Compared to standard computational methods based on infinite periodic systems, our protocol represents a new purely theoretical approach to investigate the solvation properties of ions. It is applicable to the solvation of organic ions, which are pivotal in environmental, industrial, and biophysical chemistry but have been little investigated theoretically up to the present.

Revisiting a many-body model for water based on a single polarizable site: from gas phase clusters to liquid and air/liquid water systems.

We present a revised version of the water many-body model TCPE [M. Masella and J.-P. Flament, J. Chem. Phys. 107, 9105 (1997)], which is based on a static three charge sites and a single polarizable site to model the molecular electrostatic properties of water, and on an anisotropic short range many-body energy term specially designed to accurately model hydrogen bonding in water. The parameters of the revised model, denoted TCPE/2013, are here developed to reproduce the ab initio energetic and geometrical properties of small water clusters (up to hexamers) and the repulsive water interactions occurring in cation first hydration shells. The model parameters have also been refined to reproduce two liquid water properties at ambient conditions, the density and the vaporization enthalpy. Thanks to its computational efficiency, the new model range of applicability was validated by performing simulations of liquid water over a wide range of temperatures and pressures, as well as by investigating water liquid/vapor interfaces over a large range of temperatures. It is shown to reproduce several important water properties at an accurate enough level of precision, such as the existence liquid water density maxima up to a pressure of 1000 atm, the water boiling temperature, the properties of the water critical point (temperature, pressure, and density), and the existence of a "singularity" temperature at about 225 K in the supercooled regime. This model appears thus to be particularly well-suited for characterizing ion hydration properties under different temperature and pressure conditions, as well as in different phases and interfaces.

A multiscale coarse-grained polarizable solvent model for handling long tail bulk electrostatics.

A multiscale coarse-grained approach able to handle efficiently the solvation of microscopic solutes in extended chemical environment is described. That approach is able to compute readily and efficiently very long-range solute/solvent electrostatic microscopic interactions, up to the 1-µm scale, by considering a reduced amount of computational resources. All the required parameters are assigned to reproduce available data concerning the solvation of single ions. Such a strategy makes it possible to reproduce with good accuracy the solvation properties concerning simple ion pairs in solution (in particular, the asymptotic behavior of the ion pair potentials of mean force). This new method represents an extension of the polarizable pseudoparticle solvent model, which has been recently improved to account for the main features of hydrophobic effects in liquid water (Masella et al., J. Comput. Chem. 2011, 32, 2664). This multiscale approach is well suited to be used for computing the impact of charge changes in free energy computations, in terms of both accuracy and efficiency.

Further insights in the ability of classical nonadditive potentials to model actinide ion-water interactions.

Pursuing our efforts on the development of accurate classical models to simulate radionuclides in complex environments (Réal et al., J. Phys. Chem. A 2010, 114, 15913; Trumm et al. J. Chem. Phys. 2012, 136, 044509), this article places a large emphasis on the discussion of the influence of models/parameters uncertainties on the computed structural, dynamical, and temporal properties. Two actinide test cases, trivalent curium and tetravalent thorium, have been studied with three different potential energy functions, which allow us to account for the polarization and charge-transfer effects occurring in hydrated actinide ion systems. The first type of models considers only an additive energy term for modeling ion/water charge-transfer effects, whereas the other two treat cooperative charge-transfer interactions with two different analytical expressions. Model parameters are assigned to reproduce high-level ab initio data concerning only hydrated ion species in gas phase. For the two types of cooperative charge-transfer models, we define two sets of parameters allowing or not to cancel out possible errors inherent to the force field used to model water/water interactions at the ion vicinity. We define thus five different models to characterize the solvation of each ion. For both ions, our cooperative charge-transfer models lead to close results in terms of structure in solution: the coordination number is included within 8 and 9, and the mean ion/water oxygen distances are 2.45 and 2.49 Å, respectively, for Th(IV) and Cm(III).

Modeling the hydration of mono-atomic anions from the gas phase to the bulk phase: the case of the halide ions F-, Cl-, and Br-.

In this work, we investigate the hydration of the halide ions fluoride, chloride, and bromide using classical molecular dynamics simulations at the 10 ns scale and based on a polarizable force-field approach, which treats explicitly the cooperative bond character of strong hydrogen bond networks. We have carried out a thorough analysis of the ab initio data at the MP2 or CCSD(T) level concerning anion/water clusters in gas phase to adjust the force-field parameters. In particular, we consider the anion static polarizabilities computed in gas phase using large atomic basis sets including additional diffuse functions. The information extracted from trajectories in solution shows well structured first hydration shells formed of 6.7, 7.0, and 7.6 water molecules at about 2.78 Å, 3.15 Å, and 3.36 Å for fluoride, chloride, and bromide, respectively. These results are in excellent agreement with the latest neutron- and x-ray diffraction studies. In addition, our model reproduces several other properties of halide ions in solution, such as diffusion coefficients, description of hydration processes, and exchange reactions. Moreover, it is also able to reproduce the electrostatic properties of the anions in solution (in terms of anion dipole moment) as reported by recent ab initio quantum simulations. All the results show the ability of the proposed model in predicting data, as well as the need of accounting explicitly for the cooperative character of strong hydrogen bonds to reproduce ab initio potential energy surfaces in a mean square sense and to build up a reliable force field.