AMBER: cis - trans isomerization

From: Nina Fischer <nfischer.informatik.uni-tuebingen.de>
Date: Thu, 27 Nov 2008 15:39:21 +0100

Dear all,

I am trying to perform a cis->trans isomerization of PRO during a MD
simulation using AMBER 10.
Unfortunately I don't end up with the trans form of PRO in the end. I
also tried amongst others a stepwise procedure (e.g. incrementally
increasing the angle by 2° using 88 segments of 5 ps) which was also not
working. Could anyone give me a hint?

Thanks a lot in advance!

Sincerely,
  Nina


My input file (similar to the example file in the AMBER 10 manual p.109)

nstlim=100000, cut=99.0, igb=1, saltcon=0.1,
ntpr=100, ntwr=100000, ntt=3, gamma_ln=5.0,
ntx=1, irest=0, ig = 256251,
ntc=2, ntf=2, tol=0.000001,
dt=0.002, ntb=0, tempi=300., temp0=300.,
jar=1,
/
&wt type='DUMPFREQ', istep1=1, /
&wt type='END', /
DISANG=tors.RST
DUMPAVE=tors_vs_t.dump
LISTIN=POUT
LISTOUT=POUT


My tors.RST file
(produced by makeANG_RST with the following constraint: 107 PRO OMEGA
170.0 190.0):
  &rst iat = 1620, 1635, 1637, 1647,
           r1 = 169.0, r2 = 170.0, r3 = 190.0, r4 = 191.0,
           rk2 = 2.0, rk3 = 2.0, &end


My tors_vs_t.dump (output):

    170.00000 1.38157 0.20546 0.00000
    170.00000 1.38157 0.20546 -0.00035
    169.99830 0.47077 0.20656 -0.00070
    169.99660 0.08517 0.20703 -0.00105
    169.99490 0.34210 0.20672 -0.00140
    169.99320 1.25996 0.20560 -0.00175
    169.99150 2.12485 0.20454 -0.00210
    169.98980 3.49226 0.20287 -0.00245
    169.98810 4.91085 0.20114 -0.00279
    169.98640 6.09679 0.19969 -0.00313
    169.98470 7.01486 0.19857 -0.00347
    169.98300 7.60174 0.19786 -0.00381
    169.98130 7.87793 0.19752 -0.00414
    169.97960 7.62205 0.19783 -0.00448
    169.97790 6.49825 0.19919 -0.00482
    169.97620 4.31581 0.20185 -0.00516
    169.97450 1.44021 0.20535 -0.00550
    169.97280 358.57785 -0.22981 -0.00548
    169.97110 355.02320 -0.22548 -0.00510
    169.96940 350.97144 -0.22055 -0.00472
    169.96770 347.28227 -0.21605 -0.00435
    169.96600 344.66392 -0.21286 -0.00398
    169.96430 342.96674 -0.21080 -0.00362
    169.96260 342.63136 -0.21039 -0.00326
    169.96090 343.02748 -0.21088 -0.00291
    169.95920 344.48603 -0.21266 -0.00255
    169.95750 346.34113 -0.21492 -0.00218
    169.95580 348.20425 -0.21719 -0.00181
        .
        .
        .
      0.01700 13.30929 -0.01620 -13.77487
      0.01530 12.20643 -0.01485 -13.77484
      0.01360 11.51945 -0.01402 -13.77482
      0.01190 11.69342 -0.01423 -13.77479
      0.01020 12.56779 -0.01530 -13.77477
      0.00850 13.41669 -0.01634 -13.77474
      0.00680 14.71493 -0.01792 -13.77471
      0.00510 16.05183 -0.01955 -13.77468
      0.00340 16.97904 -0.02068 -13.77465
      0.00170 17.85735 -0.02176 -13.77461



-- 
Eberhard Karls University Tübingen
Wilhelm Schickard Institut for Computer Science
Division for Simulation of Biological Systems
Room C314, Sand 14, D-72076 Tübingen
phone:	+49-7071-29-70456
fax:	  +49-7071-29-5152
email:	nfischer.informatik.uni-tuebingen.de
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo.scripps.edu
Received on Fri Dec 05 2008 - 17:59:40 PST
Custom Search