CHARMM Development Project Forums CHARMM Community Script Archive Sample scripts for basic tasks in CHARMM

Thanks.

Quote:

---

If I had picked any other value you could still have asked the same question, so perhaps the number 30 is not so special.
You can change it to any other value that better suits your mood, or remove that part of the script if you want.

---

Respected Sir,
I have a problem using your stream file. I am trying to generate a truncated octahedron using the stream file set-up-truncated-octahedron.str file. But it is giving me the following output error message.

CHARMM> image byresidue select .not. SOLUTE end

***** LEVEL -2 WARNING FROM <CHARMM> *****
***** IMAGES NEED TO BE PRESENT
******************************************
BOMLEV ( 0) IS REACHED - TERMINATING. WRNLEV IS 5


I dont understand the error. What kind of error is it ? where and how do i get the image.

Please try to guide me. I will be highly thankful to you.

Indrajit Deb.
Research Scholar.
Dept. of Biophysics, molecular biology & bioinformatics.
Calcutta University.
India.

The stream file sets up the PBC with images. Did you modify it, or were there other errors prior to the one you report?

Respected Sir,
Thanks for reply, there are no prior error report.
Do charmm have to read any .img file like tips.img. I am waiting for your help.
Thanks

It would help if you posted more of the output; as I said the script does normally setup PBC and images for you, so something unexpected must have happened before the error you posted. Does your system fulfill the stated assumptions for using the script, eg "PSF and coordinates from a system at least as large as the truncated octahedron have to be present"? One way to get this is to solvate your system in a sphere (solvate-sphere.str).

Hi,

I am a beginner using charmm.
I have a question about running script "gen-prot.inp". When reading sequence from PDB file, why do we use 4pti-a.pdb not 4pti.pdb ? Is it because of the number of peptide chain ?
Then, can we identify the number of peptide chain from the number of chain ID ?

Finally, how do you make 4pti-a.pdb from 4pti.pdb ?
Is there any script that change from 4pti.pdb to 4pti-a.pdb or do I need to change them by hand?

Thank you.

Jun

Yes it is because there is more than one chain in the original file; read the comments in the script.
Yes.
Editor.
There are probably lots scripts out there to do this. Note that you would already be finihsed with this task if you had simply done it instead of having to wait for an answer here.

Hey Lennart,

I am using your solvent-box script, the generated crd file contains 263 water molecules and the psf generated contains 252.

Is that normal or is their something weird that is going on somewhere?

That should not happen.

Hi,
Rick
I am learner about CHARMM
When I the above gen_alpha_helix.inp
I got hte problem ! Unix environment variable CHM_HOME should point to the charmm installation directory
I am using the C43b1
what is the actual problem ?
what sequence should i follow for molecular dynamics simulation of the given above tar file examples
Regards,
ED