Re: [AMBER] NMR restraints

From: Debarati DasGupta <debarati_dasgupta.hotmail.com>
Date: Sat, 31 Oct 2020 14:46:15 +0000

Hello Prof. Case,
I just ran a 100 step minimization to see how the NMR restraints are applied.
So in the min.out file; after the part where the TI masks are read, I found this section on restraints.
Is this getting printed out on the out file because I have set LISTIN=POUT?

I see the 3 distance restraints as mentioned in the file output.disang; but I am still not sure if they are being read properly, why is the NMR restraints bond,angle,torsion all 0 in all the steps of minimization?
Also this is my output.disang file contents:
[cid:image002.png.01D6AF73.0CEDA8D0]

This is the section of my minimization output:


Ewald parameters:
     verbose = 0, ew_type = 0, nbflag = 1, use_pme = 1
     vdwmeth = 1, eedmeth = 1, netfrc = 0
     Box X = 84.482 Box Y = 81.908 Box Z = 87.672
     Alpha = 90.000 Beta = 90.000 Gamma = 90.000
     NFFT1 = 90 NFFT2 = 90 NFFT3 = 90
     Cutoff= 8.000 Tol =0.100E-04
     Ewald Coefficient = 0.34864
     Interpolation order = 4
     TI Mask 1 :IPA; matches 12 atoms
     TI Mask 2 matches 0 atoms
     TI region 1: 50200 atoms
     TI region 2: 50188 atoms
     SC Mask 1 :IPA; matches 12 atoms

--------------------------------------------------------------------------------
   3. ATOMIC COORDINATES AND VELOCITIES
--------------------------------------------------------------------------------

default_name
begin time read from input coords = 0.000 ps



           Begin reading energy term weight changes/NMR restraints
WEIGHT CHANGES:
DUMPFREQ 1 0 0.000000 0.000000 0 0
                         ** No weight changes given **

RESTRAINTS:
Requested file redirections:
  DISANG = output.disang
  DUMPAVE = dist_0.00922.dat
  LISTIN = POUT
  LISTOUT = POUT
Restraints will be read from file: output.disang
Here are comments from the DISANG input file:

******
C2 ( 4396)-CA ( 204) NSTEP1= 0 NSTEP2= 0
R1 = 5.608 R2 = 6.108 R3 = 6.108 R4 = 6.608 RK2 = 30.000 RK3 = 30.000
Rcurr: 6.313 Rcurr-(R2+R3)/2: 0.205 MIN(Rcurr-R2,Rcurr-R3): 0.205
******
C2 ( 4396)-CA ( 2132) NSTEP1= 0 NSTEP2= 0
R1 = 6.333 R2 = 6.833 R3 = 6.833 R4 = 7.333 RK2 = 30.000 RK3 = 30.000
Rcurr: 6.883 Rcurr-(R2+R3)/2: 0.050 MIN(Rcurr-R2,Rcurr-R3): 0.050
******
C2 ( 4396)-CA ( 1387) NSTEP1= 0 NSTEP2= 0
R1 = 5.166 R2 = 5.666 R3 = 5.666 R4 = 6.166 RK2 = 30.000 RK3 = 30.000
Rcurr: 5.641 Rcurr-(R2+R3)/2: 0.025 MIN(Rcurr-R2,Rcurr-R3): 0.025
                       Number of restraints read = 3

                  Done reading weight changes/NMR restraints



Number of triangulated 3-point waters found: 15264
Number of shake restraints removed in TI region 1 : 0
Number of shake restraints removed in TI region 2 : 0




*********************************************************************************************************************************************************************

From: David A Case<mailto:david.case.rutgers.edu>
Sent: 30 October 2020 21:22
To: AMBER Mailing List<mailto:amber.ambermd.org>
Subject: Re: [AMBER] NMR restraints

On Fri, Oct 30, 2020, Debarati DasGupta wrote:

>I was trying to set up a minimization (TI run) and I see that NMR bond
>angle Torsion energies are all 0.

>LISTIN=POUT

To add to what Carlos said: if you set LISTIN=POUT, you should see a
detailed printout of all the NMR calculations on the initial step. You
can check to see if they are defined the way you want, and see how the
energies are computed.

[If for some reason this info doesn't show up in parallel, just run a
single step calculation in serial mode.]

....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber



_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber

FE11AF737B1C4870B092E7D9DCFB9142.png
(image/png attachment: FE11AF737B1C4870B092E7D9DCFB9142.png)

Received on Sat Oct 31 2020 - 08:00:02 PDT
Custom Search