Hi James,
What you've done seems correct to me; I have prepared an example gist here:
https://gist.github.com/ptosco/40c2530c67d9c0930b8efbc8c92da0be
which indeed shows correct matches.
One observation here is that Open3DALIGN was conceived to generate
alignments for 3D-QSAR; it really doesn't a great job at computing a 3D
similarity score, as it was not made for that purpose. So the attempts
to score an existing overlay "in place" might be disappointing.
Cheers,
p.
On 26/11/2019 17:00, James Davidson wrote:
Dear All (especially Paolo!),
I have a strong suspicion I have already asked this at some point in
the past – so apologies in advance (but I can’t seem to find the answer)…
I am interested in taking an existing overlay of two RDKit molecules
in 3D and scoring the overlay using Open3DAlign scoring scheme (eg
with MMFF atom-types), but **without** trying to optimise the
alignment or score.
I thought setting maxIters=0 in the call to AllChem.GetO3A() would do
the trick (I even tried setting options=3 to “trigger local
optimization”). Eg
o3a = AllChem.GetO3A(prb_mol, ref_mol, maxIters=0, options=3)
o3a.Matches() # Show the matches
But while the options setting certainly changes the matching atoms
(and the score), the matches don’t seem to correspond well to my
starting alignment…
Any advice is greatly appreciated (including, of course, simply
pointing me to the old answer that I am likely missing!)
Kind regards
James
______________________________________________________________________
PLEASE READ - This email is confidential and may be privileged. It is
intended for the named addressee(s) only and access to it by anyone
else is unauthorised. If you are not an addressee, any disclosure or
copying of the contents of this email or any action taken (or not
taken) in reliance on it is unauthorised and may be unlawful. If you
have received this email in error, please notify the sender or
postmas...@vernalis.com. Email is not a secure method of communication
and the Company cannot accept responsibility for the accuracy or
completeness of this message or any attachment(s). Please check this
email for virus infection for which the Company accepts no
responsibility. If verification of this email is sought then please
request a hard copy. Unless otherwise stated, any views or opinions
presented are solely those of the author and do not represent those of
the Company.
Vernalis Limited (Company no. 2304992), Vernalis (R&D) Limited (no.
1985479) and Vernalis Development Limited (no. 2600483)
Granta Park
Great Abington
Cambridge
CB21 6GB, UK
Tel: +44 (0)1223 895 555
______________________________________________________________________
_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
