Models solvent

Specific solute-solvent interactions involving the first solvation shell only can be treated in detail by discrete solvent models. The various approaches like point charge models, siipennoleciilar calculations, quantum theories of reactions in solution, and their implementations in Monte Carlo methods and molecular dynamics simulations like the Car-Parrinello method are discussed elsewhere in this encyclopedia. Here only some points will be briefly mentioned that seem of relevance for later sections. 

Aguilar M A and Olivares del Valle F J 1989 A computation procedure for the dispersion component of the interaction energy in continuum solute solvent models Ohem. Rhys. 138 327-36... 

Sitkoff, D., Sharp, K. A., Honig, B. Accurate calculation of hydration free energies using macroscopic solvent models. J. Phys. Chem. 98 (1994) 1978-1988... 

It is often the case that the solvent acts as a bulk medium, which affects the solute mainly by its dielectric properties. Therefore, as in the case of electrostatic shielding presented above, explicitly defined solvent molecules do not have to be present. In fact, the bulk can be considered as perturbing the molecule in the gas phase , leading to so-called continuum solvent models [14, 15]. To represent the electrostatic contribution to the free energy of solvation, the generalized Bom (GB) method is widely used. Wilhin the GB equation, AG equals the difference between and the vacuum Coulomb energy (Eq. (38)) ... 

The explicit definition of water molecules seems to be the best way to represent the bulk properties of the solvent correctly. If only a thin layer of explicitly defined solvent molecules is used (due to hmited computational resources), difficulties may rise to reproduce the bulk behavior of water, especially near the border with the vacuum. Even with the definition of a full solvent environment the results depend on the model used for this purpose. In the relative simple case of TIP3P and SPC, which are widely and successfully used, the atoms of the water molecule have fixed charges and fixed relative orientation. Even without internal motions and the charge polarization ability, TIP3P reproduces the bulk properties of water quite well. For a further discussion of other available solvent models, readers are referred to Chapter VII, Section 1.3.2 of the Handbook. Unfortunately, the more sophisticated the water models are (to reproduce the physical properties and thermodynamics of this outstanding solvent correctly), the more impractical they are for being used within molecular dynamics simulations. 

Sitkoff D, K A Sharp and B Honig 1994. Accurate Calculation of Hydration Free Energies Usin Macroscopic Solvent Models. Journal of Physical Chemistry 98 1978-1988. 

C. Explicit Solvent Models and the Importance of Balancing the External Interactions... 

In this chapter we provide an introductory overview of the imphcit solvent models commonly used in biomolecular simulations. A number of questions concerning the formulation and development of imphcit solvent models are addressed. In Section II, we begin by providing a rigorous fonmilation of imphcit solvent from statistical mechanics. In addition, the fundamental concept of the potential of mean force (PMF) is introduced. In Section III, a decomposition of the PMF in terms of nonpolar and electrostatic contributions is elaborated. Owing to its importance in biophysics. Section IV is devoted entirely to classical continuum electrostatics. For the sake of completeness, other computational... 

As a first step, it is important to estabUsh implicit solvent models on fundamental principles. For the sake of concreteness, let us consider a solute u immersed in a bulk solution V. The configuration of the solute is represented by the vector X = xj, Xo,.... All other degrees of freedom of the bulk solution surrounding the solute, which may include solvent... 

In the following sections, we describe an implicit solvent model based on this free energy decomposition that is widely used in biophysics. It consists in representing the nonpolar free energy contributions on the basis of the solvent-accessible surface area... 

SASA), a concept introduced by Lee and Richards [9], and the electrostatic free energy contribution on the basis of the Poisson-Boltzmann (PB) equation of macroscopic electrostatics, an idea that goes back to Born [10], Debye and Htickel [11], Kirkwood [12], and Onsager [13]. The combination of these two approximations forms the SASA/PB implicit solvent model. In the next section we analyze the microscopic significance of the nonpolar and electrostatic free energy contributions and describe the SASA/PB implicit solvent model. 

In Section III we described an approximation to the nonpolar free energy contribution based on the concept of the solvent-accessible surface area (SASA) [see Eq. (15)]. In the SASA/PB implicit solvent model, the nonpolar free energy contribution is complemented by a macroscopic continuum electrostatic calculation based on the PB equation, thus yielding an approximation to the total free energy, AVP = A different implicit... 

B Roux, T Simonson, eds. Implicit Solvent Models for Biomolecular Simulations. Special Issue of Biophys Chem Amsterdam Elsevier, 1999. 

For solvent models where the cavity/dispersion interaction is parameterized by fitting to experimental solvation energies, the use of a few explicit solvent molecules for the first solvation sphere is not recommended, as the parameterization represents a best fit to experimental data without any explicit solvent present. 

A classical description of M can for example be a standard force field with (partial) atomic charges, while a quantum description involves calculation of the electronic wave function. The latter may be either a semi-empirical model, such as AMI or PM3, or any of the ab initio methods, i.e. HF, MCSCF, CISD, MP2 etc. Although the electrostatic potential can be derived directly from the electronic wave function, it is usually fitted to a set of atomic charges or multipoles, as discussed in Section 9.2, which then are used in the actual solvent model. 

The mixed solvent models, where the first solvation sphere is accounted for by including a number of solvent molecules, implicitly include the solute-solvent cavity/ dispersion terms, although the corresponding tenns between the solvent molecules and the continuum are usually neglected. Once discrete solvent molecules are included, however, the problem of configuration sampling arises. Nevertheless, in many cases the first solvation shell is by far the most important, and mixed models may yield substantially better results than pure continuum models, at the price of an increase in computational cost. 

CHEMICAL REACTIONS IN THE GAS PHASE AND IN SIMPLE SOLVENT MODELS... 

Exercise 2.4. Repeat Exercise 1.6 with the LD solvent model instead of the external charge. 

---