[Biojava-l] Problem finding Restriction Enzyme sites

[Richard Holland]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hello.

Your code is throwing an exception saying that it doesn't know what the
definition for ECOR1 is. This is because the molbio package doesn't
contain any restriction enzyme definitions by default, and you haven't
loaded any into it yet (at least, not in the code fragment you quote).

Before you can use the methods in this package, you need to have made a
call to RestrictionEnzymeManager.loadEnzymeFile() and pass it a
reference to a file containing some restriction enzyme definitions in
REBASE format #31.

Once you've used this method to load a REBASE file containing a
definition for the ECOR1 enzyme, you'll be able to use the code.

cheers,
Richard

Parsa Hosseini wrote:
> Hi,
> 
> I've been using the 'molbio' package to do some Restriction enzyme analysis. I've come across errors when I want to find the corresponding recognition site. I'd be very thankful if I could be lead down the path to find were I am going wrong.
> 
> Thank you in advance,
> 
> Parsa
> 
> phossein at umd.edu
> 
> ------------------------
> 
> try {
> 			
> SymbolList siteOfEnzyme = DNATools.createDNA("GAATC");
> 			Sequence sequence = DNATools.createDNASequence("AAAAGAATCTTC", "mySequence");
> 			RestrictionEnzyme re1 = new RestrictionEnzyme("ECOR1", siteOfEnzyme, 0, 0);
> SimpleThreadPool tr = new SimpleThreadPool();
> RestrictionMapper reMapper = new RestrictionMapper(tr);
> reMapper.addEnzyme(re1);			System.out.println(reMapper.annotate(sequence));		
> 		}
> 		catch (IllegalSymbolException e) {
> 			
> 		}
> 		catch (IllegalAlphabetException e) {
> 			
> 		}
> 
> Exception in thread "Thread-1" org.biojava.bio.BioRuntimeException: Failed to complete search for ECOR1 GAATC (0/0) at org.biojava.bio.molbio.RestrictionSiteFinder.run(RestrictionSiteFinder.java:136) at org.biojava.utils.SimpleThreadPool$PooledThread.run(SimpleThreadPool.java:295)
> Caused by: java.lang.IllegalArgumentException: RestrictionEnzyme 'ECOR1' is not registered. No precompiled Pattern is available
> 	at org.biojava.bio.molbio.RestrictionEnzymeManager.getPatterns(RestrictionEnzymeManager.java:280)
> 	at org.biojava.bio.molbio.RestrictionSiteFinder.run(RestrictionSiteFinder.java:77)
> 	... 1 more
> 
> 
> 
> 
> 
> _______________________________________________
> Biojava-l mailing list  -  Biojava-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/biojava-l
> 

- --
Richard Holland (BioMart)
EMBL EBI, Wellcome Trust Genome Campus,
Hinxton, Cambridgeshire CB10 1SD, UK
Tel. +44 (0)1223 494416

http://www.biomart.org/
http://www.biojava.org/
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.2.2 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org

iD8DBQFHZ4tD4C5LeMEKA/QRApOrAKCYZ9rXbPzyKuAV4DrwaSL424s/wACfSIDC
K6TFXwVdQSrV4bTG4DGKYhI=
=Nc4w
-----END PGP SIGNATURE-----