Anthony Armstrong wrote:
Hello all,

I'd like to run REMD on a peptide system, but it turns out out my
system is a little to large, and I've determined that to get a swap
probablity of .2, a separation of temperature of only a couple degrees
between replicas would be required.  I ran across the paper of Liu et al.
(PNAS 2005 102:13749-13754) on Replica Exchange with Solute Tempering
which addresses this problem and which appears to attempt temperature
swaps based on a "modified" total potential in which solute/solute and
solute/protein interactions are first scaled.

Are you sure you want to do this? All but one of your replicas become non-physical and you can not do fun stuff like extracting kinetics info from the trajectory anymore (see David van der Spoel and M. Marvin Seibert: Protein Folding Kinetics and Thermodynamics from Atomistic Simulations, Phys. Rev. Lett. 96 pp. 238102, 2006).


My question is:

How might I go about calculating a user defined energy (in this case
E(X) = Ep(X) + A(T) * Ew(X) + B(T) * Epw -->  where X is the set of
system coordinates, subscripts p, w, and pw denote protein/protein,
water/water, and protein/water interactions, respectively, and A(T) and
B(T) are scale factors involving the target temperature for the given
replica) and then pass this energy rather than the unscaled potential
to gromacs for determining whether or not to swap temperatures.

Any help would be much appreciated.

The problem is in PME, for which you can not directly extract interaction energies. There is a solution for this, but it is not implemented. Otherwise everything is present, except for some additional parameters to determine your A and B functions.


Thanks,
Tony

Anthony A. Armstrong
Department of Biophysics and Biophysical Chemistry
Johns Hopkins University School of Medicine
725 North Wolfe Street
Baltimore, MD 21205
410-955-8715
410-955-0637 FAX
[EMAIL PROTECTED]

On Wed, 19 Jul 2006, Anthony Armstrong wrote:


Hello,

I'm hoping to get some tips on removing protein rotation inside a
triclinic water box.  Apologies if this has been covered recently, but
the most recent discussion I've been able to track down is a thread from
2003 which stated a patch was available to prevent such rotations and an
algorithm would be formally implemented into GROMACS "in the near future"
(Tsjerk Wassenaar t.a.wassenaar at chem.rug.nl Tue May 20 13:52:00 CEST
2003).  Has the implementation occurred?  If I use comm-grps = Protein and
comm-mode = angular?  Of course I'm in the process of trying it, but
thought some kind person from the Gromacs community might be able to save
me some time if it is going to fail.  Thanks in advance.

Tony

Anthony A. Armstrong
Department of Biophysics and Biophysical Chemistry
Johns Hopkins University School of Medicine
725 North Wolfe Street
Baltimore, MD 21205
410-955-8715
410-955-0637 FAX
[EMAIL PROTECTED]


_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


--
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,          75124 Uppsala, Sweden
phone:  46 18 471 4205          fax: 46 18 511 755
[EMAIL PROTECTED]       [EMAIL PROTECTED]   http://folding.bmc.uu.se
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to