From: Ashkan Shekaari (shekaari_at_email.kntu.ac.ir)
Date: Wed Jul 14 2021 - 05:50:01 CDT
One more thing,
Because the ligand is presumably much smaller than your protein of interest,
you may have to use QMMM.
Good luck with your calculations ...
---
,Best
Ashkan Shekaari
PhD Candidate in Solid State Physics
Department of Physics
K. N. Toosi University of Technology, Tehran, 15875-4416, Iran [ mailto:shekaari_at_email.kntu.ac.ir ] [ https://urldefense.com/v3/__https://scholar.google.com/citations?hl=en&user=LDzbXB4AAAAJ&view_op=list_works&sortby=pubdate__;!!DZ3fjg!t5MpHIp-EC7S9EzJNG0W5y1kmPlRQM9sKfdsCwEUeGLHd1D9Jyh6vuKFo7DwYt8nfg$ ]
From: "shekaari" <shekaari_at_email.kntu.ac.ir>
To: "namd-l" <namd-l_at_ks.uiuc.edu>, "M." <matabay2010_at_gmail.com>
Sent: Wednesday, July 14, 2021 3:12:02 PM
Subject: Re: namd-l: question about free energy calculation through FEP method in NAMD
Hi Maryam,
Give it a try using
1- SMD and then the Jarzynski's equality
or
2- replica-exchange umbrella sampling.
Both are implemented in NAMD.
---
,Best
Ashkan Shekaari
PhD Candidate in Solid State Physics
Department of Physics
K. N. Toosi University of Technology, Tehran, 15875-4416, Iran [ mailto:shekaari_at_email.kntu.ac.ir ] [ https://urldefense.com/v3/__https://scholar.google.com/citations?hl=en&user=LDzbXB4AAAAJ&view_op=list_works&sortby=pubdate__;!!DZ3fjg!sqBe6JihE0aFs8dIlpjvF3p5aILDvQuzCcsnkUAAgrEyl9Gst96rqy7nfCWBxxkPQA$ ]
From: "M. A" <matabay2010_at_gmail.com>
To: "namd-l" <namd-l_at_ks.uiuc.edu>, "Jérôme Hénin" <jerome.henin_at_ibpc.fr>
Sent: Wednesday, July 14, 2021 11:05:05 AM
Subject: namd-l: question about free energy calculation through FEP method in NAMD
Dear Jérôme,
Hi
Do you have any suggestions to solve this problem?
I want to calculate the free energy ( ΔG ) between a ligand and protein through FEP method in NAMD. In my structure, the ligand is located on the surface of the protein, not inside the protein, lt means a part of the ligand that is closed to the protein surface has interaction with protein, and the other part of the ligand which is in contact with water molecules, has interaction with these molecules. So, when I change the B column of atoms of ligand to -1 (in solvate.fep file), the energy (VdW and Elc) turns off for all atoms of ligand (all atoms in contact with protein and water molecules both) while I need to eliminate the only interaction between the ligand and protein. How can I specify that it eliminates just the interaction between the ligand and protein and does not consider the interaction of ligand and water molecules.
I appreciate all your help.
best regards
Maryam
-- Maryam Atabay Ph.D in Nano Computational Chemistry Department of Physics, Sharif University of Technology
This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:11 CST