File Download
There are no files associated with this item.
Links for fulltext
(May Require Subscription)
- Publisher Website: 10.1002/prot.10279
- Scopus: eid_2-s2.0-0037441479
- PMID: 12557187
- WOS: WOS:000180691300007
- Find via
Supplementary
-
Bookmarks:
- CiteULike: 7
- Citations:
- Appears in Collections:
Article: Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine "dipeptides" (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution
Title | Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine "dipeptides" (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution |
---|---|
Authors | |
Keywords | Alanine dipeptide Dynamics simulation Glycine dipeptide Peptide backbone QM/MM model Ramachandran plot Solution conformation |
Issue Date | 2003 |
Citation | Proteins: Structure, Function And Genetics, 2003, v. 50 n. 3, p. 451-463 How to Cite? |
Abstract | We compare the conformational distributions of Ace-Ala-Nme and Ace-Gly-Nme sampled in long simulations with several molecular mechanics (MM) force fields and with a fast combined quantum mechanics/molecular mechanics (QM/MM) force field, in which the solute's intramolecular energy and forces are calculated with the self-consistent charge density functional tight binding method (SCCDFTB), and the solvent is represented by either one of the well-known SPC and TIP3P models. All MM force fields give two main states for Ace-Ala-Nme, β and α separated by free energy barriers, but the ratio in which these are sampled varies by a factor of 30, from a high in favor of β of 6 to a low of 1/5. The frequency of transitions between states is particularly low with the amber and charmm force fields, for which the distributions are noticeably narrower, and the energy barriers between states higher. The lower of the two barriers lies between α and β at values of ψ near 0 for all MM simulations except for charmm22. The results of the QM/MM simulations vary less with the choice of MM force field; the ratio β/α varies between 1.5 and 2.2, the easy pass lies at ψ near 0, and transitions between states are more frequent than for amber and charmm, but less frequent than for cedar. For Ace-Gly-Nme, all force fields locate a diffuse stable region around φ = π and ψ = 0, whereas the amber force field gives two additional densely sampled states near φ = ±100° and ψ = 0, which are also found with the QM/MM force field. For both solutes, the distribution from the QM/MM simulation shows greater similarity with the distribution in highresolution protein structures than is the case for any of the MM simulations. © 2003 Wiley-Liss, Inc. |
Persistent Identifier | http://hdl.handle.net/10722/167783 |
ISSN | 2023 Impact Factor: 3.2 2023 SCImago Journal Rankings: 1.086 |
ISI Accession Number ID | |
References |
DC Field | Value | Language |
---|---|---|
dc.contributor.author | Hu, H | en_US |
dc.contributor.author | Elstner, M | en_US |
dc.contributor.author | Hermans, J | en_US |
dc.date.accessioned | 2012-10-08T03:11:31Z | - |
dc.date.available | 2012-10-08T03:11:31Z | - |
dc.date.issued | 2003 | en_US |
dc.identifier.citation | Proteins: Structure, Function And Genetics, 2003, v. 50 n. 3, p. 451-463 | en_US |
dc.identifier.issn | 0887-3585 | en_US |
dc.identifier.uri | http://hdl.handle.net/10722/167783 | - |
dc.description.abstract | We compare the conformational distributions of Ace-Ala-Nme and Ace-Gly-Nme sampled in long simulations with several molecular mechanics (MM) force fields and with a fast combined quantum mechanics/molecular mechanics (QM/MM) force field, in which the solute's intramolecular energy and forces are calculated with the self-consistent charge density functional tight binding method (SCCDFTB), and the solvent is represented by either one of the well-known SPC and TIP3P models. All MM force fields give two main states for Ace-Ala-Nme, β and α separated by free energy barriers, but the ratio in which these are sampled varies by a factor of 30, from a high in favor of β of 6 to a low of 1/5. The frequency of transitions between states is particularly low with the amber and charmm force fields, for which the distributions are noticeably narrower, and the energy barriers between states higher. The lower of the two barriers lies between α and β at values of ψ near 0 for all MM simulations except for charmm22. The results of the QM/MM simulations vary less with the choice of MM force field; the ratio β/α varies between 1.5 and 2.2, the easy pass lies at ψ near 0, and transitions between states are more frequent than for amber and charmm, but less frequent than for cedar. For Ace-Gly-Nme, all force fields locate a diffuse stable region around φ = π and ψ = 0, whereas the amber force field gives two additional densely sampled states near φ = ±100° and ψ = 0, which are also found with the QM/MM force field. For both solutes, the distribution from the QM/MM simulation shows greater similarity with the distribution in highresolution protein structures than is the case for any of the MM simulations. © 2003 Wiley-Liss, Inc. | en_US |
dc.language | eng | en_US |
dc.relation.ispartof | Proteins: Structure, Function and Genetics | en_US |
dc.subject | Alanine dipeptide | - |
dc.subject | Dynamics simulation | - |
dc.subject | Glycine dipeptide | - |
dc.subject | Peptide backbone | - |
dc.subject | QM/MM model | - |
dc.subject | Ramachandran plot | - |
dc.subject | Solution conformation | - |
dc.subject.mesh | Alanine - Chemistry | en_US |
dc.subject.mesh | Computer Simulation | en_US |
dc.subject.mesh | Databases, Protein | en_US |
dc.subject.mesh | Dipeptides - Chemistry | en_US |
dc.subject.mesh | Glycine - Chemistry | en_US |
dc.subject.mesh | Models, Molecular | en_US |
dc.subject.mesh | Protein Conformation | en_US |
dc.subject.mesh | Protein Folding | en_US |
dc.subject.mesh | Proteins - Chemistry | en_US |
dc.subject.mesh | Quantum Theory | en_US |
dc.subject.mesh | Solutions - Chemistry | en_US |
dc.subject.mesh | Water - Chemistry | en_US |
dc.title | Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine "dipeptides" (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution | en_US |
dc.type | Article | en_US |
dc.identifier.email | Hu, H:haohu@hku.hk | en_US |
dc.identifier.authority | Hu, H=rp00707 | en_US |
dc.description.nature | link_to_subscribed_fulltext | en_US |
dc.identifier.doi | 10.1002/prot.10279 | en_US |
dc.identifier.pmid | 12557187 | - |
dc.identifier.scopus | eid_2-s2.0-0037441479 | en_US |
dc.relation.references | http://www.scopus.com/mlt/select.url?eid=2-s2.0-0037441479&selection=ref&src=s&origin=recordpage | en_US |
dc.identifier.volume | 50 | en_US |
dc.identifier.issue | 3 | en_US |
dc.identifier.spage | 451 | en_US |
dc.identifier.epage | 463 | en_US |
dc.identifier.isi | WOS:000180691300007 | - |
dc.identifier.scopusauthorid | Hu, H=7404097564 | en_US |
dc.identifier.scopusauthorid | Elstner, M=7004554960 | en_US |
dc.identifier.scopusauthorid | Hermans, J=7201896483 | en_US |
dc.identifier.citeulike | 1007761 | - |
dc.identifier.issnl | 0887-3585 | - |