[gmx-users] Re: Error during executing the protein ligand tutorial
Mark Abraham
Mark.Abraham at anu.edu.au
Tue Jun 21 08:44:06 CEST 2011
On 21/06/2011 4:33 PM, bharat gupta wrote:
> So, after adding 1 NA ion, I started with energy mimization, but I am
> getting the following error after md run command :-
Please search for help first.
http://www.gromacs.org/Documentation/Errors#LINCS.2fSETTLE.2fSHAKE_warnings
Mark
> Step= 1, Dmax= 1.0e-02 nm, Epot= 1.36889e+09 Fmax= 1.97591e+11,
> atom= 1248
> Step= 2, Dmax= 1.2e-02 nm, Epot= 5.46814e+07 Fmax= 3.57239e+09,
> atom= 972
> Step= 3, Dmax= 1.4e-02 nm, Epot= 8.68244e+06 Fmax= 5.80655e+08,
> atom= 1253
> Step= 4, Dmax= 1.7e-02 nm, Epot= 3.33619e+06 Fmax= 9.98000e+07,
> atom= 1248
> Step= 5, Dmax= 2.1e-02 nm, Epot= 2.30163e+06 Fmax= 1.19057e+07,
> atom= 1253
> Step= 6, Dmax= 2.5e-02 nm, Epot= 2.07254e+06 Fmax= 3.20568e+06,
> atom= 1304
> Step= 7, Dmax= 3.0e-02 nm, Epot= 1.97667e+06 Fmax= 1.59513e+06,
> atom= 1045
> Step= 8, Dmax= 3.6e-02 nm, Epot= 1.85877e+06 Fmax= 8.25518e+05,
> atom= 3616
> Step= 9, Dmax= 4.3e-02 nm, Epot= 1.66148e+06 Fmax= 7.69483e+05,
> atom= 3616
> Step= 10, Dmax= 5.2e-02 nm, Epot= 1.43976e+06 Fmax= 6.47087e+05,
> atom= 3616
> Step= 11, Dmax= 6.2e-02 nm, Epot= 1.16869e+06 Fmax= 6.39015e+05,
> atom= 3616
> Step= 12, Dmax= 7.4e-02 nm, Epot= 8.84664e+05 Fmax= 5.75494e+05,
> atom= 3616
> Step= 13, Dmax= 8.9e-02 nm, Epot= 6.62825e+05 Fmax= 1.01984e+07,
> atom= 47956
>
> step 14: Water molecule starting at atom 47956 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 15, Dmax= 5.3e-02 nm, Epot= 5.67330e+05 Fmax= 5.36987e+05,
> atom= 36166
> Step= 16, Dmax= 6.4e-02 nm, Epot= 3.51341e+05 Fmax= 4.19647e+05,
> atom= 3619
> Step= 17, Dmax= 7.7e-02 nm, Epot= 6.34731e+04 Fmax= 1.56555e+06,
> atom= 3616
> Step= 18, Dmax= 9.2e-02 nm, Epot= -2.03364e+04 Fmax= 4.60916e+05,
> atom= 3616
> Step= 19, Dmax= 1.1e-01 nm, Epot= -3.50504e+05 Fmax= 1.53090e+06,
> atom= 213514
>
> step 20: Water molecule starting at atom 213514 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 20, Dmax= 1.3e-01 nm, Epot= -4.11692e+05 Fmax= 4.49357e+05,
> atom= 3619
> Step= 21, Dmax= 1.6e-01 nm, Epot= -4.87207e+05 Fmax= 3.39700e+07,
> atom= 3629
> Step= 22, Dmax= 1.9e-01 nm, Epot= -7.85503e+05 Fmax= 3.56083e+05,
> atom= 3619
>
> step 23: Water molecule starting at atom 206593 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
>
> step 23: Water molecule starting at atom 212776 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Wrote pdb files with previous and current coordinates
> Step= 23, Dmax= 2.3e-01 nm, Epot= 1.95863e+09 Fmax= 6.76309e+11,
> atom= 3633
> step 24: Water molecule starting at atom 212776 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 25, Dmax= 5.8e-02 nm, Epot= -9.28980e+05 Fmax= 2.81692e+05,
> atom= 362976
>
> step 26: Water molecule starting at atom 42463 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 27, Dmax= 3.5e-02 nm, Epot= -1.03111e+06 Fmax= 2.76677e+05,
> atom= 36293
> Step= 28, Dmax= 4.1e-02 nm, Epot= -1.14733e+06 Fmax= 2.66398e+05,
> atom= 3629
> Step= 29, Dmax= 5.0e-02 nm, Epot= -1.26821e+06 Fmax= 2.51779e+05,
> atom= 3629
> Step= 30, Dmax= 6.0e-02 nm, Epot= -1.42873e+06 Fmax= 2.31393e+05,
> atom= 3629
> Step= 31, Dmax= 7.2e-02 nm, Epot= -1.59623e+06 Fmax= 3.53379e+06,
> atom= 3633
>
> step 32: Water molecule starting at atom 292885 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 33, Dmax= 4.3e-02 nm, Epot= -1.62300e+06 Fmax= 2.45780e+05,
> atom= 361985
> Step= 34, Dmax= 5.2e-02 nm, Epot= -1.75124e+06 Fmax= 2.87155e+05,
> atom= 213958
> Step= 35, Dmax= 6.2e-02 nm, Epot= -1.86014e+06 Fmax= 6.76483e+05,
> atom= 3634
>
> step 36: Water molecule starting at atom 256024 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 36, Dmax= 7.4e-02 nm, Epot= -1.91639e+06 Fmax= 2.49523e+05,
> atom= 3619
> Step= 37, Dmax= 8.9e-02 nm, Epot= -2.11154e+06 Fmax= 4.87466e+05,
> atom= 3636
> Step= 38, Dmax= 1.1e-01 nm, Epot= -2.20858e+06 Fmax= 1.92308e+05,
> atom= 3629
> Step= 40, Dmax= 6.4e-02 nm, Epot= -2.34036e+06 Fmax= 3.07614e+06,
> atom= 3634
> Step= 41, Dmax= 7.7e-02 nm, Epot= -2.36181e+06 Fmax= 1.95100e+05,
> atom= 3634
> Step= 42, Dmax= 9.2e-02 nm, Epot= -2.55367e+06 Fmax= 1.66837e+05,
> atom= 3629
> Step= 43, Dmax= 1.1e-01 nm, Epot= -2.67052e+06 Fmax= 1.01067e+07,
> atom= 3618
> Step= 44, Dmax= 1.3e-01 nm, Epot= -2.76975e+06 Fmax= 2.25183e+05,
> atom= 3618
> Step= 45, Dmax= 1.6e-01 nm, Epot= -2.95070e+06 Fmax= 5.66364e+05,
> atom= 3614
>
> step 46: Water molecule starting at atom 249877 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
>
> step 46: Water molecule starting at atom 85966 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Wrote pdb files with previous and current coordinates
> Step= 46, Dmax= 1.9e-01 nm, Epot= -2.97755e+06 Fmax= 9.94270e+06,
> atom= 249877
> Step= 47, Dmax= 2.3e-01 nm, Epot= -2.99471e+06 Fmax= 2.40481e+05,
> atom= 3614
> Step= 48, Dmax= 2.8e-01 nm, Epot= -3.24066e+06 Fmax= 2.63430e+06,
> atom= 3608
>
> step 49: Water molecule starting at atom 82687 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
>
> step 49: Water molecule starting at atom 129469 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Wrote pdb files with previous and current coordinates
> Step= 49, Dmax= 3.3e-01 nm, Epot= 1.19344e+08 Fmax= 1.97429e+10,
> atom= 3608
> step 50: Water molecule starting at atom 129469 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 51, Dmax= 8.3e-02 nm, Epot= -3.28700e+06 Fmax= 3.42115e+05,
> atom= 3620
> Step= 52, Dmax= 9.9e-02 nm, Epot= -3.35221e+06 Fmax= 5.71868e+05,
> atom= 3620
>
> step 53: Water molecule starting at atom 82687 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 53, Dmax= 1.2e-01 nm, Epot= -3.37935e+06 Fmax= 2.23519e+05,
> atom= 3620
> Step= 54, Dmax= 1.4e-01 nm, Epot= -3.48138e+06 Fmax= 1.24343e+06,
> atom= 3620
>
> step 55: Water molecule starting at atom 82687 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 55, Dmax= 1.7e-01 nm, Epot= -3.47365e+06 Fmax= 4.00718e+07,
> atom= 82687
> step 56: Water molecule starting at atom 82687 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 56, Dmax= 8.6e-02 nm, Epot= -3.49998e+06 Fmax= 1.67258e+05,
> atom= 3620
> Step= 57, Dmax= 1.0e-01 nm, Epot= -3.59531e+06 Fmax= 9.54825e+04,
> atom= 3634
> Step= 59, Dmax= 6.2e-02 nm, Epot= -3.67520e+06 Fmax= 4.07963e+05,
> atom= 3622
> Step= 60, Dmax= 7.4e-02 nm, Epot= -3.69432e+06 Fmax= 1.64972e+05,
> atom= 3622
> Step= 61, Dmax= 8.9e-02 nm, Epot= -3.75825e+06 Fmax= 2.68322e+05,
> atom= 3613
> Step= 62, Dmax= 1.1e-01 nm, Epot= -3.78424e+06 Fmax= 2.62259e+05,
> atom= 3613
> Step= 63, Dmax= 1.3e-01 nm, Epot= -3.84569e+06 Fmax= 2.41290e+05,
> atom= 3622
> Step= 64, Dmax= 1.5e-01 nm, Epot= -3.84643e+06 Fmax= 2.38963e+06,
> atom= 3622
> Step= 65, Dmax= 1.8e-01 nm, Epot= -3.87919e+06 Fmax= 5.87824e+05,
> atom= 3622
> Step= 66, Dmax= 2.2e-01 nm, Epot= -3.94983e+06 Fmax= 8.25176e+04,
> atom= 3634
>
> step 67: Water molecule starting at atom 170884 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 69, Dmax= 6.6e-02 nm, Epot= -4.00963e+06 Fmax= 4.49243e+05,
> atom= 3607
> Step= 70, Dmax= 8.0e-02 nm, Epot= -4.02678e+06 Fmax= 1.17113e+05,
> atom= 3607
> Step= 71, Dmax= 9.6e-02 nm, Epot= -4.07699e+06 Fmax= 7.52519e+05,
> atom= 3614
>
> step 72: Water molecule starting at atom 212809 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 72, Dmax= 1.1e-01 nm, Epot= -4.07756e+06 Fmax= 1.16158e+06,
> atom= 3614
> Step= 73, Dmax= 1.4e-01 nm, Epot= -4.10951e+06 Fmax= 1.30165e+05,
> atom= 3613
> Step= 74, Dmax= 1.7e-01 nm, Epot= -4.14116e+06 Fmax= 4.07543e+06,
> atom= 3613
>
> step 75: Water molecule starting at atom 212809 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 75, Dmax= 2.0e-01 nm, Epot= -4.18383e+06 Fmax= 8.70248e+05,
> atom= 212809
> Step= 76, Dmax= 2.4e-01 nm, Epot= -4.20933e+06 Fmax= 1.39623e+05,
> atom= 3613
> Step= 78, Dmax= 1.4e-01 nm, Epot= -4.26151e+06 Fmax= 9.65736e+05,
> atom= 3613
>
> step 79: Water molecule starting at atom 212809 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 79, Dmax= 1.7e-01 nm, Epot= -4.26416e+06 Fmax= 9.92323e+05,
> atom= 212809
> Step= 80, Dmax= 2.1e-01 nm, Epot= -4.28752e+06 Fmax= 1.65125e+05,
> atom= 3613
> Step= 81, Dmax= 2.5e-01 nm, Epot= -4.35502e+06 Fmax= 8.28132e+05,
> atom= 3613
>
> step 82: Water molecule starting at atom 212809 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 82, Dmax= 3.0e-01 nm, Epot= -1.93670e+06 Fmax= 2.75398e+08,
> atom= 3607
> step 83: Water molecule starting at atom 212809 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
> Step= 83, Dmax= 1.5e-01 nm, Epot= -4.36885e+06 Fmax= 3.77499e+09,
> atom= 212809
> Step= 84, Dmax= 1.8e-01 nm, Epot= -4.36949e+06 Fmax= 4.70130e+05,
> atom= 3607
> Step= 85, Dmax= 2.1e-01 nm, Epot= -4.39550e+06 Fmax= 6.53132e+05,
> atom= 3607
>
> step 86: Water molecule starting at atom 129571 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
>
>
>
> On Tue, Jun 21, 2011 at 3:15 PM, bharat gupta
> <bharat.85.monu at gmail.com <mailto:bharat.85.monu at gmail.com>> wrote:
>
> I fixed it ... but now after using the command : grompp -f em.mdp
> -c solv.gro -p topol.top -o ions.tpr
>
> I am getting a charge of -9.999971e-01
>
>
> Since the charge has to be in whole number, what shall I do in
> this case. (The ligand that I am using is phosphotyrosine)
>
>
>
> On Tue, Jun 21, 2011 at 2:13 PM, Mark Abraham
> <Mark.Abraham at anu.edu.au <mailto:Mark.Abraham at anu.edu.au>> wrote:
>
> On 21/06/2011 3:03 PM, bharat gupta wrote:
>
> Now I changed the .top file in this way and I am getting
> this error now :-
>
> change :- ; Include forcefield parameters
> #include "charmm27.ff/forcefield.itp"
>
> ;Include ligand topology
> #include "PTR.itp"
>
> [ moleculetype ]
> ; Name nrexcl
> Protein_chain_A 3
>
> =====================
>
> error:-
>
>
> Ignoring obsolete mdp entry 'title'
>
> Back Off! I just backed up mdout.mdp to ./#mdout.mdp.8#
> Generated 23220 of the 23220 non-bonded parameter combinations
> Generating 1-4 interactions: fudge = 1
> Generated 20092 of the 23220 1-4 parameter combinations
> Excluding 3 bonded neighbours molecule type 'Protein_chain_A'
>
> -------------------------------------------------------
> Program grompp, VERSION 4.5.4
> Source code file: toppush.c, line: 1987
>
> Fatal error:
> No such moleculetype LIG
> For more information and tips for troubleshooting, please
> check the GROMACS
> website at http://www.gromacs.org/Documentation/Errors
> -------------------------------------------------------
>
>
> You have to make an attempt to solve your own problems :-) We
> have better things to do than oversee every move you make.
> Follow that link and read about this issue.
>
>
> Mark
> --
> gmx-users mailing list gmx-users at gromacs.org
> <mailto:gmx-users at gromacs.org>
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/Search before
> posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org
> <mailto:gmx-users-request at gromacs.org>.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
>
>
>
> --
> Bharat
> Ph.D. Candidate
> Room No. : 7202A, 2nd Floor
> Biomolecular Engineering Laboratory
> Division of Chemical Engineering and Polymer Science
> Pusan National University
> Busan -609735
> South Korea
> Lab phone no. - +82-51-510-3680, +82-51-583-8343
> Mobile no. - 010-5818-3680
> E-mail : monu46010 at yahoo.com <mailto:monu46010 at yahoo.com>
>
>
>
>
> --
> Bharat
> Ph.D. Candidate
> Room No. : 7202A, 2nd Floor
> Biomolecular Engineering Laboratory
> Division of Chemical Engineering and Polymer Science
> Pusan National University
> Busan -609735
> South Korea
> Lab phone no. - +82-51-510-3680, +82-51-583-8343
> Mobile no. - 010-5818-3680
> E-mail : monu46010 at yahoo.com <mailto:monu46010 at yahoo.com>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110621/7fb79d07/attachment.html>
More information about the gromacs.org_gmx-users
mailing list