Hi guys, could someone who knows chemistry (I don't) tell me if I'm completely off the mark here?
Attached sdf is D-alanine downloaded from pubchem. Also from pubchem,
copied from the same "cid 71080" page is the inchi string:
InChI=1S/C3H7NO2/c1-2(4)3(5)6/h2H,4H2,1H3,(H,5,6)/t2-/m1/s1
(it's also in the file).
My understanding is that by CIP rules L-alanine is S-alanine and
D-alanine is R-alanine. So when I run (python)
mol = rdkit.Chem.SupplierFromFilename( "71080.sdf" ).next()
for atom in mol.GetAtoms() :
print ("%s%d" % (atom.GetSymbol(),(atom.GetIdx() + 1))),
if atom.GetChiralTag() ==\
rdkit.Chem.rdchem.ChiralType.CHI_TETRAHEDRAL_CW : print ": R"
elif atom.GetChiralTag() ==\
rdkit.Chem.rdchem.ChiralType.CHI_TETRAHEDRAL_CCW : print ": S"
elif atom.GetChiralTag() ==\
rdkit.Chem.rdchem.ChiralType.CHI_UNSPECIFIED : print ": N"
I should get a carbon with an "R". What I get is:
O1 : N
O2 : N
N3 : N
C4 : S
C5 : N
C6 : N
-- with or without calling rdkit.Chem.FindMolChiralCenters( mol )
Is this what's supposed to happen?
TIA
--
Dimitri Maziuk
Programmer/sysadmin
BioMagResBank, UW-Madison -- http://www.bmrb.wisc.edu
71080.sdf
Description: application/vnd.kinar
signature.asc
Description: OpenPGP digital signature
------------------------------------------------------------------------------ Slashdot TV. Videos for Nerds. Stuff that Matters. http://pubads.g.doubleclick.net/gampad/clk?id=160591471&iu=/4140/ostg.clktrk
_______________________________________________ Rdkit-discuss mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
