Dear Franz,
Thanks a lot for your answer. In fact, in view of the data values you reported to be normal in your case, I double-checked everything and I realized that the high values I was obtainted were not for "Esw" but for "Eww" (both absolute and normalized).
I've been trying different GIST options and in fact you can obtain extremely high values for the Esw, but just if you do not remove the ions of your trajectory (as it is actually mentioned in the manual). Hence, if ions are not stripped from the trajectory file before running GIST through cpptraj, they are considered as part of the solute and you get that crazy numbers.
However, the question of whether the very high numbers I get for Eww are expected or not is still open...
Thanks again,
Best,
Reyes
-----Mensaje original-----
De: Franz Waibl <franz.waibl.uibk.ac.at>
Enviado el: 09 September 2021 10:46
Para: amber.ambermd.org
Asunto: Re: [AMBER] GIST question
CAUTION: This email originated from outside the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe.
Dear Reyes,
An Esw-norm value of 125576000 is pretty strange, especially if it is in a high-density region. For instance, in one of my GIST calculations, the Esw values ranged from -27.7 to 10.8. You should double check whether this value is correct, and whether you get the same value with a different version of the GIST code (i.e.: different cpptraj version, or the PME-GIST or GIGIST implementations). Which force field / water model are you using?
Regarding the values for consecutive voxels, it is expected that they can have different energy values. The default grid spacing for GIST is
0.5 Å, and 0.5 Å can have a large impact on the water energy, especially at the protein surface.
Sorry if this doesn't really answer your questions. I figured I'd reply anyway since the GIST user-base is rather small. I have worked with GIST a lot, but mainly used ff14SB and TIP3P, and didn't ever notice such high Esw values.
Best regards,
Franz
On 08.09.21 18:24, Reyes Nuñez Franco wrote:
> Dear Amber users,
>
>
> I am using GIST to analyze the structure and thermodynamics of water in the vicinity of the binding site of a complex (protein:ligand complex). I wanted to analyze the solvation differences in complexes formed by different proteins (with a very similar Carbohydrate Recognition Domain) bound to the same ligand.
>
> I was able to do it and visualize different properties but still, there are some things I don't fully understand.
>
> I represented the gist-g0.dx file. In order to assign some properties to the more localized water bridge regions, I treated the gist.out file. I ordered it based on the g_O values (higher to lower).
>
> Let's imagine the simplest case: two different proteins (binding site practically identical) bound to a ligand (the same one for both systems).
>
> I don't understand why for one of the systems the voxels with higher g_O values (top 30) are characterized by an Esw-norm ranging from -2.1 to -0.15 and an Eww-norm ranging from -0.15 to 0.46. However, for the other system, the voxels with higher g_O values are characterized by an Esw-norm ranging from -1.24 to 125576000 and an Eww-norm ranging from -0.10 to 0.30.
>
> Could someone help me understand these numbers? Conditions for the GIST calculations are exactly the same. What is more curious, there are cases in which two consecutive voxels (pointing to the same water region), display so different Esw-norm values (boxel 43573, Esw-norm=946.391 vs boxel 43572, Esw-norm = -0.128).
>
> I would really appreciate if some experts in the field could help me.
>
> Thanks in advance!
>
> Best,
>
> Reyes
>
>
>
>
>
> Reyes Núñez Franco (PhD student)
>
> CIC bioGUNE
>
>
>
>
>
> [ASOCIACION CENTRO DE INVESTIGACION COOPERATIVA EN BIOCIENCIAS]
>
> All information and attachments included in this email are confidential and intended for the original recipient only. You must not share any part of this message with any third party. If you have received this message by mistake, please reply to let us know immediately, so we can make sure such a mistake does not happen again and delete this message from your system. We have put our efforts into ensuring that this message is free of errors and viruses but, despite our efforts, you should always scan all emails for any threats with proper software, as the sender does not accept liability for any damage inflicted by viewing the content of this email. We also would like to inform that all personal data included in this email or in the attached documentation, could be processed by ASOCIACION CENTRO DE INVESTIGACION COOPERATIVA EN BIOCIENCIAS with the aim of managing the legal relationship between the parties. If you would like to access, modify or oppose to the treatment of your data, you can contact us in writing.
>
> Please consider the environment before printing this email
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
[ASOCIACION CENTRO DE INVESTIGACION COOPERATIVA EN BIOCIENCIAS]
All information and attachments included in this email are confidential and intended for the original recipient only. You must not share any part of this message with any third party. If you have received this message by mistake, please reply to let us know immediately, so we can make sure such a mistake does not happen again and delete this message from your system. We have put our efforts into ensuring that this message is free of errors and viruses but, despite our efforts, you should always scan all emails for any threats with proper software, as the sender does not accept liability for any damage inflicted by viewing the content of this email. We also would like to inform that all personal data included in this email or in the attached documentation, could be processed by ASOCIACION CENTRO DE INVESTIGACION COOPERATIVA EN BIOCIENCIAS with the aim of managing the legal relationship between the parties. If you would like to access, modify or oppose to the treatment of your data, you can contact us in writing.
Please consider the environment before printing this email
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Sep 13 2021 - 03:30:02 PDT
