<div>I believe BioPython project can help you. It load a PDB file and you have methods which help you selecting atoms.</div><div><br>Please see [1]...</div><div> </div><div>[1] <a href="http://biopython.org/wiki/Main_Page">http://biopython.org/wiki/Main_Page</a></div>
<div> </div><div>Best regards,</div><div> </div><div>--<br>Rodrigo Antonio Faccioli<br>Ph.D Student in Electrical Engineering<br>University of Sao Paulo - USP<br>Engineering School of Sao Carlos - EESC<br>Department of Electrical Engineering - SEL<br>
Intelligent System in Structure Bioinformatics<br><a href="http://laips.sel.eesc.usp.br" target="_blank">http://laips.sel.eesc.usp.br</a><br>Phone: 55 (16) 3373-9366 Ext 229<br>Curriculum Lattes - <a href="http://lattes.cnpq.br/1025157978990218" target="_blank">http://lattes.cnpq.br/1025157978990218</a><br>
Public Profile - <a href="http://br.linkedin.com/pub/rodrigo-faccioli/7/589/a5" target="_blank">http://br.linkedin.com/pub/rodrigo-faccioli/7/589/a5</a><br>
<br><br></div><div class="gmail_quote">On Thu, Apr 19, 2012 at 5:25 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br><blockquote style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid" class="gmail_quote">
<div class="im"><br>
<br>
patrick wintrode wrote:<br>
<blockquote style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid" class="gmail_quote">
Hi.<br>
<br>
For each alpha carbon in my protein, I want to write out the x, y and z coordinates as separate time series.<br>
<br>
If I use g_traj with the flags -oxt -x (or y orz) and -n along with an index file selecting the appropriate alpha carbon, will that do the trick? Does anyone know of a less cumbersome way of doing this?<br>
<br>
</blockquote>
<br></div>
You don't necessarily have to do each alpha carbon separately (though you could script all of it, which is relatively easy). If you use the default C-alpha group, all of the chosen x/y/z coordinates are printed in order of the atom number, which can then be parsed with something like awk or perl.<br>
<br>
-Justin<br>
<br>
-- <br>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" target="_blank" value="+15402319080">(540) 231-9080</a><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<span class="HOEnZb"><font color="#888888"><br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</font></span></blockquote></div><br>