No abstract

An implicit solvent model ISM for use in molecular mechanics, MonteCarlo, and molecular dynamics simulations on proteins and nucleic acids is proposed that is based on describing the electrostatic component with a screened Coulomb Ž . potential SCP . The SCP has been extended so that both the electrostatic interaction energies and the self-energy terms are included in the solvation energy where the latter has been calculated from the integral Born equation. In addition, the SCP is generalized to allow an accounting of the positional dependence of the interactions between charges in a dielectric. To test the potential of the method to provide a reliable and fast implicit description of solvent, a parameter set has been determined to calculate hydration free energies of a group of small-molecule analogs for neutral amino acid residues. Comparison of the calculated with experimental hydration energies of 18 molecules yields a correlation greater than 0.99, demonstrating the feasibility of the method to form the basis of an ISM. A parameter set based solely on hydration energies is not sufficient to account for the considerably different physical conditions found for a solute in the protein ''solvent.'' Methods are explored for further generalizing the parameters to account for macromolecular structure, and it is shown that it may be possible to find a positional function of the coordinates that correlates well with the fraction of each amino acid residue buried in the macromolecule. Such a function would reduce computing time by replacing the need for repetitive calculations of the solvent-accessible surface area Ž . and its derivatives in the case of molecular dynamics simulations with a simple

A continuum hydration model based upon the atomic charges provided with the CFF91 force field [A. B. Schmidt and R. M. Fine (1994) Molecular Simulation, 13, 347-365] has been extended to the octanol-water transfer. The electrostatic component of the transfer free energy is calculated using the finite-difference solution to the Poisson-Boltzmann equation while the nonpolar contributions are assumed to be proportional to the solute-excluded volume in water. All atomic charges and radii besides the aromatic carbon radius are equal in both solvents. The octanol dielectric constant and the probe radius are the main fitting parameters defining the octanol phase. The model has been tested for 38 organic molecules related to the amino acid residues and generally provides a high accuracy. In particular, the mean unsigned error for N-acetyl amino acid amides is 0.5 kcal/mol.

No abstract

No abstract