Re: [AMBER] Box info not found in prmtop file

From: Amit Sharma (Asstt. Prof., MCARS) <"Amit>
Date: Fri, 11 Dec 2020 20:22:32 +0530

My tleap input is as follows

> source leaprc.protein.ff14SBonlysc
> source leaprc.gaff2
Joung-Chetham monovalent ions for> source leaprc.water.opc
This gives (nothing unusual)



> loadAmberParams FMC.frcmod
> loadoff FMC.lib
> x = loadpdb protein.pdb (where protein.pdb is protein containing FMC
(FMN))
> solvateoct x OPCBOX 8.0
This gives
"Scaling up box by a factor of 1.309835 to meet diagonal cut criterion
  Solute vdw bounding box: 52.851 41.765 42.392
  Total bounding box for atom centers: 73.808 73.808 73.808
      (box expansion for 'iso' is 37.1%)
  Solvent unit box: 18.865 18.478 19.006
The number of boxes: x= 4 y= 4 z= 4
  Volume: 209323.503 A^3 (oct)
  Total mass 109655.332 amu, Density 0.870 g/cc
  Added 4991 residues."


and then on checking the molecule x gives non integral, non zero charge of
-5.844897. This becomes +0.155103 after I add 6 Na+ ions.

then I save rst7 and prmtop files. And, at this stage I do not get box info
in rst7.

Amit







On Fri, Dec 11, 2020 at 7:01 PM Carlos Simmerling <
carlos.simmerling.gmail.com> wrote:

> I suggest looking in your leap.log file to see what leap responded when you
> asked to solvate the system. if you still have trouble, sharing that part
> of the log file along with your leap input may make it easier for people to
> help you debug.
>
> On Fri, Dec 11, 2020 at 8:14 AM Amit Sharma (Asstt. Prof., MCARS) <
> asharma4.jmi.ac.in> wrote:
>
> > Hi,
> >
> > I solvated my protein+chromophore system in an OPC water model using 8
> > Angstrom OPCBOX.
> >
> > But, last line of rst7 file does not contain box info. Hence, initial
> > minimization (with ntb = 1) is not working. Although, the issue does not
> > arise when using a tip3p water box (but, i need to use opc water).
> >
> > Then, I tried minimization after copying the box info from the pdb, but
> > this aswell stops with the following out info
> >
> > "NFFT1 = 72 NFFT2 = 375 NFFT3 = 375
> > Cutoff= 8.000 Tol =0.100E-04
> > Ewald Coefficient = 0.34864
> > Interpolation order = 4
> >
> > Ewald PARAMETER RANGE CHECKING:
> > parameter alpha: (unit cell angle) has value 0.31095E+04
> > This is outside the legal range
> > Lower limit: 0.30000E+02 Upper limit: 0.15000E+03
> > Check ew_legal.h"
> >
> > Any suggestions on how to resolve this issue?
> >
> > Thanks,
> >
> > Amit
> >
> > --
> > _
> > _
> >
> >
> >
> > *Please consider the environment before printing this email.*
> > The
> > information contained in this electronic message and any *attachments to
> > this message are intended for the exclusive use of addressee(s) and may
> > contain proprietary, confidential or privileged information. If you are
> > not
> > the intended recipient, you should not disseminate, distribute **or copy
> > this e-mail. Please notify the
> > sender immediately and destroy all copies of
> > this message and any attachments.***
> > **
> > **WARNING***: The recipient should
> > check this email and any attachments for
> > the presence of viruses. JMI
> > accepts no liability for any damage caused by any
> > virus transmitted by this
> > email.*
> > _______________________________________________
> > 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
>

-- 
_
_
*Please consider the environment before printing this email.*
The 
information contained in this electronic message and any *attachments to 
this message are intended for the exclusive use of addressee(s) and may 
contain proprietary, confidential or privileged information. If you are not 
the intended recipient, you should not disseminate, distribute **or copy 
this e-mail. Please notify the
sender immediately and destroy all copies of 
this message and any attachments.***
**
**WARNING***: The recipient should 
check this email and any attachments for
the presence of viruses. JMI 
accepts no liability for any damage caused by any
virus transmitted by this 
email.*
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Dec 11 2020 - 07:00:03 PST
Custom Search