<div>Thanks Yang Ye! I will look into this further. Thanks for the program, I'm sure that will save me much time.</div>
<div> </div>
<div>Russell Green<br><br> </div>
<div><span class="gmail_quote">On 6/11/07, <b class="gmail_sendername">Yang Ye</b> <<a href="mailto:leafyoung@yahoo.com">leafyoung@yahoo.com</a>> wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div>
<div style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">
<div style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman,new york,times,serif">As what we normally do with dPCA, after obtaining the 2D projection of first and second principle, we make it into a 3D plot of free energy surface. The awk program is provided below to do such conversion. On this surface, we can located regions with low free energy and this leads us to go back to 1D projection to locate which frames are inside the regions. Clustering or averaging may apply to all the frames found in one region. So, in the end, structure are found inside the trajectory for closer inspection.
<br><br>BEGIN {<br> i=0<br> grid_num=40.0<br><br> x_max=0.0<br> y_max=0.0<br> x_min=0.0<br> y_min=0.0<br><br> temp=300<br> dkt=8.31/1000.0*temp/4.182<br>}<br><br>{<br> x[i]=$1<br> y[i]=$2<br><br> if ($1>x_max) x_max=$1
<br> if ($2>y_max) y_max=$2<br> if ($1<x_min) x_min=$1<br> if ($2<y_min) y_min=$2<br> i++<br>}<br><br>END {<br> pt_num=i<br><br> grid_x_size=(x_max-x_min)/grid_num<br> grid_y_size=(y_max-y_min)/grid_num
<br><br> dm=0<br> for(i=0;i<pt_num;i++) {<br> grid[(int((x[i]-x_min)/grid_x_size)), (int((y[i]-y_min)/grid_y_size))]++<br> g=grid[(int((x[i]-x_min)/grid_x_size)), (int((y[i]-y_min)/grid_y_size))]<br> if(dm<g) { dm=g; }
<br> }<br><br> mm=0;<br> for(i=0;i<grid_num;i++) {<br> for(j=0;j<grid_num;j++) {<br> if(grid[i,j]!="") { grid[i,j]=dkt*log(dm/grid[i,j]); if (grid[i,j]>mm) {mm=grid[i,j];} }<br> }
<br> }<br><br> for(i=0;i<grid_num;i++) {<br> for(j=0;j<grid_num;j++) {<br> if(grid[i,j]=="") { grid[i,j]=mm+0.01 }<br> print x_min+(i+0.5)*grid_x_size " " y_min+(j+0.5)*grid_y_size " " grid[i,j]
<br> }<br> }<br>}<br><br>
<div style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman,new york,times,serif">
<div><span class="e" id="q_1131c12dcf33a710_1">----- Original Message ----<br>From: Russell Green <<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:russellwgreen@gmail.com" target="_blank">russellwgreen@gmail.com
</a>><br>To: Discussion list for GROMACS users <<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>Sent: Tuesday, June 12, 2007 2:17:45 AM
<br>Subject: [gmx-users] dPCA "atoms"<br><br>
<div>Dear Gromacs community,</div>
<div> </div>
<div> I've been using the PCA analysis in dihedral space for a chain of my protein of interest. I made a script to collect the atom quadruplets corresponding to the phi and psi angles and output them to an index file. I followed the steps in the paper manual and everything worked as it should; however, is there any way to relate the dihedral PCA information to each residue that contributed to the analysis? As it stands, my dummy reference file contains 715 "atoms" but the analysis was done for about 536 residues so that's 3 atoms per residue helping to define the dihedral angles. Could someone please provide some suggestions or
<font style="BACKGROUND-COLOR: rgb(255,255,255)">more insight? I've searched the net and the mailing list but have found no information regarding my question.</font></div>
<div> </div>
<div>Thanks in advance,</div>
<div>Russell Green</div>
<div>UAB Undergraduate Researcher</div></span></div>
<div>_______________________________________________<br>gmx-users mailing list <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/search" target="_blank">
http://www.gromacs.org/search</a> before posting!<br>Please don't post (un)subscribe requests to the list. Use the <br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org" target="_blank">
gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php
</a></div></div><br> </div></div></div><br>_______________________________________________<br>gmx-users mailing list <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org
</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at
<a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>Please don't post (un)subscribe requests to the list. Use the
<br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">
http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote></div><br>