The issues of electronic polarizability in molecular dynamics simulations are discussed. We argue that the charges of ionized groups in proteins, and charges of ions in conventional non-polarizable force fields such as CHARMM, AMBER, GROMOS, etc should be scaled by a factor about 0.7. Our model explains why a neglect of electronic solvation energy, which typically amounts to about a half of total solvation energy, in non-polarizable simulations with un-scaled charges can produce a correct result; however, the correct solvation energy of ions does not guarantee the correctness of ion-ion pair interactions in many non-polarizable simulations. The inclusion of electronic screening for charged moieties is shown to result in significant changes in protein dynamics and can give rise to new qualitative results compared with the traditional non-polarizable force field simulations. The model also explains the striking difference between the value of water dipole μ∼ 3D reported in recent ab initio and experimental studies with the value μ(eff)∼ 2.3D typically used in the empirical potentials, such as TIP3P or SPC/E. It is shown that the effective dipole of water can be understood as a scaled value μ(eff) = μ/√ε(el), where ε(el) = 1.78 is the electronic (high-frequency) dielectric constant of water. This simple theoretical framework provides important insights into the nature of the effective parameters, which is crucial when the computational models of liquid water are used for simulations in different environments, such as proteins, or for interaction with solutes.

A simple model for accounting for electronic polarization in molecular dynamics (MD) simulations is discussed. In this model, called molecular dynamics electronic continuum (MDEC), the electronic polarization is treated explicitly in terms of the electronic continuum (EC) approximation, while the nuclear dynamics is described with a fixed-charge force field. In such a force-field all atomic charges are scaled to reflect the screening effect by the electronic continuum. The MDEC model is rather similar but not equivalent to the standard nonpolarizable force-fields; the differences are discussed. Of our particular interest is the calculation of the electrostatic part of solvation energy using standard nonpolarizable MD simulations. In a low-dielectric environment, such as protein, the standard MD approach produces qualitatively wrong results. The difficulty is in mistreatment of the electronic polarizability. We show how the results can be much improved using the MDEC approach. We also show how the dielectric constant of the medium obtained in a MD simulation with nonpolarizable force-field is related to the static (total) dielectric constant, which includes both the nuclear and electronic relaxation effects. Using the MDEC model, we discuss recent calculations of dielectric constants of alcohols and alkanes, and show that the MDEC results are comparable with those obtained with the polarizable Drude oscillator model. The applicability of the method to calculations of dielectric properties of proteins is discussed.

Electronic polarizability is an important factor in molecular interactions. In the conventional force fields such as AMBER or CHARMM, however, there is inconsistency in how the effect of electronic dielectric screening of Coulombic interactions, inherent for the condensed phase media, is treated. Namely, the screening appears to be accounted for via effective charges only for neutral moieties, whereas the charged residues are treated as if they were in vacuum. As a result, the electrostatic interactions between ionized groups are exaggerated in molecular simulations by the factor of about 2. The discussed here MDEC (Molecular Dynamics in Electronic Continuum) model provides a theoretical framework for modification of the standard non-polarizable force fields to make them consistent with the idea of uniform electronic screening of partial atomic charges. The present theory states that the charges of ionized groups and ions should be scaled; i.e. reduced by a factor of about 0.7. In several examples, including the interaction between Na+ ions, which is of interest for ion-channel simulations, and the dynamics of an important salt-bride in Cytochrome c Oxidase, we compared the standard non-polarizable MD simulations with MDEC simulations, and demonstrated that MDEC charge scaling procedure results in more accurate interactions. The inclusion of electronic screening for charged moieties is shown to result in significant changes in protein dynamics and can give rise to new qualitative results compared with the traditional non-polarizable force fields simulations.

Employing the continuum dielectric model for electronic polarizability, we have developed a new consistent procedure for parameterization of the effective nonpolarizable potential of liquid water. The model explains the striking difference between the value of water dipole moment μ~3D reported in recent ab initio and experimental studies with the value μeff~2.3D typically used in the empirical potentials, such as TIP3P or SPC/E. It is shown that the consistency of the parameterization scheme can be achieved if the magnitude of the effective dipole of water is understood as a scaled value μeff=μ∕εel, where εel =1.78 is the electronic (high-frequency) dielectric constant of water, and a new electronic polarization energy term, missing in the previous theories, is included. The new term is evaluated by using Kirkwood - Onsager theory. The new scheme is fully consistent with experimental data on enthalpy of vaporization, density, diffusion coefficient, and static dielectric constant. The new theoretical framework provides important insights into the nature of the effective parameters, which is crucial when the computational models of liquid water are used for simulations in different environments, such as proteins, or for interaction with solutes.

Although a great number of computational models of water are available today, the majority of current biological simulations are done with simple models, such as TIP3P and SPC, developed almost thirty years ago and only slightly modified since then. The reason is that the non-polarizable force fields that are mostly used to describe proteins and other biological molecules are incompatible with more sophisticated modern polarizable models of water. The issue is electronic polarizability: in liquid state, in protein, and in vacuum the water molecule is polarized differently, and therefore has different properties; thus the only way to describe all these different media with the same model is to use a polarizable water model. However, to be compatible with the force field of the rest of the system, e.g. a protein, the latter should be polarizable as well. Here we describe a novel model of water that is in effect polarizable, and yet compatible with the standard non-polarizable force fields such as AMBER, CHARMM, GROMOS, OPLS, etc. Thus the model resolves the outstanding problem of incompatibility.