/* * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$) * Copyright (C) $$Year-Rel$$ The Jalview Authors * * This file is part of Jalview. * * Jalview is free software: you can redistribute it and/or * modify it under the terms of the GNU General Public License * as published by the Free Software Foundation, either version 3 * of the License, or (at your option) any later version. * * Jalview is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty * of MERCHANTABILITY or FITNESS FOR A PARTICULAR * PURPOSE. See the GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with Jalview. If not, see . * The Jalview Authors are detailed in the 'AUTHORS' file. */ package jalview.ws.dbsources; import static org.testng.AssertJUnit.assertEquals; import static org.testng.AssertJUnit.assertFalse; import static org.testng.AssertJUnit.assertNotNull; import static org.testng.AssertJUnit.assertNull; import jalview.datamodel.PDBEntry; import jalview.datamodel.SequenceI; import jalview.datamodel.xdb.uniprot.UniprotEntry; import jalview.datamodel.xdb.uniprot.UniprotFeature; import jalview.gui.JvOptionPane; import java.io.Reader; import java.io.StringReader; import java.util.Vector; import org.testng.Assert; import org.testng.annotations.BeforeClass; import org.testng.annotations.Test; public class UniprotTest { @BeforeClass(alwaysRun = true) public void setUpJvOptionPane() { JvOptionPane.setInteractiveMode(false); JvOptionPane.setMockResponse(JvOptionPane.CANCEL_OPTION); } // adapted from http://www.uniprot.org/uniprot/A9CKP4.xml private static final String UNIPROT_XML = "" + "" + "" + "A9CKP4" + "A9CKP5" + "A9CKP4_AGRT5" + "A9CKP4_AGRT6" + "Mitogen-activated protein kinase 13Henry" + "" + "" + "" + "" + "" + "" + "ML" + "ML" + "M" + "LLMVM" + "LLLMVML" + "LLLMVKMLML" + "MHAPL VSKDL" + ""; /** * Test the method that unmarshals XML to a Uniprot model */ @Test(groups = { "Functional" }) public void testGetUniprotEntries() { Uniprot u = new Uniprot(); Reader reader = new StringReader(UNIPROT_XML); Vector entries = u.getUniprotEntries(reader); assertEquals(1, entries.size()); UniprotEntry entry = entries.get(0); assertEquals(2, entry.getName().size()); assertEquals("A9CKP4_AGRT5", entry.getName().get(0)); assertEquals("A9CKP4_AGRT6", entry.getName().get(1)); assertEquals(2, entry.getAccession().size()); assertEquals("A9CKP4", entry.getAccession().get(0)); assertEquals("A9CKP5", entry.getAccession().get(1)); /* * UniprotSequence drops any space characters */ assertEquals("MHAPLVSKDL", entry.getUniprotSequence().getContent()); assertEquals(2, entry.getProtein().getName().size()); assertEquals("Mitogen-activated protein kinase 13", entry.getProtein() .getName().get(0)); assertEquals("Henry", entry.getProtein().getName().get(1)); /* * Check sequence features */ Vector features = entry.getFeature(); assertEquals(9, features.size()); UniprotFeature sf = features.get(0); assertEquals("signal peptide", sf.getType()); assertNull(sf.getDescription()); assertNull(sf.getStatus()); assertEquals(1, sf.getBegin()); assertEquals(18, sf.getEnd()); sf = features.get(1); assertEquals("propeptide", sf.getType()); assertEquals("Activation peptide", sf.getDescription()); assertEquals(19, sf.getPosition()); assertEquals(19, sf.getBegin()); assertEquals(20, sf.getEnd()); sf = features.get(2); assertEquals("chain", sf.getType()); assertEquals("Granzyme B", sf.getDescription()); assertEquals(21, sf.getPosition()); assertEquals(21, sf.getBegin()); assertEquals(247, sf.getEnd()); sf = features.get(3); assertEquals("sequence variant", sf.getType()); assertNull(sf.getDescription()); assertEquals(41, sf.getPosition()); assertEquals(41, sf.getBegin()); assertEquals(41, sf.getEnd()); sf = features.get(4); assertEquals("sequence variant", sf.getType()); assertEquals("Pathogenic", sf.getDescription()); assertEquals(41, sf.getPosition()); assertEquals(41, sf.getBegin()); assertEquals(41, sf.getEnd()); sf = features.get(5); assertEquals("sequence variant", sf.getType()); assertEquals("Pathogenic", sf.getDescription()); assertEquals(41, sf.getPosition()); assertEquals(41, sf.getBegin()); assertEquals(41, sf.getEnd()); sf = features.get(6); assertEquals("sequence variant", sf.getType()); assertEquals("Foo", sf.getDescription()); assertEquals(42, sf.getPosition()); assertEquals(42, sf.getBegin()); assertEquals(42, sf.getEnd()); Assert.assertEquals(Uniprot.getDescription(sf), "p.Met42Leu" + "\n" + "p.Met42LeuMetVal Foo"); sf = features.get(7); assertEquals(42, sf.getBegin()); assertEquals(43, sf.getEnd()); Assert.assertEquals(Uniprot.getDescription(sf), "p.MetLeu42LeuLeu" + "\n" + "p.MetLeu42LeuMetVal Foo"); sf = features.get(8); assertEquals(42, sf.getBegin()); assertEquals(45, sf.getEnd()); Assert.assertEquals(Uniprot.getDescription(sf), "p.MLML42LeuLeu" + "\n" + "p.MLML42LMVK Foo Too"); /* * Check cross-references */ Vector xrefs = entry.getDbReference(); assertEquals(3, xrefs.size()); PDBEntry xref = xrefs.get(0); assertEquals("2FSQ", xref.getId()); assertEquals("PDB", xref.getType()); assertEquals("X-ray", xref.getProperty("method")); assertEquals("1.40", xref.getProperty("resolution")); xref = xrefs.get(1); assertEquals("2FSR", xref.getId()); assertEquals("PDBsum", xref.getType()); assertFalse(xref.getProperties().hasMoreElements()); xref = xrefs.get(2); assertEquals("AE007869", xref.getId()); assertEquals("EMBL", xref.getType()); assertEquals("AAK85932.1", xref.getProperty("protein sequence ID")); assertEquals("Genomic_DNA", xref.getProperty("molecule type")); } @Test(groups = { "Functional" }) public void testGetUniprotSequence() { UniprotEntry entry = new Uniprot().getUniprotEntries( new StringReader(UNIPROT_XML)).get(0); SequenceI seq = new Uniprot().uniprotEntryToSequenceI(entry); assertNotNull(seq); assertEquals(6, seq.getDBRefs().length); // 2*Uniprot, PDB, PDBsum, 2*EMBL } /** * Test the method that formats the sequence id */ @Test(groups = { "Functional" }) public void testGetUniprotEntryId() { UniprotEntry entry = new Uniprot().getUniprotEntries( new StringReader(UNIPROT_XML)).get(0); /* * name formatted with Uniprot Entry name */ String expectedName = "A9CKP4_AGRT5|A9CKP4_AGRT6"; assertEquals(expectedName, Uniprot.getUniprotEntryId(entry)); } /** * Test the method that formats the sequence description */ @Test(groups = { "Functional" }) public void testGetUniprotEntryDescription() { UniprotEntry entry = new Uniprot().getUniprotEntries( new StringReader(UNIPROT_XML)).get(0); /* * recommended names concatenated with space separator */ String expectedDescription = "Mitogen-activated protein kinase 13 Henry"; assertEquals(expectedDescription, Uniprot.getUniprotEntryDescription(entry)); } @Test(groups = { "Functional" }) public void testGetDescription() { UniprotFeature uf = new UniprotFeature(); assertEquals("", Uniprot.getDescription(uf)); uf.setDescription("Hello"); assertEquals("Hello", Uniprot.getDescription(uf)); uf.setPosition(23); uf.setOriginal("K"); Vector vars = new Vector(); vars.add("y"); uf.setVariation(vars); assertEquals("p.Lys23Tyr Hello", Uniprot.getDescription(uf)); vars.clear(); vars.add("z"); // unknown variant - fails gracefully uf.setVariation(vars); assertEquals("p.Lys23z Hello", Uniprot.getDescription(uf)); uf.setVariation(null); // variant missing - is ignored assertEquals("Hello", Uniprot.getDescription(uf)); } }