[gmx-users] Umbrella Sampling tutorial

shahab shariati shahab.shariati at gmail.com
Sat Nov 9 15:38:51 CET 2013


Dear Justin

Thanks for your explanation.

My system contains lipid bilayer + drug + water molecules.

I want to calculate Potential of mean force as a function of the distance
between the centers of mass of drug and the lipid bilayer.

Box vector along which the pulling is being conducted is Z.

1) Are these issues true about my system?

lipid bilayer = reference group for pulling.
drug molecule = group to which pulling force is applied.

Then, should I use position restraining on the lipid bilayer?

2) The system you considered was the dissociation of a single peptide from
the growing end of an protofibril.

So, in your summary_distances.dat, distance between chain A and chain B was
increased.

But, I want to consider translocation of the drug molecule from water into
the lipid bilayer. On the other hand, I want to consider approaching drug
molecule to lipid bilayer.

So, in my summary_distances.dat, distance between drug molecule and lipid
bilayer will be decrease.

My mean is that my case is contrary to your case. Nonetheless, should I use
exactly Pull code section of your md_pull.mdp file?

Best wishes



More information about the gromacs.org_gmx-users mailing list