Please see the part of leap log about the saved files
>savepdb x lov_wation.pdb
Writing pdb file: lov_wation.pdb
printing CRYST1 record to PDB file with box info
Warning: Converting N-terminal residue name to PDB format: NPHE -> PHE
Warning: Converting C-terminal residue name to PDB format: CLEU -> LEU
> x = loadpdb lov_wation.pdb
Loading PDB file: ./lov_wation.pdb
total atoms in file: 22919
> saveamberparm x lov_wation.prmtop lov_wation.rst7
Checking Unit.
Warning: The unperturbed charge of the unit (0.155103) is not integral.
Warning: The unperturbed charge of the unit (0.155103) is not zero.
Note: Ignoring the error and warnings from Unit Checking.
Building topology.
Building atom parameters.
Building bond parameters.
Building angle parameters.
Building proper torsion parameters.
Building improper torsion parameters.
total 480 improper torsions applied
Building H-Bond parameters.
Incorporating Non-Bonded adjustments.
Not Marking per-residue atom chain types.
Marking per-residue atom chain types.
(Residues lacking connect0/connect1 -
these don't have chain types marked:
res total affected
CLEU 1
FMC 1
NPHE 1
WAT 5130
)
(no restraints)
Amit
On Fri, Dec 11, 2020 at 9:40 PM Carlos Simmerling <
carlos.simmerling.gmail.com> wrote:
> I'm still not seeing the part where you save your files. To be any help I
> need to see the leap inputs. Perhaps others can help out with the info that
> you've provided though.
>
> On Fri, Dec 11, 2020 at 11:01 AM Amit Sharma (Asstt. Prof., MCARS) <
> asharma4.jmi.ac.in> wrote:
>
> > last four lines of rst7 are as follows:
> >
> > "4.4530000 -1.4530000 -33.5420000 7.4510000 -4.4900000 -37.0780000
> > 6.6460000 -4.6000000 -36.7610000 7.3220000 -4.3890000 -37.9350000
> > 7.3130000 -4.4920000 -37.1570000 2.0370000 -1.4170000 -35.4280000
> > 2.5810000 -1.7580000 -34.8380000 2.0900000 -0.5590000 -35.2860000
> > 2.1260000 -1.3420000 -35.3200000"
> >
> >
> > and the starting input lines from leap are as follows:
> >
> > "Welcome to LEaP!
> > (no leaprc in search path)
> > > source leaprc.protein.ff14SBonlysc
> > ----- Source:
> >
> >
> /home/amit/Amber20_complete/tools/amber20/dat/leap/cmd/leaprc.protein.ff14SBonlysc
> > ----- Source of
> >
> >
> /home/amit/Amber20_complete/tools/amber20/dat/leap/cmd/leaprc.protein.ff14SBonlysc
> > done
> > Log file: ./leap.log
> > Loading parameters:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/parm/parm10.dat
> > Reading title:
> > PARM99 + frcmod.ff99SB + frcmod.parmbsc0 + OL3 for RNA
> > Loading parameters:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/parm/frcmod.ff14SB
> > Reading force field modification type file (frcmod)
> > Reading title:
> > ff14SB protein backbone and sidechain parameters
> > Loading parameters:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/parm/frcmod.ff99SB14
> > Reading force field modification type file (frcmod)
> > Reading title:
> > ff99SB backbone parameters (Hornak & Simmerling) with ff14SB atom types
> > Loading library:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/lib/amino12.lib
> > Loading library:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/lib/aminoct12.lib
> > Loading library:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/lib/aminont12.lib
> > > source leaprc.gaff2
> > ----- Source:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/cmd/leaprc.gaff2
> > ----- Source of
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/cmd/leaprc.gaff2 done
> > Log file: ./leap.log
> > Loading parameters:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/parm/gaff2.dat
> > Reading title:
> > AMBER General Force Field for organic molecules (Version 2.11, May 2016)
> > > source leaprc.water.opc
> > ----- Source:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/cmd/leaprc.water.opc
> > ----- Source of
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/cmd/leaprc.water.opc
> > done
> > Loading library:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/lib/atomic_ions.lib
> > Loading library:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/lib/solvents.lib
> > Loading parameters:
> > /home/amit/Amber20_complete/tools/amber20/dat/leap/parm/frcmod.opc
> > Reading force field modification type file (frcmod)
> > Reading title:
> > Parameters for OPC water (JPCL, 2014, 5 (21), pp 3863–3871)
> > Loading parameters:
> >
> >
> /home/amit/Amber20_complete/tools/amber20/dat/leap/parm/frcmod.ionsjc_tip4pew
> > Reading force field modification type file (frcmod)
> > Reading title:
> > Monovalent ion parameters for Ewald and TIP4P/EW water from Joung &
> > Cheatham JPCB (2008)
> > Loading parameters:
> >
> >
> /home/amit/Amber20_complete/tools/amber20/dat/leap/parm/frcmod.ions234lm_126_tip4pew
> > Reading force field modification type file (frcmod)
> > Reading title:
> > Li/Merz ion parameters of divalent to tetravalent ions for TIP4P/EW water
> > model (12-6 normal usage set)
> > > loadAmberParams FMC.frcmod
> > Loading parameters: ./FMC.frcmod
> > Reading force field modification type file (frcmod)
> > Reading title:
> > this was copied from FMCcys but then edited to make the protein part (CB,
> > SG etc) match ff14SB CYX atom types (c3, ss to 2C, S) since that's what
> > will bond to this residue. note that c3 showed up more than just the CB
> > atom
> > > loadoff FMC.lib
> > Loading library: ./FMC.lib"
> >
> > As you see there is no box info, I added the following box info at the
> last
> > line from the pdb
> > 67.7859143 67.7859143 67.7859143 109.4712190 109.4712190 109.471219
> >
> > and now i have set a minimization run through sander with following input
> >
> > &cntrl
> > imin = 1,
> > maxcyc = 100000,
> > ncyc = 50000,
> > ntb = 0, ! no periodicity
> > cut = 99, ! non bond cut off
> > ntr = 1,
> > restraintmask= ':1 - 145 & !(@H=) &!(:CYX|:FMC=)',
> > restraint_wt=100.0,
> > /
> >
> >
> >
> > Amit
> >
> >
> > On Fri, Dec 11, 2020 at 9:01 PM Carlos Simmerling <
> > carlos.simmerling.gmail.com> wrote:
> >
> > > I don't see anything unusual, maybe others do... but it's hard to know
> > > without seeing your actual leap input.
> > > perhaps you could also show the last few lines of the rst7 produced by
> > > leap.
> > >
> > > On Fri, Dec 11, 2020 at 9:53 AM Amit Sharma (Asstt. Prof., MCARS) <
> > > asharma4.jmi.ac.in> wrote:
> > >
> > > > 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
> > > >
> > > _______________________________________________
> > > 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
> >
> _______________________________________________
> 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 - 08:30:03 PST
