AMBER Archive (2007)

Subject: Re: RE: RE: AMBER: PB Bomb in pb_aaradi(): No radius assigned for atom 4474 F3 f

From: emilia wu (emiliawu_at_dicp.ac.cn)
Date: Mon Jan 22 2007 - 18:22:39 CST


Hi, Ray

This time I set MS=0, radiopt=0, and npopt=1, but the job still can not pass.
and the log file of mm_pbsa.pl gives out "
=>> Creating input
    Sander input
    PBSA input

=>> Calculating energy / entropy contributions
    Calc contrib for ../01/thrombin_cda_wat_com.crd.1
        Calc MM/GB/SAS
 bad atom type: f
        /home2/handicp/wuling/amber9/amber9/exe/sander -O -i sander_com.in -o sander_com.1.out -c ../01/thrombin_cda_wat_com.crd.1 -p ../01/protein.prmtop not successful"

and the sander_com.1.out stop at "
Potential function:
     ntf = 1, ntb = 0, igb = 2, nsnb = 99999
     ipol = 0, gbsa = 1, iesp = 0
     dielc = 1.00000, cut = 999.00000, intdiel = 1.00000
     saltcon = 0.00000, offset = 0.09000, gbalpha= 0.80000
     gbbeta = 0.00000, gbgamma = 2.90912, surften = 1.00000
     rdt = 0.00000, rgbmax = 25.00000
     alpb = 0
     scnb = 2.00000, scee = 1.20000

Frozen or restrained atoms:
     ibelly = 0, ntr = 0

Energy minimization:
     maxcyc = 1, ncyc = 0, ntmin = 1
     dx0 = 0.01000, drms = 0.00010
"
how do I deal with it?
Thanks in advance!
        

emilia
======= 2007-01-21 02:50:58 ÄúÔÚÀ´ÐÅÖÐдµÀ£º=======

>Emilia,
>
>I forgot to mention ... if you use radiopt=0, you also need to set npopt=1
>(the default is 2). See PB manual pages for more info.
>
>Best,
>Ray
>
>==========================================
>Ray Luo, Ph.D.
>Dept Molecular Biology & Biochemistry
>University of California, Irvine, CA 92697
>USPS: PO Box 3900
>Phone: (949) 824-9528
>Email: rluo_at_uci.edu
>Web: http://rayl0.bio.uci.edu/
>==========================================
>
>
>-----Original Message-----
>From: owner-amber_at_scripps.edu [mailto:owner-amber_at_scripps.edu] On Behalf Of
>emilia wu
>Sent: Friday, January 19, 2007 5:54 PM
>To: amber_at_scripps.edu
>Subject: Re: RE: AMBER: PB Bomb in pb_aaradi(): No radius assigned for atom
>4474 F3 f
>
>Hi, Ray
>
>Thanks for your reply. I already try radiopt=0, but it still can not pass,
>and
>the output file of pbsa_com.1.out gives out "PB Bomb in pb_read(): use of
>radi other than vdw sigma for
> np solvation dispersion/cavity is not supported!".
>so what should I do?
>
>
>emilia
>
>
>
>======= 2007-01-20 03:58:18 ÄúÔÚÀ´ÐÅÖÐдµÀ£º=======
>
>>Emilia,
>>
>>Since you are using gaff force field, as indicated by lower-case atom
>types,
>>please set radiopt = 0 to use the radii in your prmtop file.
>>
>>Best,
>>Ray
>>
>>==========================================
>>Ray Luo, Ph.D.
>>Dept Molecular Biology & Biochemistry
>>University of California, Irvine, CA 92697
>>USPS: PO Box 3900
>>Phone: (949) 824-9528
>>Email: rluo_at_uci.edu
>>Web: http://rayl0.bio.uci.edu/
>>==========================================
>>
>>
>>-----Original Message-----
>>From: owner-amber_at_scripps.edu [mailto:owner-amber_at_scripps.edu] On Behalf Of
>>emilia wu
>>Sent: Friday, January 19, 2007 5:16 AM
>>To: amber
>>Subject: AMBER: PB Bomb in pb_aaradi(): No radius assigned for atom 4474 F3
>>f
>>
>>Dear amber users,
>>
>> I try to calculate the binding free energy of an inhibitor and a protein
>>with MM/PBSA using Amber9,
>>but it stopped and the output tell me that"PB Bomb in pb_aaradi(): No
>radius
>>assigned for atom
>> 4474 F3 f". and there is a fluorin atom in the inhibitor, so how to fix
>>it, if I should add some
>>radius parameters for fluorin, which file should be modified, and what
>>parameter should I use?
>>
>>thanks in advance!
>>
>>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡emilia wu
>>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡emiliawu_at_dicp.ac.cn
>>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡2007-01-19
>>-----------------------------------------------------------------------
>>The AMBER Mail Reflector
>>To post, send mail to amber_at_scripps.edu
>>To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>>
>>-----------------------------------------------------------------------
>>The AMBER Mail Reflector
>>To post, send mail to amber_at_scripps.edu
>>To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
>= = = = = = = = = = = = = = = = = = = =
>
>
>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡ÖÂ
>Àñ£¡
>
>
>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡emilia wu
>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡emiliawu_at_dicp.ac.cn
>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡2007-01-20
>
>-----------------------------------------------------------------------
>The AMBER Mail Reflector
>To post, send mail to amber_at_scripps.edu
>To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
>-----------------------------------------------------------------------
>The AMBER Mail Reflector
>To post, send mail to amber_at_scripps.edu
>To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu

= = = = = = = = = = = = = = = = = = = =
                        

¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡ÖÂ
Àñ£¡
 
                                 
¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡emilia wu
¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡emiliawu_at_dicp.ac.cn
¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡2007-01-22

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu