Water models for biomolecular simulations

Abstract

Modern simulation and modeling approaches to investigation of biomolecular structure and function rely heavily on a variety of methods—water models—to approximate the influence of solvent. We give a brief overview of several distinct classes of available water models, with the emphasis on the conceptual basis at each level of approximation. The main focus is on classes of models most widely used in atomistic simulations, including popular implicit and explicit solvent models. Among the latter, nonpolarizable N‐point models are covered in most detail, including some recent methodological advances and nuances. Notes on practical availability and usage in biomolecular simulations are included. Atomistic simulations that were hardly possible only a short while ago have revealed significant problems that can be traced to deficiencies of most commonly used N‐point water models. Recently developed models of this class approximate experimental properties of liquid water much closer than before, and show promise in practical biomolecular simulations. Obstacles to wider adoption of these more accurate water models, both technical and conceptual, are discussed. It is argued that verifying robustness of simulation results to the choice of water model can be of immediate benefit even in the absence of a clear replacement for older models; a specific strategy is proposed. The review is concluded with a discussion on how force‐field development efforts can benefit from better solvent models, and vice versa. WIREs Comput Mol Sci 2018, 8:e1347. doi: 10.1002/wcms.1347 This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

“…SHAKE or SETTLE), are used in combination for increased computational efficiency. However, more complex, and presumably more accurate, models such as TIP4P-Ew (Horn et al 2004) are being increasingly utilized in biomolecular simulations (Onufriev & Izadi, 2018). Ion parameters are intimately tied to water models (Joung & Cheatham III, 2008), and all the force fields representing system components should be compatible with one another in order to provide an accurate balance of protein-water-lipid-ion interactions that allows for a complete description of ion permeation and channel gating.…”

Intuitive, reproducible high-throughput molecular dynamics in Galaxy: a tutorial

“…Apart from the force field, a water model also needs to be selected to model the solvent; a wide range of models exist for this purpose. Here we are using the common TIP3P model, which is an example of a 'three-site model' -socalled because the molecule is modeled using three points, corresponding to the three atoms of water (Four-and five-site models include additional 'dummy atoms' representing the negative charges of the lone pairs of the oxygen atom) (22).…”

Force‐field parametrization based on radial and energy distribution functions

“…Parametrization of FFs is considered essential, and various methods and variations of parameters are still being developed actively. For example, the water molecule is the most active target of parameterization . A study regarding the SPC/E water model reported that a 1% increase in the length of the OH bond led to accurate reproducibility of the translational and rotational diffusion of pure water and proteins.…”

“…For example, the water molecule is the most active target of parameterization. [6] A study regarding the SPC/E water model [7] reported that a 1% increase in the length of the OH bond led to accurate reproducibility of the translational and rotational diffusion of pure water and proteins. This newly developed model was named SPC/E b .…”