Abstract:Protein hydration plays a vital role in many biological functions and molecular simulations are frequently used to study the effect of protein hydration at the atomic level. However, the accuracy of these simulations has often been highly sensitive to the water model used, a phenomenon best known in the case of intrinsically disordered proteins. In the present study, we investigated the extent to which the choice of water model alters the behavior of complex networks of transport tunnels. Tunnels are essential because they allow substrates and products to access and exit the active sites of enzymes that are otherwise deeply embedded within the protein structure. The ability of these tunnels to regulate access directly affects enzyme efficiency and selectivity, making their study crucial for understanding enzyme function and inhibition at a mechanistic level. By performing all-atom molecular dynamics simulations of the wild-type haloalkane dehalogenase LinBWT and its two variants, LinB32 and LinB86, with synthetically engineered tunnel networks in TIP3P and OPC water models, we investigated the effects of these models on the overall tunnel topology. We also analyzed the properties of the main tunnels, such as their conformation, bottleneck dimensions, sampling efficiency, and duration of the tunnel opening. Our data demonstrate that all three proteins exhibited similar conformational behavior in both water models and differed in the geometrical characteristics of their auxiliary tunnels, in line with experimental observations. Interestingly, the results indicate that the stability of the open tunnels is sensitive to the water model and the system under question. Our findings suggest that the 3-point TIP3P model can provide comparable inference on the overall topology of the networks of primary tunnels and their geometry, and thus may be a desirable choice when computational resources are limited or when compatibility issues impede usage of OPC with certain protein force fields. However, when a more thorough investigation is performed, such as the calculation of ligand unbinding rates via such tunnel networks, where precision and intricate details are paramount, the more costly 4-point OPC model would be more suited.

Accuracy of TIP4P/2005 and SPC/Fw Water Models

Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations.

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

Modeling Solution Behavior of Poly(N-isopropylacrylamide): A Comparison between Water Models

The shear viscosities of common water models by non-equilibrium molecular dynamics simulations

A Comparative Study of TIP4P-2005, SPC/E, SPC, and TIP3P-Ew Models for Predicting Water Transport Coefficients Using EMD and NEMD Simulations

Seven-Site Effective Pair Potential for Simulating Liquid Water

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

Faster Sampling in Molecular Dynamics Simulations with TIP3P-F Water

Fine Tuning the Intermolecular Interactions of Water Clusters Using the Dispersion-Corrected Density Functional Theory

On the composition dependence of the microscopic structure, thermodynamic, dynamic and dielectric properties of water-dimethyl formamide model mixtures. Molecular dynamics simulation results

On the composition dependence of thermodynamic, dynamic and dielectric properties of water-methanol model mixtures. Molecular dynamics simulation results

Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Getting the intermolecular forces correct: introducing the ASTA strategy for a water model

On the Choice of Different Water Model in Molecular Dynamics Simulations of Nanopore Transport Phenomena

Impact of water models on structure and dynamics of enzyme tunnels

On the transferability of three water models developed by adaptive force matching

Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Influence of water models on water movement through AQP1

Polarizable Mean-Field Model of Water for Biological Simulations with Amber and Charmm force fields