Hi Geoff,

Thank you.  I think you're right that obrotate would be better suited to my 
needs if I removed the SMARTS matching feature.  I can hack the .cpp myself and 
include an option for users to choose whether to use a SMARTS pattern or the 
existing atom numbers in the .xyz file.

Thanks,

- Lee-Ping

----- Original Message -----
From: "Geoffrey Hutchison" <geo...@pitt.edu>
To: "Lee-Ping Wang" <leep...@stanford.edu>
Cc: openbabel-discuss@lists.sourceforge.net
Sent: Friday, April 27, 2012 6:15:22 PM
Subject: Re: [Open Babel] Questionable behavior in obrotate (or babel)

> I believe that 'obrotate' is determining the atom numbers from the SMART 
> pattern, and this becomes a problem when the atom numbers in 'obrotate' don't 
> match the atom numbers from the XYZ generated by babel.  In the first example 
> (isoleucine), there are two explicit hydrogens that obrotate is not counting. 
>  In the second example, there are two identical CH3 groups; 'obrotate' will 
> rotate the second one even though I explicitly specified the first one.

I apologize no one else responded to your message, I've been a bit busy 
finishing off the semester!

> I am able to work around this manually by trial and error and 'guessing' 
> which atom numbers are requested by 'obrotate', but I think it would be much 
> more straightforward if the atom numbers could be taken from the XYZ file and 
> not determined from the SMART pattern.  If the atom numbers must be 
> determined from the SMART pattern, then it should be done in a manner 
> consistent with the numbering scheme in 'babel'.

If I understand, you actually really want a tool which changes the torsion of a 
set of bonded atoms defined by atom numbers you supply -- you really don't care 
about the SMARTS matching feature of obrotate at all, correct?

If so, I'd suggest you just hack obrotate.cpp to your needs. I don't have time 
tonight, but it's not a huge ordeal.

-Geoff

---
Prof. Geoffrey Hutchison
Department of Chemistry
University of Pittsburgh
tel: (412) 648-0492
email: geo...@pitt.edu
web: http://hutchison.chem.pitt.edu/


------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
OpenBabel-discuss mailing list
OpenBabel-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/openbabel-discuss

Reply via email to