AMBER Archive (2005)Subject: Re: AMBER: generate charmm paramater file from amber parm file
From: David E. Konerding DSD staff (dekonerding_at_lbl.gov)
Date: Wed Feb 16 2005 - 12:13:23 CST
HL Eastwood wrote:
>Dear amber users
>
>I have an amber derived pdb file and corresponding amber paramater file. Is it
>possible for me to generate a charmm parameter file from this data?
>
>An entry in the amber archive posted by Roass Walker suggests using ptraj:
>
>"...you would want the following input file for ptraj (for amber to charmm):
>
>trajin inputtrajectory.traj.gz
>trajout outputtrajectory charmm "
>
>I have done this and it generates a binary trajectory file for charmm. However,
>I still require the paramater file. Any suggestions on how to get this?
>
>
While you could write a program to extract the AMBER parameter data from
the prmtop file, and generate a CHARMm input file (after futzing with
all the values to make them more CHARMm-compatible), it would probably
be easier to minimally reformat the
PDB file so that CHARMm or NAMD's psfgen could convert it to a .psf
file. Unfortunately,
I think NAMD's psfgen has slightly different definitions for protein and
nucleic acid termiii,
and it sometimes adds atoms; you have to be very careful in your atom
and residue naming to ensure everything works the same (IE, you end up
with a PSF file that has a one to one mapping
to the coordinate identity and ordering as the prmtop and AMEBR coords)
Here's an example psfgen input file that I used with my PDB file(s)
starting from 1A31 or 1A35.pdb, after I did a bunch of reformatting:
topology top_all27_prot_na.inp
alias residue HIS HSD
alias atom ILE CD1 CD
alias atom PHE OXT OT1
segment topo {
pdb topoI.pdb
}
coordpdb topoI.pdb topo
guesscoord
alias residue A ADE
alias residue G GUA
alias residue C CYT
alias residue T THY
segment dnaC {
first none
last none
pdb dnaC.pdb
}
coordpdb dnaC.pdb dnaC
guesscoord
segment dnaD {
first none
last none
pdb dnaD.pdb
}
coordpdb dnaD.pdb dnaD
guesscoord
writepdb topoI_dna.psf.pdb
writepsf topoI_dna.psf
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|