Intraprotein Electrostatics Derived from First Principles: Divide-and-Conquer Approaches for QM/MM Calculations

Pablo A. Molina, Hui Li, Jan H. Jensen

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of Nζ of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pKa of Lys55, and compared to other computational methodologies in which force field charges are employed.

Original languageEnglish (US)
Pages (from-to)1971-1979
Number of pages9
JournalJournal of Computational Chemistry
Volume24
Issue number16
DOIs
StatePublished - Dec 1 2003

Fingerprint

Divide and conquer
First-principles
Electrostatics
Charge
Force Field
Renormalization
Inhibitor
Affine transformation
Carrier concentration
Protons
Electron
Proteins
Protein
Integer
Methodology
DAQ

ASJC Scopus subject areas

  • Chemistry(all)
  • Computational Mathematics

Cite this

Intraprotein Electrostatics Derived from First Principles : Divide-and-Conquer Approaches for QM/MM Calculations. / Molina, Pablo A.; Li, Hui; Jensen, Jan H.

In: Journal of Computational Chemistry, Vol. 24, No. 16, 01.12.2003, p. 1971-1979.

Research output: Contribution to journalArticle

@article{5620bcc49d0643519a2d81a9ce26617b,
title = "Intraprotein Electrostatics Derived from First Principles: Divide-and-Conquer Approaches for QM/MM Calculations",
abstract = "Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of Nζ of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pKa of Lys55, and compared to other computational methodologies in which force field charges are employed.",
author = "Molina, {Pablo A.} and Hui Li and Jensen, {Jan H.}",
year = "2003",
month = "12",
day = "1",
doi = "10.1002/jcc.10333",
language = "English (US)",
volume = "24",
pages = "1971--1979",
journal = "Journal of Computational Chemistry",
issn = "0192-8651",
publisher = "John Wiley and Sons Inc.",
number = "16",

}

TY - JOUR

T1 - Intraprotein Electrostatics Derived from First Principles

T2 - Divide-and-Conquer Approaches for QM/MM Calculations

AU - Molina, Pablo A.

AU - Li, Hui

AU - Jensen, Jan H.

PY - 2003/12/1

Y1 - 2003/12/1

N2 - Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of Nζ of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pKa of Lys55, and compared to other computational methodologies in which force field charges are employed.

AB - Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of Nζ of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pKa of Lys55, and compared to other computational methodologies in which force field charges are employed.

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

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

U2 - 10.1002/jcc.10333

DO - 10.1002/jcc.10333

M3 - Article

C2 - 14531051

AN - SCOPUS:0242494830

VL - 24

SP - 1971

EP - 1979

JO - Journal of Computational Chemistry

JF - Journal of Computational Chemistry

SN - 0192-8651

IS - 16

ER -