Question: Implementation Of Π-Π And Aromatic Interactions In Protein-Ligand Interactions
gravatar for Flow
8.2 years ago by
Flow1.5k wrote:

Do you know about any docking software, or any other approach, that can take into account efficiently π-π, aromatic and cation-π interactions in protein-ligand interactions?

PS: could you give bibliographic references?

interaction scoring docking • 2.4k views
ADD COMMENTlink modified 8.2 years ago by João Rodrigues2.5k • written 8.2 years ago by Flow1.5k

you could also use "pure" quantum methods such Gaussian, GAMESS and the like

ADD REPLYlink written 8.2 years ago by Hans Laxness0
gravatar for Vladimir Chupakhin
8.2 years ago by
Toledo, Spain
Vladimir Chupakhin520 wrote:

Mostly all of the docking software include pi-pi interactions but treat them as vdw interations. If you want to include true pi-X use more sophisticated scoring functions for example QM-Polarized Ligand Docking from Schrodinger - but this one should be time consuming.

ADD COMMENTlink written 8.2 years ago by Vladimir Chupakhin520

Some standard molecular modelling programs (GROMACS) support QM/MM methods so that would be probably worth looking into. These usually interface standard QM software such as GAMESS, MOPAC, GAUSSIAN, the first two available for free.

ADD REPLYlink written 8.2 years ago by João Rodrigues2.5k

I knew about Maestro, the thing is, do you know about any open source option?

ADD REPLYlink written 8.2 years ago by Flow1.5k
gravatar for João Rodrigues
8.2 years ago by
João Rodrigues2.5k
Stanford University, U
João Rodrigues2.5k wrote:

This is by no means an extensive analysis of the subject, but I was curious to know this as well and I googled a bit.

This message about the AMBERff seems to be a good starting point. If you can access the forcefield parameters, you can see that most ring carbons have a slightly negative charge indeed. I also checked the forcefield that we use in HADDOCK and the atom type of ring carbons in tryptophan, despite having neutral charge, is that of a carbonyl carbon. In general, I would say this is represented by the ring having a slightly negative charge..

Nevertheless, from what I read, most potentials do not include this effect explicitly, and might fail at identifying it, or at least giving it its proper strength. As chupvl suggested, you will have to go to the QM level to accurately describe this. I found a paper comparing bit this effect in current forcefields, here.

The CHARMM forums might give more insight, but they are currently inaccessible..

ADD COMMENTlink written 8.2 years ago by João Rodrigues2.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1802 users visited in the last hour