Re: [AMBER] Analysis of minimization stage

From: Valentina Romano <valentina.romano.unibas.ch>
Date: Mon, 12 May 2014 16:26:44 +0000

Analyzing the minimization output (the frames i obtained after the minimization), I realized that the 2 atoms involved in the previous error message belonged to the ligands and were not restrained during the minimization. Thus I changed this part in the minimization.in file:

Initial minimisation PknG-Adenine complex: minimization solvent + ions
 &cntrl
  imin=1,
  maxcyc=1000,
  ntb=1,
  ntr=1,
  cut=10,
  restraint_wt = 10.0,
  restraintmask=':1-246.CA,C,O,N | .3834-3848',
  ntpr=1, ntwx=1, ntwr=100

Atoms from 3834 to 3848 are ligand's atoms (with H atoms).

I analyzed the output frames and the ligand was much more stable then before.
(The energy was flat and no overlaps were detected).

This structure was used as strating point for a short MD (to relax solvent molecules):

15ps MD PknG-Adenine complex: restraints on PknG-Ade residues while the solvent is leeting free
 to relax
 &cntrl
  imin=0,
  irest=0,
  ntx=1,
  ig=-1,
  ntb=1,
  ntr=1,
  cut=10,
  ntc=2,
  ntf=2,
  tempi=300.0,
  temp0=300.0,
  ntt=3,
  gamma_ln=5.0,
  nstlim=15000, dt=0.001,
  ntpr=1, ntwx=1, ntwr=1000
  restraint_wt = 10.0,
  restraintmask=':1-246.CA,C,O,N | .3834-3848'
 /
Sander stopped again at the begining:

vlimit exceeded for step 1; vmax = 82.0601

     Coordinate resetting (SHAKE) cannot be accomplished,
     deviation is too large
     NITER, NIT, LL, I and J are : 0 3 1923 3841 3848

     Note: This is usually a symptom of some deeper
     problem with the energetics of the system.

If I restrained all ligand's atoms, why atoms 3841 and 3848 are still causing a problem?

Analizing the output file i saw a message related to the charge of the whole system:

     Sum of charges from parm topology file = 0.00099988
     Forcing neutrality...

Could it be a cause of my problem?

I ran the MD as:

sander -O -i PknGAde-wt-md.in -o PknGAde-wt-md.out -p ../PknGAde_params/PknGHAdeH_ion_wt.prmtop -c PknGAde-wt-min.rst -r PknGAde-wt-md.rst -ref PknGAde-wt-min.rst -x PknGAde-wt-md.mdcrd &

Thus the set of initial coord and the ref coord are the same (both are the minimized structure), is it correct?

Vale

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Valentina Romano | PhD Student | Biozentrum, University of Basel & SIB Swiss Institute of Bioinformatics
Klingelbergstrasse 61 | CH-4056 Basel |

Phone: +41 61 267 15 80


________________________________________
From: Jason Swails [jason.swails.gmail.com]
Sent: Monday, May 12, 2014 4:18 PM
To: amber.ambermd.org
Subject: Re: [AMBER] Analysis of minimization stage

On Mon, 2014-05-12 at 09:59 -0400, David A Case wrote:
> On Mon, May 12, 2014, Valentina Romano wrote:
> >
> > CPPTRAJ: Trajectory Analysis. V13.22
> > ___ ___ ___ ___
> > | \/ | \/ | \/ |
> > _|_/\_|_/\_|_/\_|_
> > AmberParm Title: [default_name]
> > Radius Set: modified Bondi radii (mbondi)
> > INPUT: Reading Input from file PknGAde-wt-min.check_overlap
> > [trajin PknGAde-wt-min.rst]
> > [PknGAde-wt-min.rst] contains 1 frames.
> > [checkoverlap]
>
> Try adding the "reportfile overlaps.dat" to the checkoverlap command.
> I'm always uncertain about when cpptraj writes things to stdout, and when
> they go to one of the data files.

Whenever I've used checkoverlaps before the report gets printed to
stdout...

--
Jason M. Swails
BioMaPS,
Rutgers University
Postdoctoral Researcher
_______________________________________________
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
Received on Mon May 12 2014 - 09:30:02 PDT
Custom Search