Previous Thread
Next Thread
Print Thread
How to deal with more than two proteins?
#962 02/21/04 01:59 PM
Joined: Nov 2003
Posts: 76
freesky Offline OP
Forum Member
OP Offline
Forum Member
Joined: Nov 2003
Posts: 76
Hi all,
I want to simulate the folding process of two helix proteins in water, but I do not know how to deal with two proteins in CHARMM. For example, there will be only one set of coordinates from the PDB file for each protein, and how will I locate the two proteins? And also, what about some other problems needed to be noticed. Your suggestions are welcome.

---- State Key Lab of Solid State Microstructure Nanjing University P.R.C
Re: How to deal with more than two proteins?
freesky #963 02/22/04 04:58 PM
Joined: Sep 2003
Posts: 8,479
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,479
In general, folding on the biological or laboratory timescale is beyond the reach of solvated MD simulations using current technology. If the helix is known to associate as a dimer in solution, you may be able to test the relative stability of a series of models using different contact faces between the 2 helices.

Treat each peptide as a separate SEGID, created with independent GENErate commands in CHARMM. Each can be translated and rotated independently by using atom selections with COOR TRANslate and COOR ROTAte commands. For a protein with multiple subunits, it's best to split the PDB file into pieces for each subunit. In this case, because the sequence and initial coordinates for both are the same, one can use duplication commands instead of reading each file separately. A simple example, with comments:

open unit 2 read card name helix.pdb ! OPEN FILE
read sequ pdb unit 2 ! READ AA SEQUENCE
gener A warn ! MAKE FIRST SEGID, A
rewind unit 2 ! REWIND FILE
read coor pdb unit 2 ! READ COORDS
close unit 2 ! CLOSE FILE
hbuild sele atom A * H* end ! BUILD ALL H ATOM COORDS
coor orient ! ALIGN ALONG X AXIS
coor dupl sele segid A end sele segid B end ! COPY COORDS
coor tran ydir 15.0 sele segid A end ! MOVE SEGID A +15 A
coor tran ydir -15.0 sele segid A end ! MOVE SEGID B -15 A

Rick Venable
computational chemist

Moderated by  lennart, rmv 

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.009s Queries: 18 (0.005s) Memory: 0.8939 MB (Peak: 0.9605 MB) Data Comp: Off Server Time: 2020-08-06 16:14:31 UTC
Valid HTML 5 and Valid CSS