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.

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.

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.

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.

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.

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.

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.

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).

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.

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.

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.

Combining a polarizable force-field and a coarse-grained polarizable solvent model. II. Accounting for hydrophobic effects.

A revised and improved version of our efficient polarizable force-field/coarse grained solvent combined approach (Masella, Borgis, and Cuniasse, J. Comput. Chem. 2008, 29, 1707) is described. The polarizable pseudo-particle solvent model represents the macroscopic solvent polarization by induced dipoles placed on mobile pseudo-particles. In this study, we propose a new formulation of the energy term handling the nonelectrostatic interactions among the pseudo-particles. This term is now able to reproduce the energetic and structural response of liquid water due to the presence of a hydrophobic spherical cavity. Accordingly, the parameters of the energy term handling the nonpolar solute/solvent interactions have been refined to reproduce the free-solvation energy of small solutes, based on a standard thermodynamic integration scheme. The reliability of this new approach has been checked for the properties of solvated methane and of the solvated methane dimer, as well as by performing 10 × 20 ns molecular dynamics (MD) trajectories for three solvated proteins. A long-time stability of the protein structures along the trajectories is observed. Moreover, our method still provides a measure of the protein solvation thermodynamic at the same accuracy as standard Poisson-Boltzman continuum methods. These results show the relevance of our approach and its applicability to massively coupled MD schemes to accurately and intensively explore solvated macromolecule potential energy surfaces.

Structural and atoms-in-molecules analysis of hydrogen-bond network around nitroxides in liquid water.

In this study, we investigated the hydrogen-bond network patterns involving the NO moieties of five small nitroxides in liquid water by analyzing nanosecond scale molecular dynamics trajectories. To this end, we implemented two types of hydrogen-bond definitions, based on electronic structure, using Bader's atoms-in-molecules analysis and based on geometric criteria. In each definition framework, the nitroxide/water hydrogen-bond networks appear very variable from a nitroxide to another. Moreover, each definition clearly leads to a different picture of nitroxide hydration. For instance, the electronic structure-based definition predicts a number of hydrogen bonds around the nitroxide NO moiety usually larger than geometric structure-based ones. One particularly interesting result is that the strength of a nitroxide/water hydrogen bond does not depend on its linearity, leading us to question the relevance of geometric definition based on angular cutoffs to study this type of hydrogen bond. Moreover, none of the hydrogen-bond definitions we consider in the present study is able to quantitatively correlate the strength of nitroxide/water hydrogen-bond networks with the aqueous nitroxide spin properties. This clearly exhibits that the hydrogen-bonding concept is not reliable enough to draw quantitative conclusions concerning such properties.

Combining a polarizable force-field and a coarse-grained polarizable solvent model: application to long dynamics simulations of bovine pancreatic trypsin inhibitor.

The dynamic coupling between a polarizable protein force field and a particle-based implicit solvent model is described. The polarizable force field, TCPEp, developed recently to simulate protein systems, is characterized by a reduced number of polarizable sites, with a substantial gain in efficiency for an equal chemical accuracy. The Polarizable Pseudo-Particle (PPP) solvent model represents the macroscopic solvent polarization by induced dipoles placed on mobile Lennard-Jones pseudo-particles. The solvent-induced dipoles are sensitive to the solute electric field, but not to each other, so that the computational cost of solvent-solvent interactions is basically negligible. The solute and solvent induced dipoles are determined self-consistently and the equations of motion are solved using an efficient iterative multiple time step procedure. The solvation cost with respect to vacuum simulations is shown to decrease with solute size: the estimated multiplicative factor is 2.5 for a protein containing about 1000 atoms, and as low as 1.15 for 8000 atoms. The model is tested for six 20 ns molecular dynamics trajectories of a traditional benchmark system: the hydrated Bovine Pancreatic Trypsin Inhibitor (BPTI). Even though the TCPEp parameters have not been refined to be used with the solvent PPP model, we observe a good conservation of the BPTI structure along the trajectories. Moreover, our approach is able to provide a description of the protein solvation thermodynamic at the same accuracy as the standard Poisson-Boltzman continuum methods. It provides in addition a good description of the microscopic structural aspects concerning the solute/solvent interaction.

Prediction of nitroxide hyperfine coupling constants in solution from combined nanosecond scale simulations and quantum computations.

We present a combined theoretical approach based on analyzing molecular dynamics trajectories (at the nanosecond scale) generated by use of classical polarizable force fields and on quantum calculations to compute averaged hyperfine coupling constants. That method is used to estimate the constant of a prototypical nitroxide: the dimethylnitroxide. The molecule is embedded during the simulations in a cubic box containing about 500 water molecules and the molecular dynamics is generated using periodic conditions. Once the trajectories are achieved, the nitroxide and its first hydration shell molecules are extracted, and the coupling constants are computed by considering the latter aggregates by means of quantum computations. However, all the water molecules of the bulk are also accounted for during those computations by means of the electrostatic potential fitted method. Our results exhibit that in order to predict accurate and reliable coupling constants, one needs to describe carefully the out-of-plane motion of the nitroxide nitrogen and to sample trajectories with a time interval of 400 fs at least to generate an uncorrelated large set of nitroxide structures. Compared to Car-Parrinello molecular dynamics techniques, our approach can be used readily to compute hyperfine coupling constants of large systems, such as nitroxides of great size interacting with macromolecules such as proteins or polymers.

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.