[] = OffPrintProtein(vector) OffPrintProtein: Prints the configuration encoded by vector Fixing the lengths of the edges, each angle defines the position of the next residue respective the two previous ones. For reference on the Offline HP model see: -- H. P. Hsu, V. Mehra and P. Grassberger (2003) Structure optimization in an off-lattice protein model. -- Phys Rev E Stat Nonlin Soft Matter Phys. 2003 Sep;68(3 Pt 2):037703. Epub 2003 Sep 30. -- http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=PLEEE8000068000003037703000001&idtype=cvips&gifs=yes INPUTS vector: Sequence of residues (Angles (continuous values) defining the position of the next residue respective the two previous ) Last version 8/26/2008. Roberto Santana (roberto.santana@ehu.es)
0001 function[] = OffPrintProtein(vector) 0002 % [] = OffPrintProtein(vector) 0003 % OffPrintProtein: Prints the configuration encoded by vector 0004 % Fixing the lengths of the edges, each angle defines the position of the next residue respective the two previous ones. 0005 % For reference on the Offline HP model see: 0006 %-- H. P. Hsu, V. Mehra and P. Grassberger (2003) Structure optimization in an off-lattice protein model. 0007 %-- Phys Rev E Stat Nonlin Soft Matter Phys. 2003 Sep;68(3 Pt 2):037703. Epub 2003 Sep 30. 0008 %-- http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=PLEEE8000068000003037703000001&idtype=cvips&gifs=yes 0009 % INPUTS 0010 % vector: Sequence of residues (Angles (continuous values) defining the position of the next residue respective the two previous ) 0011 % 0012 % Last version 8/26/2008. Roberto Santana (roberto.santana@ehu.es) 0013 0014 0015 0016 % InitConf has been initialized with a Fibbonacci sequence by the RunOffLineProtein.m routine 0017 global HPInitConf; 0018 0019 0020 % OffFindPos translates the vector of angles to the positions into the lattice. 0021 [Pos] = OffFindPos(vector); 0022 0023 sizeChain = size(HPInitConf,2) 0024 0025 figure 0026 hold on 0027 0028 for i=1:sizeChain 0029 if(HPInitConf(i) == 0) 0030 plot(Pos(i,1),Pos(i,2),'*'); 0031 else 0032 plot(Pos(i,1),Pos(i,2),'o'); 0033 end 0034 end 0035 plot(Pos(:,1),Pos(:,2),'b-'); 0036