Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water

Comparisons with IR and vibrational circular dichroism spectra

Kijeong Kwac, Kyung Koo Lee, Jae Bum Han, Kwang Im Oh, Minhaeng Cho

Research output: Contribution to journalArticle

62 Citations (Scopus)

Abstract

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed α -helix and extended Β and PII conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the PII conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is αR. Polarizable CHARMM force field gives almost exclusively PII conformation and other force fields predict that both α -helical and extended (Β and PII) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the PII conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.

Original languageEnglish
Article number105106
JournalJournal of Chemical Physics
Volume128
Issue number10
DOIs
Publication statusPublished - 2008 Mar 25

Fingerprint

Dipeptides
alanine
Alanine
dichroism
Conformations
Molecular dynamics
molecular dynamics
Water
Computer simulation
field theory (physics)
water
simulation
solutes
Solvation
Dipole moment
Distribution functions
Circular Dichroism
radial distribution
Hydrogen bonds
helices

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water : Comparisons with IR and vibrational circular dichroism spectra. / Kwac, Kijeong; Lee, Kyung Koo; Han, Jae Bum; Oh, Kwang Im; Cho, Minhaeng.

In: Journal of Chemical Physics, Vol. 128, No. 10, 105106, 25.03.2008.

Research output: Contribution to journalArticle

@article{ea56789789f44c9287417be0c9263ab0,
title = "Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: Comparisons with IR and vibrational circular dichroism spectra",
abstract = "We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed α -helix and extended Β and PII conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the PII conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is αR. Polarizable CHARMM force field gives almost exclusively PII conformation and other force fields predict that both α -helical and extended (Β and PII) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the PII conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.",
author = "Kijeong Kwac and Lee, {Kyung Koo} and Han, {Jae Bum} and Oh, {Kwang Im} and Minhaeng Cho",
year = "2008",
month = "3",
day = "25",
doi = "10.1063/1.2837461",
language = "English",
volume = "128",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "10",

}

TY - JOUR

T1 - Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water

T2 - Comparisons with IR and vibrational circular dichroism spectra

AU - Kwac, Kijeong

AU - Lee, Kyung Koo

AU - Han, Jae Bum

AU - Oh, Kwang Im

AU - Cho, Minhaeng

PY - 2008/3/25

Y1 - 2008/3/25

N2 - We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed α -helix and extended Β and PII conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the PII conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is αR. Polarizable CHARMM force field gives almost exclusively PII conformation and other force fields predict that both α -helical and extended (Β and PII) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the PII conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.

AB - We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed α -helix and extended Β and PII conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the PII conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is αR. Polarizable CHARMM force field gives almost exclusively PII conformation and other force fields predict that both α -helical and extended (Β and PII) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the PII conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.

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

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

U2 - 10.1063/1.2837461

DO - 10.1063/1.2837461

M3 - Article

VL - 128

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 10

M1 - 105106

ER -