[MMTK] atom position error
Michael Sorich
mike_sorich@hotmail.com
Fri, 19 Apr 2002 16:32:52 +0930
This is a multi-part message in MIME format.
------=_NextPart_000_0001_01C1E7BF.D9D3C9F0
Content-Type: text/plain;
charset="Windows-1252"
Content-Transfer-Encoding: quoted-printable
Hi,
=20
I have just started to use the MMTK and I am not sure where I have gone
wrong. I have a small molecule in a pdb file and I am trying to extract
the co-ordinates of each atom. An example of the code I have used is
below:
=20
from MMTK import *
from MMTK.PDB import PDBConfiguration
=20
pdb_conf1 =3D PDBConfiguration("=85\\bilr.pdb")
moleculeCollection =3D pdb_conf1.createMolecules()
bilirubin =3D moleculeCollection.objectList()[0]
=20
print bilirubin.atoms[0].position()
print bilirubin.atoms[1].position()
=20
=20
This outputs:
[-1.1456999999999999, 0.074399999999999994, 0.1062]
[-0.50519999999999987, -0.27289999999999998, 0.58729999999999993]
=20
However the co-ordinates in my pdb file are ten times greater, i.e.
REMARK 4 Bili COMPLIES WITH FORMAT V. 2.0=20
HETATM 1 N1 UNK 0 -11.457 0.744 1.062 1.00 0.00
N =20
HETATM 2 O2 UNK 0 -5.052 -2.729 5.873 1.00 0.00
O
=85
=20
Any help would be appreciated
=20
Michael Sorich
PhD Student
School of Pharmaceutical, Molecular and Biomedical Sciences
University of South Australia
Email: michael.sorich@postgrads.unisa.edu.au
mike_sorich@hotmail.com
=20
---
Outgoing mail is certified Virus Free.
Checked by AVG anti-virus system (http://www.grisoft.com).
Version: 6.0.346 / Virus Database: 194 - Release Date: 10/04/2002
=20
------=_NextPart_000_0001_01C1E7BF.D9D3C9F0
Content-Type: text/html;
charset="Windows-1252"
Content-Transfer-Encoding: quoted-printable
<html xmlns:o=3D"urn:schemas-microsoft-com:office:office" =
xmlns:w=3D"urn:schemas-microsoft-com:office:word" =
xmlns:st1=3D"urn:schemas-microsoft-com:office:smarttags" =
xmlns=3D"http://www.w3.org/TR/REC-html40">
<head>
<META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
charset=3DWindows-1252">
<meta name=3DProgId content=3DWord.Document>
<meta name=3DGenerator content=3D"Microsoft Word 10">
<meta name=3DOriginator content=3D"Microsoft Word 10">
<link rel=3DFile-List href=3D"cid:filelist.xml@01C1E7BF.D9164650">
<o:SmartTagType =
namespaceuri=3D"urn:schemas-microsoft-com:office:smarttags"
name=3D"PlaceName"/>
<o:SmartTagType =
namespaceuri=3D"urn:schemas-microsoft-com:office:smarttags"
name=3D"PlaceType"/>
<o:SmartTagType =
namespaceuri=3D"urn:schemas-microsoft-com:office:smarttags"
name=3D"place"/>
<!--[if gte mso 9]><xml>
<o:OfficeDocumentSettings>
<o:DoNotRelyOnCSS/>
</o:OfficeDocumentSettings>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:WordDocument>
<w:SpellingState>Clean</w:SpellingState>
<w:GrammarState>Clean</w:GrammarState>
<w:DocumentKind>DocumentEmail</w:DocumentKind>
<w:EnvelopeVis/>
<w:Compatibility>
<w:BreakWrappedTables/>
<w:SnapToGridInCell/>
<w:WrapTextWithPunct/>
<w:UseAsianBreakRules/>
</w:Compatibility>
<w:BrowserLevel>MicrosoftInternetExplorer4</w:BrowserLevel>
</w:WordDocument>
</xml><![endif]--><!--[if !mso]>
<style>
st1\:*{behavior:url(#default#ieooui) }
</style>
<![endif]-->
<style>
<!--
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{mso-style-parent:"";
margin:0cm;
margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:12.0pt;
font-family:"Times New Roman";
mso-fareast-font-family:"Times New Roman";}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline;
text-underline:single;}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline;
text-underline:single;}
span.EmailStyle17
{mso-style-type:personal-compose;
mso-style-noshow:yes;
mso-ansi-font-size:10.0pt;
mso-bidi-font-size:10.0pt;
font-family:Arial;
mso-ascii-font-family:Arial;
mso-hansi-font-family:Arial;
mso-bidi-font-family:Arial;
color:windowtext;}
span.SpellE
{mso-style-name:"";
mso-spl-e:yes;}
span.GramE
{mso-style-name:"";
mso-gram-e:yes;}
@page Section1
{size:595.3pt 841.9pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;
mso-header-margin:35.4pt;
mso-footer-margin:35.4pt;
mso-paper-source:0;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 10]>
<style>
/* Style Definitions */=20
table.MsoNormalTable
{mso-style-name:"Table Normal";
mso-tstyle-rowband-size:0;
mso-tstyle-colband-size:0;
mso-style-noshow:yes;
mso-style-parent:"";
mso-padding-alt:0cm 5.4pt 0cm 5.4pt;
mso-para-margin:0cm;
mso-para-margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:10.0pt;
font-family:"Times New Roman";}
</style>
<![endif]-->
</head>
<body lang=3DEN-AU link=3Dblue vlink=3Dpurple =
style=3D'tab-interval:36.0pt'>
<div class=3DSection1>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Hi,<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>I have just started to use the MMTK and I am not sure =
where
I have gone wrong. I have a small molecule in a <span =
class=3DSpellE>pdb</span>
file and I am trying to extract the co-ordinates of each atom. An =
example of
the code I have used is below:<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>from</span></font></span><fo=
nt
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> MMTK import
*<o:p></o:p></span></font></p>
<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>from</span></font></span><fo=
nt
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> MMTK.PDB
import <span =
class=3DSpellE>PDBConfiguration</span><o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>pdb_conf1 =3D <span class=3DSpellE><span =
class=3DGramE>PDBConfiguration</span></span><span
class=3DGramE>(</span>"=85\\bilr.pdb")<o:p></o:p></span></font>=
</p>
<p class=3DMsoNormal><span class=3DSpellE><span class=3DGramE><font =
size=3D2
face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'>moleculeCollection</span></f=
ont></span></span><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> =3D
pdb_conf1.createMolecules()<o:p></o:p></span></font></p>
<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>bilirubin</span></font></spa=
n><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> =3D <span
class=3DSpellE>moleculeCollection.objectList</span>()[0]<o:p></o:p></span=
></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>print</span></font></span><f=
ont
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> <span
class=3DSpellE>bilirubin.atoms</span>[0].position()<o:p></o:p></span></fo=
nt></p>
<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>print</span></font></span><f=
ont
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> <span
class=3DSpellE>bilirubin.atoms</span>[1].position()<o:p></o:p></span></fo=
nt></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>This =
outputs</span></font></span><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'>:<o:p></o:p></span></font></=
p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>[-1.1456999999999999, 0.074399999999999994, =
0.1062]<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>[-0.50519999999999987, -0.27289999999999998,
0.58729999999999993]<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>However the co-ordinates in my <span =
class=3DSpellE>pdb</span>
file are ten times greater, i.e.<o:p></o:p></span></font></p>
<p class=3DMsoNormal =
style=3D'mso-layout-grid-align:none;text-autospace:none'><font
size=3D2 face=3D"Courier New"><span =
style=3D'font-size:10.0pt;font-family:"Courier New"'>REMARK<span
style=3D'mso-spacerun:yes'> </span>4 <span =
class=3DSpellE>Bili</span>
COMPLIES WITH FORMAT V. 2.0 <o:p></o:p></span></font></p>
<p class=3DMsoNormal =
style=3D'mso-layout-grid-align:none;text-autospace:none'><font
size=3D2 face=3D"Courier New"><span =
style=3D'font-size:10.0pt;font-family:"Courier New"'>HETATM<span
style=3D'mso-spacerun:yes'> </span><span =
class=3DGramE>1<span
style=3D'mso-spacerun:yes'> </span>N1</span><span
style=3D'mso-spacerun:yes'> </span>UNK<span
style=3D'mso-spacerun:yes'> </span>0<span
style=3D'mso-spacerun:yes'> </span>-11.457<span
style=3D'mso-spacerun:yes'> </span>0.744<span
style=3D'mso-spacerun:yes'> </span>1.062<span
style=3D'mso-spacerun:yes'> </span>1.00<span
style=3D'mso-spacerun:yes'> </span>0.00<span
style=3D'mso-spacerun:yes'> &nbs=
p;
</span>N<span style=3D'mso-spacerun:yes'> =
</span><o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3D"Courier New"><span =
style=3D'font-size:10.0pt;
font-family:"Courier New"'>HETATM<span
style=3D'mso-spacerun:yes'> </span><span =
class=3DGramE>2<span
style=3D'mso-spacerun:yes'> </span>O2</span><span
style=3D'mso-spacerun:yes'> </span>UNK<span
style=3D'mso-spacerun:yes'> </span>0<span
style=3D'mso-spacerun:yes'> =
</span>-5.052<span
style=3D'mso-spacerun:yes'> </span>-2.729<span
style=3D'mso-spacerun:yes'> </span>5.873<span
style=3D'mso-spacerun:yes'> </span>1.00<span
style=3D'mso-spacerun:yes'> </span>0.00<span
style=3D'mso-spacerun:yes'> &nbs=
p;
</span>O<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3D"Courier New"><span =
style=3D'font-size:10.0pt;
font-family:"Courier New"'>=85</span></font><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'><o:p></o:p></span></font></p=
>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Any help would be =
appreciated<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial;mso-no-proof:yes'>Michael =
Sorich<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial;mso-no-proof:yes'>PhD =
Student<o:p></o:p></span></font></p>
<p class=3DMsoNormal><st1:place><st1:PlaceType><font size=3D2 =
face=3DArial><span
=
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:yes'>School</spa=
n></font></st1:PlaceType><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:
yes'> of </span></font><st1:PlaceName><font size=3D2 face=3DArial><span
=
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:yes'>Pharmaceuti=
cal</span></font></st1:PlaceName></st1:place><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:
yes'>, Molecular and Biomedical Sciences<o:p></o:p></span></font></p>
<p class=3DMsoNormal><st1:place><st1:PlaceType><font size=3D2 =
face=3DArial><span
=
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:yes'>University<=
/span></font></st1:PlaceType><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:
yes'> of </span></font><st1:PlaceName><font size=3D2 face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:yes'>South =
Australia</span></font></st1:PlaceName></st1:place><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial;mso-no-proof:
yes'><o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial;mso-no-proof:yes'>Email:
michael.sorich@postgrads.unisa.edu.au<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial;mso-no-proof:yes'><span
style=3D'mso-spacerun:yes'> &nbs=
p;
</span>mike_sorich@hotmail.com<o:p></o:p></span></font></p>
<p class=3DMsoNormal><font size=3D3 face=3D"Times New Roman"><span =
style=3D'font-size:
12.0pt'><o:p> </o:p></span></font></p>
</div>
</body>
</html>
<BR>
<P><FONT SIZE=3D2>---<BR>
Outgoing mail is certified Virus Free.<BR>
Checked by AVG anti-virus system (http://www.grisoft.com).<BR>
Version: 6.0.346 / Virus Database: 194 - Release Date: 10/04/2002<BR>
</FONT> </P>
------=_NextPart_000_0001_01C1E7BF.D9D3C9F0--