Previous Thread
Next Thread
Print Thread
Running CHARMM with ligands in proteins
#176 10/21/03 08:40 AM
Joined: Oct 2003
Posts: 40
tjaart Offline OP
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 40
Hi

I used to run CHARMM only on a protein without any ligands in. How do I insert a ligand? I already have a pdb file with the ligand in the correct position. I also made .crd, .psf and .rtf files with InsightII. When CHARMM reads in the pdb sequence will it recognise the ligand or do I have to put in a specific command for CHARMM to go and fetch the psf file during minimization?

Does anyone know of a tutorial which covers this subject?

Thanks!!


Tjaart de Beer
Re: Running CHARMM with ligands in proteins
tjaart #177 10/21/03 12:22 PM
Joined: Sep 2003
Posts: 4,772
Likes: 2
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 4,772
Likes: 2
I am not familiar with the files produced by Insight, but if they are complete .psf and .crd files for your whole system all you need to do is:
* charmm script for protein+ligand
*
read rtf card name ...
read para card name ...
read psf card name my_system.pdf
read coor card name my_system.crd

If you have separate files for protein and ligand (eg, the files you get from Insight are just for the ligand) you need something like this:
* now with protein and ligand in different files
*
read rtf card name protein.rtf
read rtf card name ligand.rtf append
read param card name protein.prm
read param card name ligand. prm append
read psf card name protein.psf
read psf card name ligand.psf append
read coor card name protein.crd
read coor card name ligand.crd append
! use read coor pdb name .... if the coordinate files are in pdb format

The protein and ligand rtf and parameter files have to be compatible with each other (same atom types).


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
Re: Running CHARMM with ligands in proteins
lennart #178 10/21/03 01:22 PM
Joined: Oct 2003
Posts: 40
tjaart Offline OP
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 40
I get the following error when running CHARMM with ligand psf and crd files genrated by InsightII:

CHARMM> generate prot setu

***** ERROR in GENIC ***** Residue '^@^@^@^@ ' was not found.

What does the ***** ERROR in GENIC ***** mean? My ligand is called PTP. How do I make CHARMM accept it?

Thanks!!!!!!


Tjaart de Beer
Re: Running CHARMM with ligands in proteins
tjaart #179 10/21/03 02:08 PM
Joined: Sep 2003
Posts: 4,772
Likes: 2
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 4,772
Likes: 2
The GENIC error you get indicates that there is a residue without a name ('^@^@^@^@' indicates an empty string) in your sequence.
How did you read in the aa-sequence of the protein? Which RTF are you using?




Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
Re: Running CHARMM with ligands in proteins
lennart #180 10/21/03 02:28 PM
Joined: Oct 2003
Posts: 40
tjaart Offline OP
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 40
CHARMM gets the sequence from the pdb file. I am using the RTF for my protein and ligand produced by InsightII. In which file would this empty string be?


Thanks!!!!!


Tjaart de Beer
Re: Running CHARMM with ligands in proteins
lennart #181 10/22/03 11:11 AM
Joined: Sep 2003
Posts: 250
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 250
While on the topic it seems like a good idea to create another forum for posting rtf's and param files concerning "unsusual" and/or non-standard CHARMM residues, although the appropriate forum might be found under "Parameter Set Discussion"...


Moderated by  John Legato, lennart 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.4
(Release build 20200307)
Responsive Width:

PHP: 5.6.33-0+deb8u1 Page Time: 0.008s Queries: 26 (0.003s) Memory: 0.9322 MB (Peak: 1.0328 MB) Data Comp: Off Server Time: 2020-06-05 05:29:36 UTC
Valid HTML 5 and Valid CSS