Dear Víctor:

what you should do depends on what you want to get. I suggest that you simply cut your PMF at the point where the bilayer is "lost" and explain why it is cut. If you really want to pursue your idea, then you need to have a good grasp of exactly what you are doing by such a method. You would want to start by calculating the free energy of applying this additional restraint to a bilayer that is not interacting with your solute. The problem then is that you would really need to also calculate the free energy of releasing this restraint upon reaching the center of the bilayer, and it sounds like that not only unfeasible, but also uninstructive. My last thought is that you'll want to ensure that the destruction of the bilayer is actually the equilibrium action of your solute. Perhaps you simply don't have the equilibrated orientation, etc, that would allow this system to be stable. In that case, you might consider using your method as a way to generate the starting structure, which would then be simulated in the absence of this additional restraint for data collection.

In the end, you need to have a clear idea of what you want the PMF for and also have a clear idea of your path fro the unrestrained state to your umbrellas. Free energy is a state property, but you need to understand the full path, and if you forget to include the cost of something like adding the additional restraint to a neat bilayer, then your study will be of little value. Note that the relative free energies in the state with the additional restraints are not useful in and of themselves just as the intermediates of any alchemical transformation are not useful on their own and are only useful as a means to get to the endpoint.

Chris.


-- original message --

I am using umbrella sampling method for calculate the PMF for one
molecule in a bilayer system made of SDS and decanol. when the molecule
is in the middle of the bilayer, generate a hole and the bilayer system
is lost.

I am thinking use position restraint for one of the atoms in the SDS (sodium dodecyl sulfate). is this a good idea?

any suggestion will be well received



Thanks in advance

Víctor E. Bahamonde Padilla
Laboratorio de Fisicoquímica Molecular
Departamento de Química
Facultad de Ciencias
Universidad de Chile
Phone: 562-978-7443
vedua...@ug.uchile.cl



--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to