Previous Thread
Next Thread
Print Thread
Can Charmm compute pi stacking interaction?
#38229 03/25/21 12:17 PM
Joined: Nov 2019
Posts: 51
A
Forum Member
OP Online Content
Forum Member
A
Joined: Nov 2019
Posts: 51
For a system of protein containing phenylalanine residue, I used following command to calculate pi-pi interaction
SKIP ALL EXCL vdw Elec
interaction select resname phe .and. type C* end sele resname phe .and. type C* end
ECHO "PHEC PHEC" @X ?VDW ?ELEC

I expected to see at least one of the two terms to be negative. But both VDW and ELEC are positive.

Re: Can Charmm compute pi stacking interaction?
Allen_123 #38230 03/25/21 02:34 PM
Joined: Sep 2003
Posts: 4,817
Likes: 2
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 4,817
Likes: 2
You select only the carbon atoms, so the electrostatic interaction will come out as repulsive since all the ring-carbons have the same charge. Try with a selection that includes all the ring atoms, and no backbone atoms. Something like: sele resname phe .and. (chem ca .or. chem hp) end. You may also want to look at specific phe-phe pairs. If you want to play with it, you can easily create a system with just two phenylalanine rings and move them apart while evaluating the energies every 0.1Å or so; perhaps also with the rings in different orientations


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
Re: Can Charmm compute pi stacking interaction?
Allen_123 #38235 03/26/21 05:43 AM
Joined: Nov 2019
Posts: 51
A
Forum Member
OP Online Content
Forum Member
A
Joined: Nov 2019
Posts: 51
I selected all the carbons and hydrogens of ring atoms from PHE residue as suggested. Still, I got positive VDW and ELEC term. Since charges of CA and HP are opposite, should I simply compute the interaction between CA and HP atoms instead ?(ex: interaction select resname phe .and. chem CA end sele resname phe .and. chem HP end).

Afterwards, I did get negative VDW and ELEC by computing interaction between CA and HP. And I also computed interaction between CA and CA and between HP and HP. But the sum of positive interaction is larger than the negative CA -HP. This result still does not make sense to me since in reality the Pi stacking interaction takes account of every atom of PHE ring.

Last edited by Allen_123; 03/26/21 05:57 AM.
Re: Can Charmm compute pi stacking interaction?
Allen_123 #38236 03/26/21 06:05 AM
Joined: Nov 2019
Posts: 51
A
Forum Member
OP Online Content
Forum Member
A
Joined: Nov 2019
Posts: 51
From VMD, the PHE rings of crystal should be close enough.

Attached Files
PROcrystal.JPG (38.12 KB, 16 downloads)
Last edited by Allen_123; 03/26/21 06:08 AM.
Re: Can Charmm compute pi stacking interaction?
Allen_123 #38238 03/26/21 03:25 PM
Joined: Sep 2003
Posts: 4,817
Likes: 2
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 4,817
Likes: 2
OK. The relevant aspect of the energies computed by CHARMM is that they should always be used in a relative sense, so you should always compute the difference wrt some reference state. This means that even though the energy values you compute are all positive, there may still be an attraction between the two rings, as long as there is a minimum in the energy surface. If I remember correctly the minimum conformation for two benzenes using the CHARMM energy function is where the edge of one interacts with the flat surface of the other (more or less).

There are also interactions within each ring, which are probably positive, but not important for the interaction between the rings. To avoid this complication you can use the interaction_energy command (energy.info).

Last edited by lennart; 03/26/21 03:59 PM. Reason: added interaction energy

Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden

Moderated by  BRBrooks, lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.4
(Release build 20200307)
Responsive Width:

PHP: 5.6.33-0+deb8u1 Page Time: 0.011s Queries: 25 (0.005s) Memory: 0.9247 MB (Peak: 1.0210 MB) Data Comp: Off Server Time: 2021-04-20 22:35:56 UTC
Valid HTML 5 and Valid CSS