A new version of the CREDO database has been released today comprising 41,691 protein-ligand complexes from the PDB.
New features include:
Variations tables has been added containing mutation information from dbSNP, EnsEMBL, OMIM and COSMICLigands table now contains an ism field to store isomeric SMILES (very useful for heteropeptides)LigandRing and ResidueRing tables now contain an is_hetero_aromatic fieldFragmentProperties now contains all descriptors from OEMolPropFor instructions on how to download and install the CREDO database, please go to the CREDO website.
In [1]: from credopymol import * In [2]: s = StructureAdaptor().fetchByPDB('3BU5') In [3]: s Out[3]: <Structure('3BU5')> In [4]: s.title Out[4]: 'CRYSTAL STRUCTURE OF THE INSULIN RECEPTOR KINASE IN COMPLEX WITH IRS2 KRLB PEPTIDE AND ATP' In [5]: s.load() In [5]: s.getMutations() Out[5]: [<XRef(Residue,OMIM,610549)>, <XRef(Residue,OMIM,125853)>, <XRef(Residue,OMIM,125853)>] In [6]: mutations = s.getMutations() In [7]: mutations[0].description Out[7]: 'G/V Insulin-resistant diabetes mellitus with acanthosis nigricans type A (IRAN type A)'
In [1]: from credopymol import * In [2]: s = StructureAdaptor().fetchByPDB('2Z5X') In [3]: s.load() In [4]: xref = s.getMutations()[0] In [5]: xref Out[5]: <XRef(Residue,dbSNP,rs1803986)> In [6]: xref.description Out[6]: 'M/I NON_SYNONYMOUS_CODING' In [7]: mut_res = xref.getObject() In [8]: mut_res Out[8]: <Residue(445, MET, )> In [9]: mut_res.show() In [10]: s.Ligands[3] Out[10]: Ligand(600, FAD, A) In [11]: s.Ligands[3].showContacts()
The image shows two crystal structures of human CDK2 in complex with inhibitors (1Y8Y, 2C6L). Structural interaction fingerprints (SIFts) were clustered (based on UniProt accession of the protein) and these two ligands were found to be neighbouring leaves on the same node. The structures can be loaded into PyMOL, aligned and the interactions visualised with command below. The resulting PyMOL scene can be downloaded from the bottom of this page.
In [1]: from credopymol import * In [2]: s = StructureAdaptor().fetchByPDB('2P33') In [3]: s.load() In [4]: s.Ligands[0].showContacts() In [5]: s.Ligands[0].showFragmentContactDensity()