Development of Coarse‐Grained Force Field by Combining Multilinear Interpolation Technique and Simplex Algorithm

Development of Multiscale Force Field for Actinide (An3+) Solutions

A multiscale force field (FF) is developed for an aqueous solution of trivalent actinide cations An3+ (An = U, Np, Pu, Am, Cm, Bk, and Cf) by using a 12-6-4 Lennard-Jones type potential considering ion-induced dipole interaction. Potential parameters are rigorously and automatically optimized by the meta-multilinear interpolation parametrization (meta-MIP) algorithm via matching the experimental properties, including ion-oxygen distance (IOD) and coordination number (CN) in the first solvation shell and hydration free energy (HFE). The water solvent models incorporate an especially developed polar coarse-grained (CG) water scheme named PW32 and three widely used all-atom (AA) level SPC/E, TIP3P, and TIP4P water schemes. Each PW32 is modeled as two bonded beads to represent three neighboring water molecules, the simulation efficiency of which is 1 to 2 orders of magnitude higher than that of AA waters. The newly developed FF shows high accuracy and transferability in reproducing the IOD, CN, and HFE of An3+. The molecular structure and water exchange dynamics of the first An3+ hydration shell and the ionic (van der Waals) radii are reinvestigated in this work. Moreover, the new FF can readily be transferred to other popular FFs, as it has practicably predicted the permeability of An3+ in a graphene oxide filter within the framework of optimized potentials for liquid simulations (OPLS)-AA FF. It holds promise for applications in exploring actinide aqueous solutions with multiscale computational overhead.

Evaluation of Sampling Algorithms Used for Bayesian Uncertainty Quantification of Molecular Dynamics Force Fields

New Bayesian parameter estimation methods have the capability to enable more physically realistic and reliable molecular dynamics (MD) simulations by providing accurate estimates of uncertainties of force-field (FF) parameters and associated properties. However, the choice of which Bayesian parameter estimation algorithm to use has not been widely investigated, despite its impact on the effective exploration of parameter space. Here, using a case example of the Embedded Atom Method (EAM) FF parameters, we investigated the ramifications of several of the algorithm choices. We found that Ensemble Slice Sampling (ESS) and Affine-Invariant Ensemble Sampling (AIES) demonstrate a new level of superior performance, culminating in more accurate parameter and property estimations with tighter uncertainty bounds, compared to traditional methods such as Metropolis-Hastings (MH), Gradient Search (GS), and Uniform Random Sampler (URS). We demonstrate that Bayesian Uncertainty Quantification with ESS and AIES leads to significantly more accurate and reliable predictions of the FF parameters and properties. The results suggest that ESS and AIES should be used to obtain more accurate parameter and uncertainty estimations while providing deeper physical insights.

Structural and thermodynamic properties of bulk triglycerides and triglyceride/water mixtures reproduced using a polarizable coarse-grained model

Triglycerides (TGs) play important roles in renewable energies, food production, medicine, and metabolism in organisms. Here, we developed a novel coarse-grained (CG) force field (FF) for triglycerides to reproduce both the structural and thermodynamic properties of bulk TGs, TG/air interfaces, and TG/water mixtures using molecular dynamics (MD) simulations. We rigorously optimized the bonded and nonbonded force parameters between the CG beads of TGs and nonbonded force parameters between TG beads and polarizable CG water beads by employing an efficient meta-multilinear interpolation parameterization algorithm recently developed by us. This CG FF performs very well in reproducing the percolating network of the TG bulk phase self-assembled in water and a variety of molecular conformations predicted by all-atom MD simulations. More importantly, it also correctly reproduces multiple experimentally measurable macroscopic thermodynamic properties, including the density and surface tensions of both the TG/air and TG/water interfaces. This paves the way for studying more complicated systems involving TGs on a large scale.

A top-down and bottom-up combined strategy for parameterization of coarse-grained force fields for phospholipids

Coarse-graining (CG) molecular dynamics (MD) simulations are widely used in interpreting experimental observations and predicting assembly morphology as well as collective behaviour but also face the problem of poor accuracy. A main issue is that cross-termed interactions between different CG beads are inadequately parameterized. This work proposes a novel top-down and bottom-up combined strategy to parameterize both self- and cross-termed interactions of zwitterionic phospholipids in water solution based on a piecewise Morse potential describing nonbonded van der Waals interactions. The self-interacting force parameters were optimized by matching experimental density, heat vapourization, and surface tension in a top-down manner, while the cross-termed interactions were optimized by fitting pseudo properties obtained from atomistic simulations in a bottom-up way, including mixing density, intermolecular energy, and radial mixing coefficient. The transferability of the CG force field (FF) was confirmed by reproducing a variety of structural and thermodynamic properties of lipid membranes in both liquid and gel phases. This FF can well depict vesicle self-assembly and vesicle fusion processes. Matching pseudo properties opens a new way to develop CG FF with increased accuracy and transferability.

Is an Inductive Effect Explicit Account Required for Atomic Charges Aimed at Use within the Force Fields?

Polarization and inductive effects are the concepts that have been widely used in qualitative and even quantitative descriptions of experimentally observed properties in chemistry. The polarization effect has proven to be important in cases of biomolecular modeling though still the vast majority of molecular simulations use the classical non-polarizable force fields. In the last few decades, a lot of effort has been put into promoting the polarization effect and incorporating it into modern force fields and charge calculation methods. In contrast, the inductive effect has not attracted such attention and is effectively absent in both classic and modern force fields. Thus, a question is whether this difference corresponds to the difference in the physical significance of the effects and their explicit account, or is an artifact that should be corrected in the next generation of force fields. The significance of the electronic effects is studied in this paper through the prism of performance of specific models for atomic charge calculation that take into explicit account a nested set of effects: the formal charge, the nearest neighbors, the inductive effect, and finally the model, which takes into account all effects, which are possible to account for using atomic charges. The specific choice for the methods is the following: formal charges, MMFF94 bond charge increments, Dynamic Electronegativity Relaxation (DENR), and RESP. We propose a special scheme for the separate estimation of each particular effect contribution. By pairwise comparing the residual molecular electrostatic potential (MEP) errors of those charge models (aimed at best reproducing the quantum chemical reference MEP), we sequentially revealed how the account of each effect contributes to the better-quality MEP reproduction. The following relative importance of effects was estimated; thus, the natural hierarchy of the effects was established. First, the account of formal charges is of primordial importance. Second, the nearest neighbors account is the next in significance. Third, the explicit account of inductive effect in empirical charge calculation schemes was shown to significantly─both qualitatively and quantitatively─improve the quality of MEP reproduction. Fourth, the contribution of polarization is indirectly assessed. Surprisingly, it is of the order of magnitude of the inductive effect even for the molecular systems, for which it is anticipated to be more significant. Finally, the relative importance of anisotropic effects in neutral molecules was additionally reviewed.

Development of accurate coarse-grained force fields for weakly polar groups by an indirect parameterization strategy

Coarse-grained (CG) molecular dynamics simulations are widely used to predict morphological structures and interpret mechanisms of mesoscopic behavior between the scope of traditional experiments and all-atom simulations. However, most current CG force fields (FFs) are not precise enough, especially for polar molecules or functional groups. A main obstacle in developing accurate CG FFs for polar molecules is the freezing problem met at room temperature. In this work, we introduce an indirect parametrization strategy for weakly polar groups by considering their short-chain homologs to avoid freezing. Here, a polar group containing three to four heavy atoms is mapped into one CG bead that is connected to one alkyl bead composed of three or four carbons. The CG beads interact via 4-parameter nonbonded Morse potentials and harmonic bonded potentials. An efficient meta-multilinear interpolation parameterization algorithm, as recently developed by us, is used to rigorously optimize the force parameters. Satisfactory accuracy is witnessed in terms of the density, heat of vaporization, surface tension, and solvation free energy of the homologs of twelve polar molecules, all deviating from the experiment by less than 5%. The transferability of the current FF is indicated by the predicted density, heat of vaporization, and end-to-end distance distributions of fatty acid methyl esters composed of multiple functional groups parameterized in this work.

Development of coarse-grained force field for alcohols: an efficient meta-multilinear interpolation parameterization algorithm

Coarse-grained (CG) molecular dynamics are powerful tools to access a mesoscopic phenomenon and simultaneously record microscopic details, but currently the CG force fields (FFs) are still limited by low parameterization efficiency and poor accuracy especially for polar molecules. In this work, we developed a Meta-Multilinear Interpolation Parameterization (Meta-MIP) algorithm to optimize the CG FFs for alcohols. This algorithm significantly boosts parameterization efficiency by constructing on-the-fly local databases to cover the global optimal parameterization path. In specific, an alcohol molecule is mapped to a heterologous model composed of an OH bead and a hydrocarbon portion which consists of alkane beads representing two to four carbon atoms. Non-bonded potentials are described by soft Morse functions that have no tail-corrections but can still retain good continuities at truncation distance. Nearly all of the properties in terms of density, heat of vaporization, surface tension, and solvation free energy for alcohols predicted by the current FFs deviate from experimental values by less than 7%. This Meta-MIP algorithm can be readily applied to force field development for a wide variety of molecules or functional groups, in many situations including but not limited to CG FFs.

Accurate MP2-based force fields predict hydration free energies for simple alkanes and alcohols in good agreement with experiments

Force fields for four small molecules, methane, ethane, methanol, and ethanol, were created by force matching MP2 gradients computed with triple-zeta-quality basis sets using the Adaptive Force Matching method. Without fitting to any experimental properties, the force fields created were able to predict hydration free energies, enthalpies of hydration, and diffusion constants in excellent agreements with experiments. The root mean square error for the predicted hydration free energies is within 1 kJ/mol of experimental measurements of Ben-Naim et al. [J. Chem. Phys. 81(4), 2016-2027 (1984)]. The good prediction of hydration free energies is particularly noteworthy, as it is an important fundamental property. Similar hydration free energies of ethane relative to methane and of ethanol relative to methanol are attributed to a near cancellation of cavitation penalty and favorable contributions from dispersion and Coulombic interactions as a result of the additional methyl group.