### Abstract

The three-dimensional reference interaction site model integral equation theory (3D-RISM) combined with the ab initio molecular orbital method (3D-RISM-SCF) is applied to a solvated macromolecular system. The solvation structure around a solute molecule is obtained from the 3D-RISM integral equation under the electrostatic potential of the solute molecule, calculated by the ab initio molecular orbital theory. The electrostatic potential should be calculated on each grid point in the three-dimensional real space. Therefore, the calculation of the electrostatic potential is the most time consuming part in this method. In this article, we propose a new procedure to save the computational cost for calculating the electrostatic potential and the solvated fock matrix. The strategy of this procedure is to evaluate the electrostatic potential and the solvated fock matrix in different ways, depending on the distance between solute and solvent. Inside the repulsive cores of solute atoms, it is possible to avoid the calculation of electrostatic potential and solvated Fock matrix by assuming the potential to be infinity. In the region sufficiently far from solute, they are evaluated classically by putting the effective point charge on each atom. In the intermediate region, the electrostatic potential is evaluated directly by integrating the molecular orbitals of the solute molecule. The electronic structure and the energy gradient of Methionine-Enkephalin and solvation structure are estimated by using this procedure in aqueous solution, and are compared with the results from other procedures. The results are compared also with those from the continuum model.

Original language | English |
---|---|

Pages (from-to) | 453-462 |

Number of pages | 10 |

Journal | Journal of Computational Chemistry |

Volume | 27 |

Issue number | 4 |

DOIs | |

Publication status | Published - Mar 1 2006 |

Externally published | Yes |

### Fingerprint

### All Science Journal Classification (ASJC) codes

- Chemistry(all)
- Computational Mathematics

### Cite this

**A new method to determine electrostatic potential around a macromolecule in solution from molecular wave functions.** / Yoshida, Norio; Hirata, Fumio.

Research output: Contribution to journal › Article

*Journal of Computational Chemistry*, vol. 27, no. 4, pp. 453-462. https://doi.org/10.1002/jcc.20356

}

TY - JOUR

T1 - A new method to determine electrostatic potential around a macromolecule in solution from molecular wave functions

AU - Yoshida, Norio

AU - Hirata, Fumio

PY - 2006/3/1

Y1 - 2006/3/1

N2 - The three-dimensional reference interaction site model integral equation theory (3D-RISM) combined with the ab initio molecular orbital method (3D-RISM-SCF) is applied to a solvated macromolecular system. The solvation structure around a solute molecule is obtained from the 3D-RISM integral equation under the electrostatic potential of the solute molecule, calculated by the ab initio molecular orbital theory. The electrostatic potential should be calculated on each grid point in the three-dimensional real space. Therefore, the calculation of the electrostatic potential is the most time consuming part in this method. In this article, we propose a new procedure to save the computational cost for calculating the electrostatic potential and the solvated fock matrix. The strategy of this procedure is to evaluate the electrostatic potential and the solvated fock matrix in different ways, depending on the distance between solute and solvent. Inside the repulsive cores of solute atoms, it is possible to avoid the calculation of electrostatic potential and solvated Fock matrix by assuming the potential to be infinity. In the region sufficiently far from solute, they are evaluated classically by putting the effective point charge on each atom. In the intermediate region, the electrostatic potential is evaluated directly by integrating the molecular orbitals of the solute molecule. The electronic structure and the energy gradient of Methionine-Enkephalin and solvation structure are estimated by using this procedure in aqueous solution, and are compared with the results from other procedures. The results are compared also with those from the continuum model.

AB - The three-dimensional reference interaction site model integral equation theory (3D-RISM) combined with the ab initio molecular orbital method (3D-RISM-SCF) is applied to a solvated macromolecular system. The solvation structure around a solute molecule is obtained from the 3D-RISM integral equation under the electrostatic potential of the solute molecule, calculated by the ab initio molecular orbital theory. The electrostatic potential should be calculated on each grid point in the three-dimensional real space. Therefore, the calculation of the electrostatic potential is the most time consuming part in this method. In this article, we propose a new procedure to save the computational cost for calculating the electrostatic potential and the solvated fock matrix. The strategy of this procedure is to evaluate the electrostatic potential and the solvated fock matrix in different ways, depending on the distance between solute and solvent. Inside the repulsive cores of solute atoms, it is possible to avoid the calculation of electrostatic potential and solvated Fock matrix by assuming the potential to be infinity. In the region sufficiently far from solute, they are evaluated classically by putting the effective point charge on each atom. In the intermediate region, the electrostatic potential is evaluated directly by integrating the molecular orbitals of the solute molecule. The electronic structure and the energy gradient of Methionine-Enkephalin and solvation structure are estimated by using this procedure in aqueous solution, and are compared with the results from other procedures. The results are compared also with those from the continuum model.

UR - http://www.scopus.com/inward/record.url?scp=33644880787&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33644880787&partnerID=8YFLogxK

U2 - 10.1002/jcc.20356

DO - 10.1002/jcc.20356

M3 - Article

VL - 27

SP - 453

EP - 462

JO - Journal of Computational Chemistry

JF - Journal of Computational Chemistry

SN - 0192-8651

IS - 4

ER -