AMBER: Charges on atoms in protein :AMBER

From: Varsha Goyal <vgusp.yahoo.ca>
Date: Wed, 1 Feb 2006 22:00:09 -0500 (EST)

Hello,
I generated pdb file after runnig leap:

My leap.in is

source leaprc.gaff
source leaprc.ff99
loadoff ZNC.lib
loadoff HIC.lib
loadoff HIF.lib
loadoff wat.lib
loadAmberParams frcmod.caii
loadAmberParams frcmod.znc
loadAmberParams ligand.frcmod
loadAmberprep HIC.in
loadAmberprep wat.in
loadAmberprep ligand.prepin
loadAmberprep HIF.in
set BZO restype protein
set ZNC restype protein
set HIC restype protein
set HIF restype protein
Prot = loadPdb protein.pdb
charge Prot
addions Prot K+ 1
solvateBox Prot TIP3PBOX 12 iso
saveAmberParm Prot prmtop.top prmcrd.crd
savepdb Prot protein-mod.pdb
check Prot
charge Prot


But when i look at structure in Sybyl7.1 and try to
display charge on residues i see 0.0 charge on every
single atom. Does amber assign charge by itself or do
we assign charges to residues before running leap.

Thanks in advance,
Varsha Gupta


        

        
                
__________________________________________________________
Find your next car at http://autos.yahoo.ca
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Feb 02 2006 - 06:10:08 PST
Custom Search