CYPSI: a structure-based interface for cytochrome P450s and ligands in Arabidopsis thaliana
© Zhang et al.; licensee BioMed Central Ltd. 2012
Received: 20 May 2012
Accepted: 18 December 2012
Published: 20 December 2012
The cytochrome P450 (CYP) superfamily enables terrestrial plants to adapt to harsh environments. CYPs are key enzymes involved in a wide range of metabolic pathways. It is particularly useful to be able to analyse the three-dimensional (3D) structure when investigating the interactions between CYPs and their substrates. However, only two plant CYP structures have been resolved. In addition, no currently available databases contain structural information on plant CYPs and ligands. Fortunately, the 3D structure of CYPs is highly conserved and this has made it possible to obtain structural information from template-based modelling (TBM).
The CYP Structure Interface (CYPSI) is a platform for CYP studies. CYPSI integrated the 3D structures for 266 A. thaliana CYPs predicted by three TBM methods: BMCD, which we developed specifically for CYP TBM; and two well-known web-servers, MUSTER and I-TASSER. After careful template selection and optimization, the models built by BMCD were accurate enough for practical application, which we demonstrated using a docking example aimed at searching for the CYPs responsible for ABA 8′-hydroxylation. CYPSI also provides extensive resources for A. thaliana CYP structure and function studies, including 400 PDB entries for solved CYPs, 48 metabolic pathways associated with A. thaliana CYPs, 232 reported CYP ligands and 18 A. thaliana CYPs docked with ligands (61 complexes in total). In addition, CYPSI also includes the ability to search for similar sequences and chemicals.
CYPSI provides comprehensive structure and function information for A. thaliana CYPs, which should facilitate investigations into the interactions between CYPs and their substrates. CYPSI has a user-friendly interface, which is available athttp://bioinfo.cau.edu.cn/CYPSI.
KeywordsCytochrome P450 CYP Structure Interface Template-based modelling BMCD ABA 8′-hydroxylation CYP707A
Cytochrome P450s (CYPs) are heme containing monooxygenases and are found in all eukaryotes. They catalyse various chemical reactions, e.g. hydroxylations, epoxidations, ring extensions and carbon-carbon bond cleavages, and have potential pharmacological and agronomic applications[1–4]. In terrestrial plants, CYPs play important roles in response to biotic and abiotic stimuli by metabolizing a wide range of small organic compounds[5–8]. CYPs are also involved in the biosynthesis of many structural components[9–13].
The three-dimensional (3D) structures of CYPs may provide valuable information that could be used to investigate the interactions between CYPs and ligands. To date, there are more than 5,100 annotated plant CYPs sequences[3, 14], but only two have resolved 3D structures (CYP74A and CYP74A2)[15, 16]. CYP structures are difficult to determine by standard X-ray or NMR analysis because most of them are membrane-bound proteins. Template-based modeling (TBM) could be a feasible alternative method for obtaining CYP structure information because the 3D structure is highly conserved. There are many choices for CYP TBM, e.g. the class-dependent sequence alignment strategy for CYP TBM, SWISS-MODEL, MUSTER and I-TASSER. I-TASSER was found to be the most accurate in a recent Critical Assessment of Techniques for Protein Structure Prediction (CASP7-9)[21–23]. However, the models generated by these web-servers have no heme, the position of which is important when investigating the interaction between CYPs and substrates. We developed a pipeline BMCD specifically for CYP TBM (abbreviation of the softwares used: PSI-B LAST, M USCLE, C OMPASS and D iscovery Studio 2.1)[24–27].
Most current CYP related resources focus on gene annotation, e.g. the Cytochrome P450 Homepage, the CYP engineering database (CYPED) and the Fungal Cytochrome P450 Database (FCPD). Although some databases collect CYP structure information, e.g. CYPED presents all available 3D CYP structures from the Protein Data Bank and SuperCYP collect many drug-drug interactions and the theoretical models for human CYPs, neither of them provide further information about the interactions between ligands and CYPs[29, 31].
Our study has developed the CYP Structure Interface (CYPSI), a platform that provides comprehensive structure and function information on all 266 A. thaliana CYPs. The models for these CYPs were predicted using the BMCD pipeline and the web-servers: MUSTER and I-TASSER. CYPSI also provides extensive resources for CYPs, including 400 PDB entries for solved CYPs, 48 metabolic pathways associated with A. thaliana CYPs, 232 reported CYPs ligands and 18 A. thaliana CYPs docked with ligands (61 complexes in total). To demonstrate the quality and utility of the 3D structures in CYPSI, this paper discusses a case study which searches for the candidate CYPs responsible for abscisic acid (ABA) 8′-hydroxylation. With the implementation of sequence alignment, the BMCD service for template selection and a structure similarity search facility for small molecules, CYPSI is a comprehensive tool for the investigation of plant CYP structures and functions.
Construction and Content
The solved CYP structures were collected from the Protein Data Bank (http://www.rcsb.org/). Up to December 2011, there were 400 PDB entries associated with 76 CYPs (see Additional file1 for details).
A total of 290 A. thaliana CYPs isoforms from 272 CYP genes distributed in 47 CYP families were collected from TAIR10 (http://www.arabidopsis.org/) andhttp://www.p450.kvl.dk/[33, 34], including functional annotations, protein sequences, c od ing s equences (CDS) and 3,000 base pairs (bp) upstream and 3,000 bp downstream of the CDS.
In addition, 48 metabolic pathways were manually collected from the PMN database and from the scientific literature. Pathways clarified in the scientific literature were marked with “y”, and those that had not been clarified were marked with “na” (see Additional file2). A total of 232 ligands in these pathways were collected from PubChem or built manually by Discovery Studio 2.1.
BMCD was specifically developed for CYP TBMs, with an emphasis on template selection and sequence alignment. First, profile-profile alignments between the sequence profiles of targets and templates were constructed using COMPASS. Next, the five templates with the smallest evolutional distances (ED) were selected for further TBM. ED was calculated as described in reference using the substitution score matrix, MIYS960102. Finally, for each target-template pair, three initial models were built using MODELLER in Discovery Studio 2.1 (Accelrys Software Inc.), using the coenzyme heme copied from the template. Of the 15 initial models created for each target, the one with the highest Profiles-3D score was retained for further refinement.
The CHARMm force field in Discovery Studio 2.1 was used by this project for all processes, including energy minimization, molecular dynamic (MD) simulation, the docking program (CDOCKER), and for interaction energy calculations (See Additional file1 for details).
Besides the BMCD, two servers: MUSTER and I-TASSER, were also used for A. thaliana CYP model generation by submitting the sequences for A. thaliana CYPs manually. The prediction results indicated that out of the 279 A. thaliana CYPs longer than 300 amino acids, 266 CYPs would have complete CYP structural domains.
A practical application of CYP 3D model
In order to demonstrate the usefulness of the CYP 3D models, a practical application search for CYPs responsible for ABA 8′-hydroxylation is presented below.
Firstly, ABA was docked to all nine CYPs candidates (11 models) proposed by Eiji Nambara et al. to be responsible for ABA 8′-hydroxylation using CDOCKER. CYP97A3, CYP97B3, CYP97C1 and CYP714A1 were excluded from further analysis because they could not bind ABA and form a suitable conformation for hydroxylation, as determined by our docking result [data not shown].
The Interaction energy (between ABA and receptors) and the Distance (between ABA C8′ and Fe) before and after MD simulation
Interaction Energy (−kcal/mol)
Key binding residues
K78, F88, F243 and L312
K78, F88, F243 and L312
K78, F88, F249 and L319
K78, F88, F243 and L312
K78, F88, F243 and L312
Y74, K78, F88, F248 and I318
Secondly, energy minimization and MD simulation were performed on the seven candidate docking complexes. We compared changes in the ABA locations in these complexes before and after MD simulation. The location of ABA in ABA-CYP704A2 changed considerably compared to ABA-CYP707A, which indicated that this complex was not stable (Table 1, Figure 3 and Additional file3).
The interaction energy between ABA and CYPs decreased significantly after energy minimization or MD simulation, which indicated that these steps were necessary if a more reliable complex was to be obtained because a lower interaction energy represents firmer binding. It should also be noted that the interaction energy for ABA-CYP704A2 was much higher than that of ABA-CYP707As (Table 1). Integration of the above results, including the binding sites, ABA location and the interaction energy, supported the hypothesis that CYP704A2 is unlikely to be ABA 8′-hydroxylase.
CYP707A4 had the lowest catalytic activity for ABA 8′-hydroxylation among the four CYP707As. Intriguingly, after MD simulation, a hydrogen bond was formed between the Tyr74 of CYP707A4 and ABA, which did not occur with the other CYP707As (Figure 2 and Additional file3), possibly because the equivalent residues for the other CYP707As were different from CYP707A4. For example, the 74th residue is Phe for CYP707A1 and CYP707A3. The residue and hydrogen bond differences at the 74th site indicated a lower catalytic activity for CYP707A4 during ABA 8′-hydroxylation, which is consistent with previous results.
In summary, the docking results suggested that many potential CYPs and key residues should be prioritised for further validation studies (Table 1) and that the results have provided valuable insights into the mechanism behind ABA 8′-hydroxylation that need further investigation.
CYPSI database construction
Hyperlinks to PDB, TAIR, UniProt and PubMed are provided. Some useful tools are also integrated into CYPSI to facilitate the browsing and search functions, including sequence alignment, a search function for chemicals with a similar structure and 3D structure animation using Jmol.
Solved CYPs structures
A. thaliana CYPs models
Besides the ability to browse the data shown above, CYPSI also provides three search capabilities: by keywords, by chemical structures and by protein sequences.
From the search box located at the upper right hand corner of the web-page, users can search for information using the keywords: Arabidopsis Genome Initiative (AGI), PDB IDs, CYP families and pathways.
In CYPSI, the BMCD server is used for template selection and sequence alignment (Additional file7). Users only need to submit the target CYP sequence and the results will feedback in a few minutes. The sequence alignments given by BMCD can be utilized directly by Discovery Studio 2.1 for TBM.
To facilitate the study of plant CYPs, we have constructed the CYPSI platform, which contains comprehensive information on CYP sequences, structures, ligands and functions. Notably, all A. thaliana CYP 3D models were predicted using the BMCD pipeline and preliminary refinements have been made, which is particularly useful when investigating CYP structures and functions. In general, there are four steps involved in TBM: template selection and sequence alignment, model construction, model refinement and model validation.
The quality of the template is a key factor that determines the quality of the predicted models. Prior to TBM, a potential template was carefully selected, taking into consideration the completeness of the structure, resolution, presence of a substrate, and the Profile-3D Score.
CYP sequences are highly diverse and it is hard to find the most suitable template and obtain the correct sequence alignment for TBM[1, 17, 47]. We developed BMCD for CYP structure modelling and used the profile-profile alignment by COMPASS and ED to evaluate the similarities between templates and targets so that the best template is selected. In addition, most models generated by BMCD are based on a single template as multiple templates may result in considerable structural errors[21, 23].
The recommended BMCD models need further refinement, which is even more difficult to control than template selection and sequence alignment[18, 21]. Energy minimization and MD simulation are the main methods used for molecular refinement. However, in general, it is difficult to improve the accuracy of the models using these methods as the force fields utilized at present are not accurate enough. For example, in the case of CYP74A modelling (Additional file8), I-TASSER utilized a special force field to refine the models. However, it performed even worse than MUSTER in terms of RMSD and TM-score. We found that many models, following energy minimization, were worse than the initial BMCD models, as evaluated by the Profile-3D ratio. Therefore, we only refined the residues around the coenzyme heme, which is essential for the study of CYP and ligand interactions.
Despite there being many defects in the field of structure modelling, the CYP models in CYPSI could still be very useful for experimental researchers. In the practical application case study, which searched for CYPs responsible for ABA 8-hydroxylation, although the sequence identities of the CYP707A-template pairs were around 30%, which is theoretically too low to build a high-quality homology model, the docking and MD simulation results coincided well with previous experimental results. These results also identified potential residues for ABA binding, which should help reveal the possible catalytic mechanism involved. However, conformational errors in these models are inevitable. Residues that are close to a ligand may affect the final docking result, so softwares that can cope with both ligand and protein flex are recommended for ligand docking, e.g. AutoDock. Further energy minimization or MD simulation methods are recommended so that more comprehensive and reliable information about the enzyme-ligand complex can be obtained.
CYPSI was constructed as a comprehensive platform, integrating sequences, structures, ligands and functional information for CYPs. In addition, it also provides useful tools and resources for CYP structural and functional investigations. The recommended models in CYPSI could be used directly for substrate docking and these enzyme-ligand complexes could provide valuable insights for experimental scientists. Further development of CYPSI will lead to the identification of more enzyme-ligand complexes.
Availability and requirements
The database is available athttp://bioinfo.cau.edu.cn/CYPSI, which is compatible with most modern web browsers. All the data in CYPSI are downloadable and freely available to the academic community.
CYPs structure interface
PSI-B LAST, M USCLE, C OMPASS, D iscovery Studio 2.1
- MD simulation:
molecular dynamic simulation
We would like to thank Ms. Wenying Xu’s critical suggestions on the CYPSI design, Dr. Zi-ding Zhang’s advice and technical support from Dr. Yi Ling. We also thank Dr. Wei-xuan Wang, Dr. Guang-bin Zhang and Dr. You-song Peng’s review and opinions on the manuscript. This work was supported by grants from the Ministry of Science and Technology of China (31171276 and 30570139).
- Rupasinghe S, Schuler MA: Homology modeling of plant cytochrome P450s. Phytochemistry Reviews 2006, 473–505.Google Scholar
- Isin EM, Guengerich FP: Complex reactions catalyzed by cytochrome P450 enzymes. Biochim Biophys Acta 2007, 1770(3):314–329. 10.1016/j.bbagen.2006.07.003View ArticlePubMedGoogle Scholar
- Nelson D, Werck-Reichhart D: A P450-centric view of plant evolution. Plant J 2011, 66(1):194–211. 10.1111/j.1365-313X.2011.04529.xView ArticlePubMedGoogle Scholar
- Werck-Reichhart D, Feyereisen R: Cytochromes P450: a success story. Genome Biol 2000, 1(6):3003. REVIEWS3003 REVIEWS3003View ArticleGoogle Scholar
- Kushiro T, Okamoto M, Nakabayashi K, Yamagishi K, Kitamura S, Asami T, Hirai N, Koshiba T, Kamiya Y, Nambara E: The Arabidopsis cytochrome P450 CYP707A encodes ABA 8′-hydroxylases: key enzymes in ABA catabolism. EMBO J 2004, 23(7):1647–1656. 10.1038/sj.emboj.7600121PubMed CentralView ArticlePubMedGoogle Scholar
- Pan G, Zhang X, Liu K, Zhang J, Wu X, Zhu J, Tu J: Map-based cloning of a novel rice cytochrome P450 gene CYP81A6 that confers resistance to two different classes of herbicides. Plant Mol Biol 2006, 61(6):933–943. 10.1007/s11103-006-0058-zView ArticlePubMedGoogle Scholar
- Robineau T, Batard Y, Nedelkina S, Cabello-Hurtado F, LeRet M, Sorokine O, Didierjean L, Werck-Reichhart D: The chemically inducible plant cytochrome P450 CYP76B1 actively metabolizes phenylureas and other xenobiotics. Plant Physiol 1998, 118(3):1049–1056. 10.1104/pp.118.3.1049PubMed CentralView ArticlePubMedGoogle Scholar
- Koster J, Thurow C, Kruse K, Meier A, Iven T, Feussner I, Gatz C: Xenobiotic- and Jasmonic Acid-Inducible Signal Transduction Pathways have Become Interdependent at the Arabidopsis thaliana CYP81D11 Promoter. Plant Physiol 2012.Google Scholar
- Fujita S, Ohnishi T, Watanabe B, Yokota T, Takatsuto S, Fujioka S, Yoshida S, Sakata K, Mizutani M: Arabidopsis CYP90B1 catalyses the early C-22 hydroxylation of C27, C28 and C29 sterols. Plant J 2006, 45(5):765–774. 10.1111/j.1365-313X.2005.02639.xView ArticlePubMedGoogle Scholar
- Humphreys JM, Hemm MR, Chapple C: New routes for lignin biosynthesis defined by biochemical characterization of recombinant ferulate 5-hydroxylase, a multifunctional cytochrome P450-dependent monooxygenase. Proc Natl Acad Sci U S A 1999, 96(18):10045–10050. 10.1073/pnas.96.18.10045PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer K, Shirley AM, Cusumano JC, Bell-Lelong DA, Chapple C: Lignin monomer composition is determined by the expression of a cytochrome P450-dependent monooxygenase in Arabidopsis. Proc Natl Acad Sci U S A 1998, 95(12):6619–6623. 10.1073/pnas.95.12.6619PubMed CentralView ArticlePubMedGoogle Scholar
- Morikawa T, Mizutani M, Ohta D: Cytochrome P450 subfamily CYP710A genes encode sterol C-22 desaturase in plants. Biochem Soc Trans 2006, 34(Pt 6):1202–1205.View ArticlePubMedGoogle Scholar
- Schoch G, Goepfert S, Morant M, Hehn A, Meyer D, Ullmann P, Werck-Reichhart D: CYP98A3 from Arabidopsis thaliana is a 3′-hydroxylase of phenolic esters, a missing link in the phenylpropanoid pathway. J Biol Chem 2001, 276(39):36566–36574. 10.1074/jbc.M104047200View ArticlePubMedGoogle Scholar
- Mizutani M, Ohta D: Diversification of P450 genes during land plant evolution. Annu Rev Plant Biol 2010, 61: 291–315. 10.1146/annurev-arplant-042809-112305View ArticlePubMedGoogle Scholar
- Li L, Chang Z, Pan Z, Fu ZQ, Wang X: Modes of heme binding and substrate access for cytochrome P450 CYP74A revealed by crystal structures of allene oxide synthase. Proc Natl Acad Sci U S A 2008, 105(37):13883–13888. 10.1073/pnas.0804099105PubMed CentralView ArticlePubMedGoogle Scholar
- Lee DS, Nioche P, Hamberg M, Raman CS: Structural insights into the evolutionary paths of oxylipin biosynthetic enzymes. Nature 2008, 455(7211):363–368. 10.1038/nature07307View ArticlePubMedGoogle Scholar
- Baudry J, Rupasinghe S, Schuler MA: Class-dependent sequence alignment strategy improves the structural and functional modeling of P450s. Protein Eng Des Sel 2006, 19(8):345–353. 10.1093/protein/gzl012View ArticlePubMedGoogle Scholar
- Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res 2003, 31(13):3381–3385. 10.1093/nar/gkg520PubMed CentralView ArticlePubMedGoogle Scholar
- Wu S, Zhang Y: MUSTER: Improving protein sequence profile-profile alignments by using multiple sources of structure information. Proteins 2008, 72(2):547–556. 10.1002/prot.21945PubMed CentralView ArticlePubMedGoogle Scholar
- Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc 2010, 5(4):725–738. 10.1038/nprot.2010.5PubMed CentralView ArticlePubMedGoogle Scholar
- Xu D, Zhang J, Roy A, Zhang Y: Automated protein structure modeling in CASP9 by I-TASSER pipeline combined with QUARK-based ab initio folding and FG-MD-based structure refinement. Proteins 2011, 79(Suppl 10):147–160.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y: I-TASSER: fully automated protein structure prediction in CASP8. Proteins 2009, 77(Suppl 9):100–113.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y: Template-based modeling and free modeling by I-TASSER in CASP7. Proteins 2007, 69(Suppl 8):108–117.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25(17):3389–3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004, 32(5):1792–1797. 10.1093/nar/gkh340PubMed CentralView ArticlePubMedGoogle Scholar
- Sadreyev RI, Grishin NV: Quality of alignment comparison by COMPASS improves with inclusion of diverse confident homologs. Bioinformatics 2004, 20(6):818–828. 10.1093/bioinformatics/btg485View ArticlePubMedGoogle Scholar
- Sali A, Potterton L, Yuan F, van Vlijmen H, Karplus M: Evaluation of comparative protein modeling by MODELLER. Proteins 1995, 23(3):318–326. 10.1002/prot.340230306View ArticlePubMedGoogle Scholar
- Nelson DR: The cytochrome p450 homepage. Hum Genomics 2009, 4(1):59–65.PubMed CentralPubMedGoogle Scholar
- Sirim D, Wagner F, Lisitsa A, Pleiss J: The cytochrome P450 engineering database: Integration of biochemical properties. BMC Biochem 2009, 10: 27. 10.1186/1471-2091-10-27PubMed CentralView ArticlePubMedGoogle Scholar
- Park J, Lee S, Choi J, Ahn K, Park B, Park J, Kang S, Lee YH: Fungal cytochrome P450 database. BMC Genomics 2008, 9: 402. 10.1186/1471-2164-9-402PubMed CentralView ArticlePubMedGoogle Scholar
- Preissner S, Kroll K, Dunkel M, Senger C, Goldsobel G, Kuzman D, Guenther S, Winnenburg R, Schroeder M, Preissner R: SuperCYP: a comprehensive database on Cytochrome P450 enzymes including a tool for analysis of CYP-drug interactions. Nucleic Acids Res 2010, 38: 237–243. 10.1093/nar/gkp970View ArticleGoogle Scholar
- Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, Burkhardt K, Feng Z, Gilliland GL, Iype L, Jain S, et al.: The Protein Data Bank. Acta Crystallogr D Biol Crystallogr 2002, 58(Pt 6 No 1):899–907.View ArticlePubMedGoogle Scholar
- Paquette SM, Bak S, Feyereisen R: Intron-exon organization and phylogeny in a large superfamily, the paralogous cytochrome P450 genes of Arabidopsis thaliana. DNA Cell Biol 2000, 19(5):307–317. 10.1089/10445490050021221View ArticlePubMedGoogle Scholar
- Rhee SY, Beavis W, Berardini TZ, Chen G, Dixon D, Doyle A, Garcia-Hernandez M, Huala E, Lander G, Montoya M, Miller N, Mueller LA, Mundodi S, Reiser L, Tacklind J, Weems DC, Wu Y, Xu I, Yoo D, Yoon J, Zhang P: The Arabidopsis Information Resource (TAIR): a model organism database providing a centralized, curated gateway to Arabidopsis biology, research materials and community. Nucleic Acids Res 2003, 31(1):224–228. 10.1093/nar/gkg076View ArticlePubMedGoogle Scholar
- Zhang P, Foerster H, Tissier CP, Mueller L, Paley S, Karp PD, Rhee SY: MetaCyc and AraCyc. Metabolic pathway databases for plant research. Plant Physiol 2005, 138(1):27–37.PubMedGoogle Scholar
- Li Q, Cheng T, Wang Y, Bryant SH: PubChem as a public resource for drug discovery. Drug Discov Today 2010, 15(23–24):1052–1057.PubMed CentralView ArticlePubMedGoogle Scholar
- Vicatos S, Reddy BV, Kaznessis Y: Prediction of distant residue contacts with the use of evolutionary information. Proteins 2005, 58(4):935–949. 10.1002/prot.20370View ArticlePubMedGoogle Scholar
- Miyazawa S, Jernigan RL: Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. J Mol Biol 1996, 256(3):623–644. 10.1006/jmbi.1996.0114View ArticlePubMedGoogle Scholar
- Eisenberg D, Luthy R, Bowie JU: VERIFY3D: assessment of protein models with three-dimensional profiles. Methods Enzymol 1997, 277: 396–404.View ArticlePubMedGoogle Scholar
- Dessau RB: Pipper CB: [“R”–project for statistical computing]. Ugeskr Laeger 2008, 170(5):328–330.PubMedGoogle Scholar
- Wu G, Robertson DH, Brooks CL 3rd, Vieth M: Detailed analysis of grid-based molecular docking: A case study of CDOCKER-A CHARMm-based MD docking algorithm. J Comput Chem 2003, 24(13):1549–1562. 10.1002/jcc.10306View ArticlePubMedGoogle Scholar
- Magrane M, Consortium U: UniProt Knowledgebase: a hub of integrated protein data. Database (Oxford) 2011, 2011: 009.View ArticleGoogle Scholar
- Herraez A: Biomolecules in the computer: Jmol to the rescue. Biochem Mol Biol Educ 2006, 34(4):255–261. 10.1002/bmb.2006.494034042644View ArticlePubMedGoogle Scholar
- Sobolev V, Sorokine A, Prilusky J, Abola EE, Edelman M: Automated analysis of interatomic contacts in proteins. Bioinformatics 1999, 15(4):327–332. 10.1093/bioinformatics/15.4.327View ArticlePubMedGoogle Scholar
- Hooft RW, Sander C, Scharf M, Vriend G: The PDBFINDER database: a summary of PDB, DSSP and HSSP information with added value. Comput Appl Biosci 1996, 12(6):525–529.PubMedGoogle Scholar
- Cao Y, Charisi A, Cheng LC, Jiang T, Girke T: ChemmineR: a compound mining framework for R. Bioinformatics 2008, 24(15):1733–1734. 10.1093/bioinformatics/btn307PubMed CentralView ArticlePubMedGoogle Scholar
- Friedman FK, Robinson RC, Dai R: Molecular modeling of mammalian cytochrome P450s. Front Biosci 2004, 9: 2796–2806. 10.2741/1437View ArticlePubMedGoogle Scholar
- Zhang Y, Skolnick J: TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res 2005, 33(7):2302–2309. 10.1093/nar/gki524PubMed CentralView ArticlePubMedGoogle Scholar
- Norgan AP, Coffman PK, Kocher JP, Katzmann DJ, Sosa CP: J Cheminform. 2011, 3(1):12.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.