Re: AMBER: Values of Spring Constants

From: David A. Case <case.scripps.edu>
Date: Wed, 1 Feb 2006 08:38:13 -0800

On Wed, Feb 01, 2006, Sai Kumar Ramadugu wrote:

> Actually I was trying to understand the functioning of Umbrella Smapling.
> For the same I took a protein of 40 aminoacids and did a GB simulation with
> restraints on phi angle of 2nd amino acid which is Alanine.
> My restraints input file is as follows:
> # 2 ALA PHI: (1 ASP C)-(2 ALA N)-(2 ALA CA)-(2 ALA C) 120.0 180.0
> &rst iat = 13, 15, 17, 23,
> r1 = 119.0, r2 = 120.0, r3 = 180.0, r4 = 181.0,
> rk2 = 2.0, rk3 = 2.0, &end

To use WHAM, you need to define a harmonic potential, where r2=r3 and r1 and
r4 are far away from r2 and r3, as described on page 138 of the Users' Manual.
The force constant will be twice the values you give for rk2=rk3 (because
Amber historically does not include a factor of 1/2 in these constraint
energies.)

...good luck...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Feb 02 2006 - 06:10:06 PST
Custom Search