Building Water Models, A Different Approach

Abstract

Simplified, classical models of water are an integral part of atomistic molecular simulations, especially in biology and chemistry where hydration effects are critical. Yet, despite several decades of effort, these models are still far from perfect. Presented here is an alternative approach to constructing point charge water models - currently, the most commonly used type. In contrast to the conventional approach, we do not impose any geometry constraints on the model other than symmetry. Instead, we optimize the distribution of point charges to best describe the "electrostatics" of the water molecule, which is key to many unusual properties of liquid water. The search for the optimal charge distribution is performed in 2D parameter space of key lowest multipole moments of the model, to find best fit to a small set of bulk water properties at room temperature. A virtually exhaustive search is enabled via analytical equations that relate the charge distribution to the multipole moments. The resulting "optimal" 3-charge, 4-point rigid water model (OPC) reproduces a comprehensive set of bulk water properties significantly more accurately than commonly used rigid models: average error relative to experiment is 0.76%. Close agreement with experiment holds over a wide range of temperatures, well outside the ambient conditions at which the fit to experiment was performed. The improvements in the proposed water model extend beyond bulk properties: compared to the common rigid models, predicted hydration free energies of small molecules in OPC water are uniformly closer to experiment, root-mean-square error < 1kcal/mol.

References in corpus (1)

• What ice can teach us about water interactions: a critical comparison of the performance of different water models

Cited by in corpus (23)

• Melting points of water models: Current situation
• Machine-learned molecular mechanics force field for the simulation of protein-ligand systems and beyond
• Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles
• On the design space between molecular mechanics and machine learning force fields
• Ephemeral ice-like local environments in classical rigid models of liquid water
• Accurate Evaluation of Charge Asymmetry in Aqueous Solvation
• Refinement of molecular dynamics ensembles using experimental data and flexible forward models
• Inclusion of Enclosed Hydration Effects in the Binding Free Energy Estimation of Dopamine D3 Receptor Complexes
• Comparison of water models for structure prediction
• A Riemannian Stochastic Representation for Quantifying Model Uncertainties in Molecular Dynamics Simulations
• How important is the dielectric constant in water modeling? Evaluation of the performance of the TIP4P/$\epsilon$ force field and its compatibility with the Joung-Cheatham NaCl model
• Benchmark free energies and entropies for saturated and compressed water
• Non-conformal coarse-grained potentials for water
• Basic stability tests of machine learning potentials for molecular simulations in computational drug discovery
• Molecular dynamics simulations with grand-canonical reweighting suggest cooperativity effects in RNA structure probing experiments
• Static Dielectric Permittivity Profiles and Coarse-graining Approaches for Water in Graphene Slit Pores
• Molecular simulations to investigate the impact of N6-methylation in RNA recognition: Improving accuracy and precision of binding free energy prediction
• MDRefine: a Python package for refining Molecular Dynamics trajectories with experimental data
• Statistical mechanical theory of liquid water
• Hydration Free Energies of Linear Alkanes: Systematic Deviations in Common Water Models and Their Correction
• Automatic learning of hydrogen-bond fixes in an AMBER RNA force field
• Dielectric permittivity of water confined in stacks of charged lipid layers: extracting profiles from molecular dynamics simulations using a modified Poisson-Boltzmann equation
• Influence of water models on water movement through AQP1