/**
* KEGG get names by cpd id.
*
* @param identifier the identifier
* @return the string
*/
public static String[] KEGGgetNameByCpd(String identifier) {
String ret = "";
try {
KEGGLocator locator = new KEGGLocator();
KEGGPortType serv = locator.getKEGGPort();
// get molecule by accession ID
String str = "cpd:" + identifier;
ret = serv.bget(str);
} catch (Exception e) {
System.err.println(e.toString());
}
// System.out.println(ret);
String names = ret.substring(ret.indexOf("NAME") + 12, ret.indexOf("FORMULA") - 1);
// System.out.println(names);
// find white spaces
Pattern p = Pattern.compile("\\s+");
Matcher m = p.matcher(names);
String cleanedNames = m.replaceAll(" ");
// System.out.println(cleanedNames);
String[] namesArr = cleanedNames.split(";");
for (int i = 0; i < namesArr.length; i++) {
Pattern p1 = Pattern.compile("^\\s+");
Matcher m1 = p1.matcher(namesArr[i]);
namesArr[i] = m1.replaceAll("");
}
return namesArr;
}
/**
* Gets the mol file from a specified KEGG ID.
*
* @param accessionID the accession id (e.g. C00509)
* @return the string
* @throws CDKException
*/
public static IAtomContainer getMol(String accessionID, String keggPath, boolean getAll)
throws CDKException {
String ret = "";
try {
KEGGLocator locator = new KEGGLocator();
KEGGPortType serv = locator.getKEGGPort();
// get molecule by accession ID
String str = "-f m cpd:" + accessionID;
boolean webservice = false;
String mol = "";
try {
// Open the file that is the first
// command line parameter
FileInputStream fstream = new FileInputStream(keggPath + accessionID + ".mol");
// Get the object of DataInputStream
DataInputStream in = new DataInputStream(fstream);
BufferedReader br = new BufferedReader(new InputStreamReader(in));
// Read File Line By Line
String line;
while ((line = br.readLine()) != null) {
// Print the content on the console
mol += line + "\n";
}
// Close the input stream
in.close();
} catch (FileNotFoundException e) {
webservice = true;
}
// if this file is not found use webservice...otherwise use local molfile
if (webservice) ret = serv.bget(str);
else ret = mol;
} catch (Exception e) {
System.err.println(e.toString());
}
MDLReader reader;
List<IAtomContainer> containersList;
reader = new MDLReader(new StringReader(ret));
ChemFile chemFile = (ChemFile) reader.read((ChemObject) new ChemFile());
containersList = ChemFileManipulator.getAllAtomContainers(chemFile);
IAtomContainer molecule = containersList.get(0);
if (getAll) return molecule;
if (!MolecularFormulaTools.isBiologicalCompound(molecule)) molecule = null;
return molecule;
}
/**
* Gets the mol file from a specified KEGG ID.
*
* @param accessionID the accession id (e.g. C00509)
* @return the string
*/
public static String KEGGgetMol(String accessionID, String keggPath) {
String ret = "";
try {
KEGGLocator locator = new KEGGLocator();
KEGGPortType serv = locator.getKEGGPort();
// get molecule by accession ID
String str = "-f m cpd:" + accessionID;
boolean webservice = false;
String mol = "";
try {
// Open the file that is the first
// command line parameter
FileInputStream fstream = new FileInputStream(keggPath + accessionID + ".mol");
// Get the object of DataInputStream
DataInputStream in = new DataInputStream(fstream);
BufferedReader br = new BufferedReader(new InputStreamReader(in));
// Read File Line By Line
String line;
while ((line = br.readLine()) != null) {
// Print the content on the console
mol += line + "\n";
}
// Close the input stream
in.close();
} catch (FileNotFoundException e) {
webservice = true;
}
// if this file is not found use webservice...otherwise use local molfile
if (webservice) ret = serv.bget(str);
else ret = mol;
} catch (Exception e) {
System.err.println(e.toString());
}
return ret;
}