Prediction techniques free energy modeling

The present chapter thus provides an overview of the current status of continuum models of solvation. We review available continuum models and computational techniques implementing such models for both electrostatic and non-electrostatic components of the free energy of solvation. We then consider a number of case studies, with particular focus on the prediction of heterocyclic tautomeric equilibria. In the discussion of the latter we center attention on the subtleties of actual chemical systems and some of the dangers of applying continuum models uncritically. We hope the reader will emerge with a balanced appreciation of the power and limitations of these methods. 

Rahaman and Hatton [152] developed a thermodynamic model for the prediction of the sizes of the protein filled and unfilled RMs as a function of system parameters such as ionic strength, protein charge, and size, Wq and protein concentration for both phase transfer and injection techniques. The important assumptions considered include (i) reverse micellar population is bidisperse, (ii) charge distribution is uniform, (iii) electrostatic interactions within a micelle and between a protein and micellar interface are represented by nonlinear Poisson-Boltzmann equation, (iv) the equilibrium micellar radii are assumed to be those that minimize the system free energy, and (v) water transferred between the two phases is too small to change chemical potential. 

Liquid/liquid partition constants within pharmaceutical chemistry have been of primary interest because of tlieir correlation with liquid/membrane partitioning behavior. A sufficiently fluid membrane may, in some sense, be regarded as a solvent. With such an outlook, tlie partitioning phenomenon may again be regarded as a liquid/liquid example, amenable to treatment with standard continuum techniques. Of course, accurate continuum solvation models typically rely on the availabihty of solvation free energies or bulk solvent properties in order to develop useful parameterizations, and such data may be sparse or non-existent for membranes. Some success, however, has been demonstrated for predicting such data either by intuitive or statistical analysis (see, for example. Chambers etal. 1999). 

In protein structure prediction, potentials are used to assign an energy-like quantity to a conformation of a protein molecule. If this quantity enables us to distinguish the native state of a protein, the potential is regarded as a reasonable model for a protein-solvent system. The rationale behind this relies on two assumptions (a) a solved protein in its native state can be described by an ensemble of closely related conformations, and (b) in this state the system is in the global minimum of free energy. Virtually all techniques designed for structure prediction are based on these principles [3,4]. 

Nearly all theories to date predict that IETS intensities should be proportional to n, the surface density of molecular scatterers. Langan and Hansma (21) used radioactively labeled chemicals to measure a surface concentration vs solution concentration curve ( Fig. 10 ) for benzoic acid on alumina using the liquid doping technique. The dashed line in Fig. 10 is a 2 parameter fit to the data using a simple statistical mechanical model by Cederberg and Kirtley (35). This model matched the free energy of the molecule on the surface with that in solution. The two parameters in this model were the surface density of binding sites ( 10" A )... 

Here, we review how the development of SAXS as a structural technique is driven by advances in computer algorithms that allow to reconstruct low-resolution electron density maps ab initio from scattering profiles. In addition, we delineate how these low-resolution models can be used in free energy electrostatics calculations. Finally, we discuss how one can exploit the hierarchical nature of RNA folding by combining the low resolution, global information provided by SAXS with local information on RNA structure, from either experiments or state-of-the-art RNA structure prediction algorithms, to further increase the resolution and quality of models obtained from SAXS. 

The above techniques have been used in numerous calculations of solute free energy profiles. Wilson and Pohorille [52] and Benjamin[53] have determined the free energy profiles for small ions at the water liquid/vapor interface and compared the results to predictions of continuum electrostatic models. The transfer of small ions to the interface involves a monotonic increase in the free energy which is in qualitative agreement with the continuum model. This behavior is consistent with the increase in the surface tension of water with the increase in the concentration of a very dilute salt solution, and it represents the fact that small ions are repelled from the liquid/vapor interface. On the other hand, calculations of the free energy profile at the water liquid/vapor interface of hydrophobic molecules, such as phenol[54] and pentyl phenol[57] and even molecules such as ethanol [58], show that these molecules are attracted to the surface region and lower the surface tension of water. In addition, the adsorption free energy of solutes at liquid/liquid interfaces[59,60] and at water/metal interfaces[61-64] have been reported. 

Because the WS mixing rule uses VLE information only at low pressure, it can also be used to make predictions at high pressure based on low-pressure prediction techniques such as UNIFAC and other group contribution methods (Orbey, Sandler, and Wong 1993). This completely predictive method using the WS and other excess free-energy-based EOS models is discussed in Chapter 5. 

---