A fast method for calculating geometry-dependent net atomic charges for polypeptides

Kwang Hwi Cho, Young Kee Kang, Kyoung Tai No, Harold A. Scheraga

Research output: Contribution to journalArticle

25 Citations (Scopus)

Abstract

A new method for calculating geometry-dependent atomic charges (GDACs) for polypeptides is presented. It overcomes the limitations of the partial equalization of orbital electronegativity (PEOE) and modified PEOE (MPEOE) methods, which depend only on connectivity, not geometry. Introduction of distance-dependent damping factors helps to include the effect of environment in determining the variation of bond distance (without explicit contribution of the correlated variation of bond angles), and thereby to reduce the number of parameters required to represent different atomic species. Since the correlation between the geometry and the dipole moments of molecules is crucial for this method, ab initio molecular orbital calculations were carried out to obtain the geometries and dipole moments with the 6-31G** basis set at the level of B3LYP theory. When bond distances are fixed prior to a charge calculation, the methodology outlined here leads to a direct calculation of the permanent molecular charge distribution represented as a set of distributed monopoles that depend on the geometry of the molecule. Hence, this method automatically accounts for the transferability of charges of small amino acid residues to build up a large polypeptide molecule, and can therefore provide an approximate description of any redistribution of charge density of large polypeptide molecules. The parameters characterizing the charge transfer in the formation of bonds were optimized by using dipole moment components and total dipole moments of 50 molecules that serve as models for the backbone and side chains of proteins. The calculated total dipole moments of these 50 molecules agree well with the ab initio results within an error of 5%. The new charge scheme has been applied to seven conformers of W-acetylalanineN′-methylamide (Ac-Ala-NHMe) with good agreement between ab initio and GDAC dipole moments. This method, however, gives poor results for conjugated systems that are larger than amides.

Original languageEnglish
Pages (from-to)3624-3634
Number of pages11
JournalJournal of Physical Chemistry B
Volume105
Issue number17
Publication statusPublished - 2001 May 3

Fingerprint

Polypeptides
polypeptides
Dipole moment
dipole moments
Peptides
Molecules
Geometry
geometry
molecules
Electronegativity
orbitals
Orbital calculations
Charge distribution
Molecular orbitals
Charge density
Amides
monopoles
charge distribution
amides
amino acids

All Science Journal Classification (ASJC) codes

  • Physical and Theoretical Chemistry
  • Surfaces, Coatings and Films
  • Materials Chemistry

Cite this

Cho, Kwang Hwi ; Kang, Young Kee ; No, Kyoung Tai ; Scheraga, Harold A. / A fast method for calculating geometry-dependent net atomic charges for polypeptides. In: Journal of Physical Chemistry B. 2001 ; Vol. 105, No. 17. pp. 3624-3634.
@article{22889280622e4a6b995b6033c90ee698,
title = "A fast method for calculating geometry-dependent net atomic charges for polypeptides",
abstract = "A new method for calculating geometry-dependent atomic charges (GDACs) for polypeptides is presented. It overcomes the limitations of the partial equalization of orbital electronegativity (PEOE) and modified PEOE (MPEOE) methods, which depend only on connectivity, not geometry. Introduction of distance-dependent damping factors helps to include the effect of environment in determining the variation of bond distance (without explicit contribution of the correlated variation of bond angles), and thereby to reduce the number of parameters required to represent different atomic species. Since the correlation between the geometry and the dipole moments of molecules is crucial for this method, ab initio molecular orbital calculations were carried out to obtain the geometries and dipole moments with the 6-31G** basis set at the level of B3LYP theory. When bond distances are fixed prior to a charge calculation, the methodology outlined here leads to a direct calculation of the permanent molecular charge distribution represented as a set of distributed monopoles that depend on the geometry of the molecule. Hence, this method automatically accounts for the transferability of charges of small amino acid residues to build up a large polypeptide molecule, and can therefore provide an approximate description of any redistribution of charge density of large polypeptide molecules. The parameters characterizing the charge transfer in the formation of bonds were optimized by using dipole moment components and total dipole moments of 50 molecules that serve as models for the backbone and side chains of proteins. The calculated total dipole moments of these 50 molecules agree well with the ab initio results within an error of 5{\%}. The new charge scheme has been applied to seven conformers of W-acetylalanineN′-methylamide (Ac-Ala-NHMe) with good agreement between ab initio and GDAC dipole moments. This method, however, gives poor results for conjugated systems that are larger than amides.",
author = "Cho, {Kwang Hwi} and Kang, {Young Kee} and No, {Kyoung Tai} and Scheraga, {Harold A.}",
year = "2001",
month = "5",
day = "3",
language = "English",
volume = "105",
pages = "3624--3634",
journal = "Journal of Physical Chemistry B Materials",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "17",

}

A fast method for calculating geometry-dependent net atomic charges for polypeptides. / Cho, Kwang Hwi; Kang, Young Kee; No, Kyoung Tai; Scheraga, Harold A.

In: Journal of Physical Chemistry B, Vol. 105, No. 17, 03.05.2001, p. 3624-3634.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A fast method for calculating geometry-dependent net atomic charges for polypeptides

AU - Cho, Kwang Hwi

AU - Kang, Young Kee

AU - No, Kyoung Tai

AU - Scheraga, Harold A.

PY - 2001/5/3

Y1 - 2001/5/3

N2 - A new method for calculating geometry-dependent atomic charges (GDACs) for polypeptides is presented. It overcomes the limitations of the partial equalization of orbital electronegativity (PEOE) and modified PEOE (MPEOE) methods, which depend only on connectivity, not geometry. Introduction of distance-dependent damping factors helps to include the effect of environment in determining the variation of bond distance (without explicit contribution of the correlated variation of bond angles), and thereby to reduce the number of parameters required to represent different atomic species. Since the correlation between the geometry and the dipole moments of molecules is crucial for this method, ab initio molecular orbital calculations were carried out to obtain the geometries and dipole moments with the 6-31G** basis set at the level of B3LYP theory. When bond distances are fixed prior to a charge calculation, the methodology outlined here leads to a direct calculation of the permanent molecular charge distribution represented as a set of distributed monopoles that depend on the geometry of the molecule. Hence, this method automatically accounts for the transferability of charges of small amino acid residues to build up a large polypeptide molecule, and can therefore provide an approximate description of any redistribution of charge density of large polypeptide molecules. The parameters characterizing the charge transfer in the formation of bonds were optimized by using dipole moment components and total dipole moments of 50 molecules that serve as models for the backbone and side chains of proteins. The calculated total dipole moments of these 50 molecules agree well with the ab initio results within an error of 5%. The new charge scheme has been applied to seven conformers of W-acetylalanineN′-methylamide (Ac-Ala-NHMe) with good agreement between ab initio and GDAC dipole moments. This method, however, gives poor results for conjugated systems that are larger than amides.

AB - A new method for calculating geometry-dependent atomic charges (GDACs) for polypeptides is presented. It overcomes the limitations of the partial equalization of orbital electronegativity (PEOE) and modified PEOE (MPEOE) methods, which depend only on connectivity, not geometry. Introduction of distance-dependent damping factors helps to include the effect of environment in determining the variation of bond distance (without explicit contribution of the correlated variation of bond angles), and thereby to reduce the number of parameters required to represent different atomic species. Since the correlation between the geometry and the dipole moments of molecules is crucial for this method, ab initio molecular orbital calculations were carried out to obtain the geometries and dipole moments with the 6-31G** basis set at the level of B3LYP theory. When bond distances are fixed prior to a charge calculation, the methodology outlined here leads to a direct calculation of the permanent molecular charge distribution represented as a set of distributed monopoles that depend on the geometry of the molecule. Hence, this method automatically accounts for the transferability of charges of small amino acid residues to build up a large polypeptide molecule, and can therefore provide an approximate description of any redistribution of charge density of large polypeptide molecules. The parameters characterizing the charge transfer in the formation of bonds were optimized by using dipole moment components and total dipole moments of 50 molecules that serve as models for the backbone and side chains of proteins. The calculated total dipole moments of these 50 molecules agree well with the ab initio results within an error of 5%. The new charge scheme has been applied to seven conformers of W-acetylalanineN′-methylamide (Ac-Ala-NHMe) with good agreement between ab initio and GDAC dipole moments. This method, however, gives poor results for conjugated systems that are larger than amides.

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

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

M3 - Article

AN - SCOPUS:0037013378

VL - 105

SP - 3624

EP - 3634

JO - Journal of Physical Chemistry B Materials

JF - Journal of Physical Chemistry B Materials

SN - 1520-6106

IS - 17

ER -