Re: FEP with fixed explicit environment?

From: Sebastian Stolzenberg (
Date: Fri Mar 13 2009 - 10:58:32 CDT

Hi Chris,

thanks for your quick reply.

There is another problem I see:
Simply fixing the environment will give me RATTLE violations, so I guess
I first need to impose weak constraints that I gradually turn into a
I see no way to calculate the free energy difference for this
transformation, or is there any?


Chris Harrison wrote:
> Sebastian,
> Simple question first: Can NAMD do this?
> Yes. There's nothing in the code that prevents it. Can't think of a reason NAMD wouldn't execute successfully.
> Hard question: Should you do this?
> During any alchemical perturbation there is the possibility that the environment dynamically responds by rearranging its conformation. If you can justify that any changes of the conformational ensemble that occur during a restrained R->A perturbation are not significantly different from conformational ensemble changes that occur during an unrestrained R->A perturbation, then you may be able to do this. So, the atoms beyond your 30 Ang radius would have to fullfill the criteria that the dynamics of those atoms you wish to restrain do not respond to the R->A perturbation and the unrestrained atoms' dynamical response(s) to the perturbation is not altered by the presence of the restrained atoms.
> Please let me know how it goes. I'm interested to know if it works successfully {meaning a) no NAMD crashes & b) you get a correct/reasonable result}.
> C.
> --
> Chris Harrison, Ph.D.
> Theoretical and Computational Biophysics Group
> NIH Resource for Macromolecular Modeling and Bioinformatics
> Beckman Institute for Advanced Science and Technology
> University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801
> Voice: 217-244-1733
> Fax: 217-244-6078
> Sebastian Stolzenberg <> writes:
>> Date: Thu, 12 Mar 2009 19:50:32 -0400
>> From: Sebastian Stolzenberg <>
>> To:
>> Subject: namd-l: FEP with fixed explicit environment?
>> Return-Path:
>> Message-ID: <>
>> X-Spam-Status: No, score=-2.2 required=5.0 tests=AWL,BAYES_00
>> autolearn=unavailable version=3.1.7-0+tcb1
>> Dear Everybody,
>> I have an equilibrated NPT structure of a protein in explicit
>> lipid/solvent with periodic boundary conditions. Let's assume I do a
>> mutation R105A with dual-topology FEP. Certainly, I will also need to
>> transform a bulk water molecule (WAT) into a sodium (SOD) to keep the
>> net charge=0. To get the final free energy difference, I will subtract
>> delta_G(WAT->SOD) that I get from a separate run.
>> The system is large, I was thinking about fixing all atoms of the system
>> except for the ones around ~30A of the R105A mutation and the WAT->SOD
>> transformation. (Of course, I will not have a fixation boundary crossing
>> covalent bonds that leads e.g. to RATTLE constraint violations).
>> Is this feasible with NAMD-FEP? Any troubles with periodic boundary
>> conditions? I know that all could be locally done with implicit
>> lipid/solvent in CHARMM, which I would like to avoid for now.
>> Thanks so much,
>> Sebastian

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:52:29 CST