Electrostatics coefficients

To summarize the attempts to refine the original LIE model, we found that an optimal equation for the binding free energy could be obtained with only one free parameter (a) and with the electrostatic coefficients (fi) derived from FEP simulations of some representative compounds in water. For the 18 compound training set that we used this model yielded a mean unsigned error of only 0.58 kcal/mol which seemed very promising. 

Warshel and coworkers have recently examined the LIE method and different versions of what they call the LRA (linear response approximation) method for the binding of a set of cyclic urea compounds to HIV protease.34 The key features of their LRA scheme is that both averages of Equation 2 are evaluated, thus requiring two extra simulations of the non-polar states (see above), that the ligand intramolecular electrostatic terms are included in the averages, and that the non-polar contribution is calculated with the PDLD method. Results of similar quality were reported with the different methods.34 However, it should be noted that the value Vi of the electrostatic coefficient was used in Ref. 34, which, as discussed above, has been shown... 

The interaction of one MTP and one PC molecules not only involves their electrostatic coefficients but also the Lennard-Jones parameters. Achieving a consistent cross-interaction requires care. 

A further complication of the interpretation of contour maps results from the requirement of electroneutrality of a molecule. Electrostatic coefficient contours may be found close to ligand atoms that have polar interactions with specific receptor atoms. However, electrostatic contours, of the opposite sign, may also surround ligand atoms not involved in any polar interaction. The charge distribution pattern is felt by CoMFA from both extremes of the molecular dipole. As a result, both elearon-rich and electron-poor regions will be included in the model. Statistically, this is a simple case of correlation between descriptors (collinearity). 

Stigter and Dill [98] studied phospholipid monolayers at the n-heptane-water interface and were able to treat the second and third virial coefficients (see Eq. XV-1) in terms of electrostatic, including dipole, interactions. At higher film pressures, Pethica and co-workers [99] observed quasi-first-order phase transitions, that is, a much flatter plateau region than shown in Fig. XV-6. 

Some electric properties of molecules are described in section Al.5.2.2 because the coefficients of the powers of Mr turn out to be related to them. The electrostatic, mduction and dispersion energies are considered m turn in section Al.5.2.3, section Al.5.2.4 and section Al.5.2.5, respectively. 

The electrostatic potential generated by a molecule A at a distant point B can be expanded m inverse powers of the distance r between B and the centre of mass (CM) of A. This series is called the multipole expansion because the coefficients can be expressed in temis of the multipole moments of the molecule. With this expansion in hand, it is... 

Wave functions can be visualized as the total electron density, orbital densities, electrostatic potential, atomic densities, or the Laplacian of the electron density. The program computes the data from the basis functions and molecular orbital coefficients. Thus, it does not need a large amount of disk space to store data, but the computation can be time-consuming. Molden can also compute electrostatic charges from the wave function. Several visualization modes are available, including contour plots, three-dimensional isosurfaces, and data slices. 

If classical Coulombic interactions are assumed among point charges for electrostatic interactions between solute and solvent, and the term for the Cl coefficients (C) is omitted, the solvated Eock operator is reduced to Eq. (6). The significance of this definition of the Eock operator from a variational principle is that it enables us to express the analytical first derivative of the free energy with respect to the nuclear coordinate of the solute molecule R ,... 

The non-bonded interaction energy, the van-der-Waals and electrostatic part of the interaction Hamiltonian are best determined by parametrizing a molecular liquid that contains the same chemical groups as the polymers against the experimentally measured thermodynamical and dynamical data, e.g., enthalpy of vaporization, diffusion coefficient, or viscosity. The parameters can then be transferred to polymers, as was done in our case, for instance in polystyrene (from benzene) [19] or poly (vinyl alcohol) (from ethanol) [20,21]. 

We denote by x the distance from the metal surface, and by n x) and rip x) the concentrations of cation vancancies and positive holes in the oxide. Let and Vp be their mobilities, and and Dp their diffusion coefficients. Let F x) be the electrostatic field in the oxide. J, the flux of cation vacancies (number crossing unit area per second), will be expressed by... 

It will be recalled that in Chapter 7, in dealing with the anions of weak acids, we found that we were able to obtain good qualitative agreement between theory and experiment, if we replaced Jmv by the electrostatic Jti. As has been pointed out more than once, the use of simple electrostatic theory necessarily implies that there is an increment in the degree of order in the solvent subject to the electrostatic field. In the -coefficients of Table 24, however, we have found that for 6 out of the... 

Turning next to the unitary part of AS0, this is given in Table 36 under the heading — N(dL/dT). It was pointed out in Secs. 90 and 106 that, to obtain the unitary part of AS0 in aqueous solution, one must subtract 16.0 e.u. for a uni-univalent solute, and 24.0 e.u. for a uni-divalent solute. In methanol solution the corresponding quantities are 14.0 and 21.0 e.u. In Table 36 it will be seen that, except for the first two solutes KBr and KC1, the values are all negative, in both solvents. It will be recalled that for KBr and KC1 the B-coefficients in viscosity are negative, and we associate the positive values for the unitary part of the entropy, shown in Table 29, with the creation of disorder in the ionic co-spheres. In every solvent the dielectric constant decreases with rise of temperature and this leads us to expect that L will increase. For KBr and KC1 in methanol solution, we see from Table 36 that dL/dT has indeed a large positive value. On the other hand, when these crystals dissolve in water, these electrostatic considerations appear to be completely overbalanced by other factors. 

Another problem arises from the presence of higher terms in the multipole expansion of the electrostatic interaction. While theoretical formulas exist for these also, they are even more approximate than those for the dipole-dipole term. Also, there is the uncertainty about the exact form of the repulsive interaction. Quite arbitrarily we shall group the higher multipole terms with the true repulsive interaction and assume that the empirical repulsive term accounts for both. The principal merit of this assumption is simplicity the theoretical and experimental coefficients of the R Q term are compared without adjustment. Since the higher multipole terms are known to be attractive and have been estimated to amount to about 20 per cent of the total attractive potential at the minimum, a rough correction for their possible effect can be made if it is believed that this is a preferable assumption. 

---