Implementation of π‐π interactions in molecular dynamics simulation

No explicit π‐π interaction term has been incorporated in the conventional molecular dynamics (MD) simulation programs in spite of its significant role in the folding of biomolecules and the clustering of organic chemicals. In this article, we propose a technique to emphasize the effect of π‐π interactions using a function of energy and implement it into an MD simulation program. Several trial calculations show that the π‐π incorporated program gives improved results consistent with experimental data on atom geometry and has no unfavorable interference with the conventional computational framework. This indicates an importance of the explicit consideration of π‐π interactions in MD simulation. © 2007 Wiley Periodicals, Inc.J Comput Chem 2007

[1]  J. Thornton,et al.  Amino/aromatic interactions in proteins: is the evidence stacked against hydrogen bonding? , 1994, Journal of molecular biology.

[2]  C. David Sherrill,et al.  Substituent Effects in π−π Interactions: Sandwich and T-Shaped Configurations , 2004 .

[3]  C. David Sherrill,et al.  Highly Accurate Coupled Cluster Potential Energy Curves for the Benzene Dimer: Sandwich, T-Shaped, and Parallel-Displaced Configurations , 2004 .

[4]  MD simulation of asymmetric phospholipid bilayers with ions and cholesterols , 2004 .

[5]  Jose Cosme,et al.  Crystal structure of human cytochrome P450 2C9 with bound warfarin , 2003, Nature.

[6]  Paul T. Matsudaira,et al.  NMR structure of the 35-residue villin headpiece subdomain , 1997, Nature Structural Biology.

[7]  N. Skelton,et al.  Tryptophan zippers: Stable, monomeric β-hairpins , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[8]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[9]  S. Tsuzuki,et al.  Origin of attraction and directionality of the pi/pi interaction: model chemistry calculations of benzene dimer interaction. , 2002, Journal of the American Chemical Society.

[10]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[11]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

[12]  J. Sussman,et al.  Structure of acetylcholinesterase complexed with E2020 (Aricept): implications for the design of new anti-Alzheimer drugs. , 1999, Structure.

[13]  Kazumasa Honda,et al.  Origin of the Attraction and Directionality of the NH/π Interaction: Comparison with OH/π and CH/π Interactions , 2000 .

[14]  Masayuki Hata,et al.  Resistant mechanism against nelfinavir of human immunodeficiency virus type 1 proteases. , 2005 .

[15]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[16]  Eric F. Johnson,et al.  The Structure of Human Microsomal Cytochrome P450 3A4 Determined by X-ray Crystallography to 2.05-Å Resolution* , 2004, Journal of Biological Chemistry.

[17]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[18]  C. Hunter,et al.  Sequence-dependent DNA structure. The role of base stacking interactions. , 1993, Journal of molecular biology.

[19]  Jose Cosme,et al.  Crystal Structures of Human Cytochrome P450 3A4 Bound to Metyrapone and Progesterone , 2004, Science.

[20]  Y. Shimohigashi,et al.  X-ray crystal structure of a dipeptide-chymotrypsin complex in an inhibitory interaction. , 1998, European journal of biochemistry.

[21]  U. Samanta,et al.  Environment of tryptophan side chains in proteins , 2000, Proteins.

[22]  Alexander D. MacKerell,et al.  CH/π interactions involving aromatic amino acids: Refinement of the CHARMM tryptophan force field , 2005, J. Comput. Chem..

[23]  C. Ampe,et al.  Solution structures of the C‐terminal headpiece subdomains of human villin and advillin, evaluation of headpiece F‐actin‐binding requirements , 2004, Protein science : a publication of the Protein Society.

[24]  D. Hung,et al.  Transition State π‐Solvation by Aromatic Rings: An Electronic Contribution to Diels–Alder Reaction Diastereoselectivity , 1987 .

[25]  M. Hata,et al.  Direct interaction between substrates and endogenous steroids in the active site may change the activity of cytochrome P450 3A4. , 2003, Biochemistry.

[26]  G A Petsko,et al.  Aromatic-aromatic interaction: a mechanism of protein structure stabilization. , 1985, Science.

[27]  A. Suenaga,et al.  Computational studies on prion proteins: effect of Ala(117)-->Val mutation. , 2002, Biophysical journal.

[28]  S. Neya,et al.  Computational analysis of the transient movement of helices in sensory rhodopsin II , 2005, Protein science : a publication of the Protein Society.

[29]  M. Machius,et al.  Pivotal role of water in the mechanism of P450BM-3. , 2001, Biochemistry.