Hi Claire,
My thought is that you are getting the core dump before your script file
ever gets down to the strip line. I would think that it would core dump the
first time you input a trajin line that doesn't coordinate to your prmtop
file. You can test this theory by inputting the lines directly into ptraj
rather than using a script file. If it blows up when you trajin a line,
then there's your problem.
If that's the case I would try this. For your overall average, it seems
like you only care about creating a structure using the Calpha atoms of the
protein and not any of the ligands. What I would do is modify your five
average structures such that they only contain the Calpha atoms (doing so
one by one in ptraj) and then create a new topology which also contained
only the Calpha atom information. Then, you should be able to rerun ptraj
with these five new pdbs and one new top file.
Hope this helps!
John
>From: "Claire Zerafa" <cporteli.waldonet.net.mt>
>Reply-To: amber.scripps.edu
>To: <amber.scripps.edu>
>Subject: AMBER: Average pdb
>Date: Thu, 23 Mar 2006 22:17:33 +0100 (W. Europe Standard Time)
>
>Hi all,
>
>I have a question please. I am trying to fit the trajectory of a protein
>bound to 5 different ligands to an average pdb file. For each protein
>ligand
>system, I have an average pdb file for the trajectory. Now I try to read
>each pdb file into ptraj and make an average pdb file out of those so that
>I
>can fix each trajectory to the average pdb file to make it more
>representative of my system. My first question is as follows:
>
>For each of my 5 systems I have a prmtop file. When I come to run ptraj to
>generate the average pdb file do I just pick one of the prmtop files? Or is
>there a way of generating an average prmtop file?
>
>Second I created a script file as follows:
>
>Trajin /DATA/CORTICOSTERONE_pca_CA_only/corticosterone_average_output.pdb
>Trajin
>/DATA/DIETHYLSTILBOESTROL_pca_CA_only/diethylstilboestrol_ptraj1_average_outp
>t.pdb
>Trajin /DATA/METRIBOLONE_pca_CA_only/metribolone_average_output.pdb
>Trajin /DATA/UNDOCKED_pca_CA_only/UNDOCKED_average_output_ca.pdb
>Trajin /DATA/VINCLOZOLIN_pca_CA_only/vinclozolin_average_output.pdb
>Rms first out /DATA/all_av_rms1.out :.CA (GENERATES rms for the CA atoms
>only)
>Strip ~:.CA (STRIPS ALL BUT CA atoms- and therefore does away with any
>differences caused by presence of different ligands in pdb files)
>Average /DATA/all_av.pdb pdb nowrap (generates a pdb file containing only
>the CA atoms of the protein)
>Go
>
>I get a core dump. I can see why this happens ptraj expects a certain
>number
>of atoms- specified by the prmtop file, and then finds a different amount-
>due to the different ligand size. I would have thought though, that since I
>am working only with the CA atoms of the protein, then this would not have
>mattered. How can I solve this please?
>
>Claire
><< BackGrnd.jpg >>
><< imstp_chubbi_en.gif >>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Mar 26 2006 - 06:10:16 PST
