Skip to content

Instantly share code, notes, and snippets.

@aschreyer
Created May 20, 2011 14:21
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save aschreyer/982989 to your computer and use it in GitHub Desktop.
Save aschreyer/982989 to your computer and use it in GitHub Desktop.
Introduction to the OpenEye PostgreSQL cartridge
SELECT openeye.enumerate_tautomers('O=C1NN=Cc2cocc12');
enumerate_tautomers
---------------------------------------
'{O=c1[nH]ncc2c1coc2,Oc1nncc2c1coc2}'
(1 row)
SELECT is_tautomeric('O=C1NN=Cc2cocc12');
is_tautomeric
---------------
t
(1 row)
SELECT match_atom_names(oeb, 'c1ccccc1')
FROM pdbchem.chem_comp_structures
WHERE het_id = 'STI';
match_atom_names
---------------------------------------------------------
'{{C14,C19,C18,C17,C16,C15},{C23,C29,C28,C27,C26,C25}}'
(1 row)
SELECT mcss(array_agg(ism), exhaustive:=true)
FROM pdbchem.chem_comps
WHERE het_id IN ('STI','NIL','PRC');
mcss
--------------------------------
'Cc1ccccc1Nc2nccc(n2)c3cccnc3'
(1 row)
SELECT murcko_scaffold(ism)
FROM pdbchem.chem_comps
WHERE het_id = 'STI';
murcko_scaffold
--------------------------------------------------------------------
'c1cc(cc(c1)NC(=O)c2ccc(cc2)C[NH+]3CC[NH2+]CC3)Nc4nccc(n4)c5cccnc5'
(1 row)
-- Generic Murcko scaffold
SELECT murcko_scaffold(ism, true)
FROM pdbchem.chem_comps
WHERE het_id = 'STI';
murcko_scaffold
---------------------------------------------------------
'CC(CC1CCCC(C1)CC2CCCC(C2)C3CCCCC3)C4CCC(CC4)CC5CCCCC5'
(1 row)
SELECT l1.ligand_id, l2.ligand_id,
rmsd(lm1.oeb, lm2.oeb, false, true)
FROM credo.ligand_molecules lm1, credo.ligand_molecules lm2,
credo.ligands l1, credo.ligands l2
WHERE l1.ligand_id = lm1.ligand_id
AND l2.ligand_id = lm2.ligand_id
AND lm1.ligand_id < lm2.ligand_id
AND l1.ligand_name = 'STI'
AND l2.ligand_name = 'STI'
ligand_id | ligand_id | rmsd
-----------+-----------+-------------------
3103 | 3919 | 0.371631226525847
3103 | 3920 | 0.401457353892573
3103 | 3921 | 0.395985127397863
3103 | 3922 | 0.395539047014552
3103 | 5333 | 2.50825681562345
3919 | 3920 | 0.167355959408002
3919 | 3921 | 0.130774619798445
3919 | 3922 | 0.182510989456054
3919 | 5333 | 2.57370057628475
3920 | 3921 | 0.153267191668418
3920 | 3922 | 0.162330607510651
3920 | 5333 | 2.46849387256747
3921 | 3922 | 0.198698782795329
3921 | 5333 | 2.55045066442146
3922 | 5333 | 2.48307566522211
(15 rows)
SELECT l1.ligand_id, l2.ligand_id,
(rocs(lm1.oeb, lm2.oeb)).combo
FROM credo.ligand_molecules lm1, credo.ligand_molecules lm2,
credo.ligands l1, credo.ligands l2
WHERE l1.ligand_id = lm1.ligand_id
AND l2.ligand_id = lm2.ligand_id
AND lm1.ligand_id < lm2.ligand_id
AND l1.ligand_name = 'STI'
AND l2.ligand_name = 'STI'
ligand_id | ligand_id | combo
-----------+-----------+----------------
3919 | 3921 | 1.97476577759
3919 | 3922 | 1.95278191566
3920 | 3921 | 1.95201945305
3920 | 3922 | 1.94943797588
3919 | 3920 | 1.94234728813
3921 | 3922 | 1.93043136597
3103 | 3919 | 1.74377202988
3103 | 3922 | 1.71013522148
3103 | 3921 | 1.69342803955
3103 | 3920 | 1.68736684322
3103 | 5333 | 0.890247941017
3921 | 5333 | 0.830210804939
3920 | 5333 | 0.830149054527
3919 | 5333 | 0.823508024216
3922 | 5333 | 0.802740335464
(15 rows)
SELECT uni_molecular_rxn('NCCN','[NH2:1]>>[Nh3+:1]');
uni_molecular_rxn
----------------------
'{C(C[NH3+])[NH3+]}'
(1 row)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment