Re: [AMBER] NaN VDWAALS Error During Minimization

From: Jatin Kashyap <jk435.njit.edu>
Date: Wed, 6 Jan 2021 00:09:59 -0500

Hi David,

I followed this procedure and I found three couples of atoms which have 0 Ang distance as shown below[1].
And here is my observation: I am feeding the complex[2] and ligand[3] pdb files as described in the Amber tutorial[4].
Everything is fine with input pdb files i.e. no overlap BUT the output complex pdb file[5] that I obtained at the end of the step “Creating topology and coordinate files for Sustiva-RT complex”
of the tutorial is giving my those three couples of atoms which have exactly same coordinates to each other. Prmtop and rst7 files obtained form same step were used for CPPTRAJ instructions[1].



[1]
1 Warning: Atoms 1:LIG_1.N and 20:LIG_1.N3 are close (0.00)
1 Warning: Atoms 4:LIG_1.C and 22:LIG_1.C13 are close (0.00)
1 Warning: Atoms 17:LIG_1.O and 18:LIG_1.O1 are close (0.57)
1 Warning: Atoms 17:LIG_1.O and 23:LIG_1.O3 are close (0.00)
.
.
...

[2]
https://drive.google.com/file/d/1ZcZWthwZgSzE6WcGsMQEMh6WKlGs_l84/view?usp=sharing

[3]
https://drive.google.com/file/d/1Viw3AZKaBZ1JeK6-ZvtXxzduUhEqmQOH/view?usp=sharing

[4]
https://ambermd.org/tutorials/basic/tutorial4b/index.php <https://ambermd.org/tutorials/basic/tutorial4b/index.php>

[5]
https://drive.google.com/file/d/1S8CDlDIRHUr-CloX3N1AcmPeRyMEpGi_/view?usp=sharing
——
Jatin Kashyap
Ph.D. Student
Dr. Dibakar Datta Group
Department of Mechanical and Industrial Engineering
New Jersey Institute of Technology (NJIT)
University Heights
Newark, NJ 07102-1982
Phone- (201)889-5783
Email- jk435.njit.edu

> On Jan 5, 2021, at 10:44 PM, David A Case <david.case.rutgers.edu> wrote:
>
> On Tue, Jan 05, 2021, Jatin Kashyap wrote:
>>
>> I am getting an error[1] upon performing minimization of drug-protein
>> complex by using the AMBER tutorial[2].
>
> You have two atoms on top of one another in your starting structure:
>>
>> NSTEP ENERGY RMS GMAX NAME NUMBER
>> 1 NaN NaN 9.2217E+09 O1 18
>>
>> BOND = 13750.1241 ANGLE = 10976.5308 DIHED = 2934.9414
>> VDWAALS = NaN EEL = Infinity EGB = 60228.2741
>> 1-4 VDW = 381037.7111 1-4 EEL = 14860.5376 RESTRAINT = 0.0000
>
> Dave Cerutti pointed you to mdgx, but (for me) it's much simpler just
> to use the cpptraj "check" command: it will print out short non-bonded
> contacts and bad bonds (the two most common problems with initial
> structures.)
>
> cpptraj <prmtop file>
> trajin <inpcrd file>
> check
> go
>
> ....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
Received on Tue Jan 05 2021 - 21:30:03 PST
Custom Search