def kegg_rest(type: str,
argument: str,
optional_argument: str = "",
sleep_time: float = .5) -> List[str]:
"""This function calls Biopython's KEGG REST function and returns the lines as a string list.
All empty lines are deleted from the list as they do not contain any information.
Arguments
----------
* type: str ~ The KEGG REST action. Can be either 'info', 'get', 'link' or 'list.
* argument: str ~ The argument for the KEGG order.
* optional_argument: str="" ~ The second argument which is necessary for 'link' and 'list' actions
to work correctly.
* sleep_time: float=10.0 ~ The time that shall be waited after a REST action is performed.
Its default value of 10.0 seconds is in accordance with the NCBI
rule that its servers shall not be contacted more often than every
10 seconds. KEGG might have lower or higher required sleep times,
but I did not find any specified time step.
"""
# Execute correct Biotpython KEGG REST function.
if type == "info":
kegg_data = REST.kegg_info(argument)
elif type == "get":
kegg_data = REST.kegg_get(argument)
elif type == "link":
kegg_data = REST.kegg_link(argument, optional_argument)
elif type == "list":
kegg_data = REST.kegg_list(argument, optional_argument)
elif type == "find":
kegg_data = REST.kegg_find(argument, optional_argument)
# Wait the sleep time doing nothing.
time.sleep(sleep_time)
# Get one string per line of the KEGG REST result.
lines: List[str] = kegg_data.read().split("\n")
# Delete empty lines.
not_empty_lines: List[str] = [i for i in lines if len(i) > 0]
return not_empty_lines
def __retrieve_info(self, dbs):
current = dict()
## Get info from http://rest.kegg.jp/info/kegg
## Note that this does not necessarily match what is returned from lists
raw_list = REST.kegg_info("kegg")
split = raw_list.read().splitlines()
current["count_info"] = dict()
for line in split:
split_line = line.split()
overlap = set(dbs) & set(split_line)
if "Release" in split_line:
current["release_short"] = float(
split_line[2].split("/")[0][:-1])
release_long = [split_line[2][:-1]] + split_line[3:]
current["release_long"] = ("_").join(release_long)
elif overlap:
current["count_info"][split_line[0]] = int(
split_line[1].replace(',', ''))
return current
def t_KEGG_Query():
"""Tests Bio.KEGG API Wrapper"""
print("Testing Bio.KEGG.query\n\n")
# info tests
resp = REST.kegg_info("kegg")
resp.read()
print(resp.url)
resp = REST.kegg_info("pathway")
resp.read()
print(resp.url)
# list tests
resp = REST.kegg_list("pathway")
resp.read()
print(resp.url)
resp = REST.kegg_list("pathway", "hsa")
resp.read()
print(resp.url)
resp = REST.kegg_list("organism")
resp.read()
print(resp.url)
resp = REST.kegg_list("hsa")
resp.read()
print(resp.url)
resp = REST.kegg_list("T01001")
resp.read()
print(resp.url)
resp = REST.kegg_list("hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_list(["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
resp = REST.kegg_list("cpd:C01290+gl:G00092")
resp.read()
print(resp.url)
resp = REST.kegg_list(["cpd:C01290", "gl:G00092"])
resp.read()
print(resp.url)
resp = REST.kegg_list("C01290+G00092")
resp.read()
print(resp.url)
resp = REST.kegg_list(["C01290", "G00092"])
resp.read()
print(resp.url)
# find tests
resp = REST.kegg_find("genes", "shiga+toxin")
resp.read()
print(resp.url)
resp = REST.kegg_find("genes", ["shiga", "toxin"])
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "C7H10O5", "formula")
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "O5C7", "formula")
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "174.05", "exact_mass")
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "300-310", "mol_weight")
resp.read()
print(resp.url)
# get tests
resp = REST.kegg_get("cpd:C01290+gl:G00092")
resp.read()
print(resp.url)
resp = REST.kegg_get(["cpd:C01290", "gl:G00092"])
resp.read()
print(resp.url)
resp = REST.kegg_get("C01290+G00092")
resp.read()
print(resp.url)
resp = REST.kegg_get(["C01290", "G00092"])
resp.read()
print(resp.url)
resp = REST.kegg_get("hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_get(["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
resp = REST.kegg_get("hsa:10458+ece:Z5100", "aaseq")
resp.read()
print(resp.url)
resp = REST.kegg_get(["hsa:10458", "ece:Z5100"], "aaseq")
resp.read()
print(resp.url)
resp = REST.kegg_get("hsa05130", "image")
resp.read()
print(resp.url)
# conv tests
resp = REST.kegg_conv("eco", "ncbi-geneid")
resp.read()
print(resp.url)
resp = REST.kegg_conv("ncbi-geneid", "eco")
resp.read()
print(resp.url)
resp = REST.kegg_conv("ncbi-gi", "hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_conv("ncbi-gi", ["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
# link tests
resp = REST.kegg_link("pathway", "hsa")
resp.read()
print(resp.url)
resp = REST.kegg_link("hsa", "pathway")
resp.read()
print(resp.url)
resp = REST.kegg_link("pathway", "hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_link("pathway", ["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
def t_KEGG_Query():
"""Tests Bio.KEGG API Wrapper"""
print("Testing Bio.KEGG.query\n\n")
# info tests
resp = REST.kegg_info("kegg")
resp.read()
print(resp.url)
resp = REST.kegg_info("pathway")
resp.read()
print(resp.url)
# list tests
resp = REST.kegg_list("pathway")
resp.read()
print(resp.url)
resp = REST.kegg_list("pathway", "hsa")
resp.read()
print(resp.url)
resp = REST.kegg_list("organism")
resp.read()
print(resp.url)
resp = REST.kegg_list("hsa")
resp.read()
print(resp.url)
resp = REST.kegg_list("T01001")
resp.read()
print(resp.url)
resp = REST.kegg_list("hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_list(["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
resp = REST.kegg_list("cpd:C01290+gl:G00092")
resp.read()
print(resp.url)
resp = REST.kegg_list(["cpd:C01290", "gl:G00092"])
resp.read()
print(resp.url)
resp = REST.kegg_list("C01290+G00092")
resp.read()
print(resp.url)
resp = REST.kegg_list(["C01290", "G00092"])
resp.read()
print(resp.url)
# find tests
resp = REST.kegg_find("genes", "shiga+toxin")
resp.read()
print(resp.url)
resp = REST.kegg_find("genes", ["shiga", "toxin"])
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "C7H10O5", "formula")
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "O5C7", "formula")
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "174.05", "exact_mass")
resp.read()
print(resp.url)
resp = REST.kegg_find("compound", "300-310", "mol_weight")
resp.read()
print(resp.url)
# get tests
resp = REST.kegg_get("cpd:C01290+gl:G00092")
resp.read()
print(resp.url)
resp = REST.kegg_get(["cpd:C01290", "gl:G00092"])
resp.read()
print(resp.url)
resp = REST.kegg_get("C01290+G00092")
resp.read()
print(resp.url)
resp = REST.kegg_get(["C01290", "G00092"])
resp.read()
print(resp.url)
resp = REST.kegg_get("hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_get(["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
resp = REST.kegg_get("hsa:10458+ece:Z5100", "aaseq")
resp.read()
print(resp.url)
resp = REST.kegg_get(["hsa:10458", "ece:Z5100"], "aaseq")
resp.read()
print(resp.url)
resp = REST.kegg_get("hsa05130", "image")
resp.read()
print(resp.url)
# conv tests
resp = REST.kegg_conv("eco", "ncbi-geneid")
resp.read()
print(resp.url)
resp = REST.kegg_conv("ncbi-geneid", "eco")
resp.read()
print(resp.url)
resp = REST.kegg_conv("ncbi-gi", "hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_conv("ncbi-gi", ["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
# link tests
resp = REST.kegg_link("pathway", "hsa")
resp.read()
print(resp.url)
resp = REST.kegg_link("hsa", "pathway")
resp.read()
print(resp.url)
resp = REST.kegg_link("pathway", "hsa:10458+ece:Z5100")
resp.read()
print(resp.url)
resp = REST.kegg_link("pathway", ["hsa:10458", "ece:Z5100"])
resp.read()
print(resp.url)
#print(mydog4)
mydog5 = REST.kegg_link("genes", mydog4).read()
return (mydog5)
def unique(list1):
x = np.array(list1)
return (np.unique(x))
def specie_kegg(dddd):
mydog6 = re.findall('[a-z]+', dddd)
return (mydog6)
a = REST.kegg_info('genes').read()
#print(a)
#print(kooo("YBR248C"))
#print(unique(specie_kegg(kooo("YBR248C"))))
def specie_kegg1(ddd):
return (specie_kegg(kooo(ddd)))
specie = pd.read_pickle('specie')
def two_words(input):
return (re.findall('^\S+[ ]*\S*', str(input))[0])
def KEGG(input1, input2):
# Perform the query
result = REST.kegg_info("kegg").read()
# Print overview
if input1 == "info" and input2 == "alt":
return print(result)
# Get all entries in the PATHWAY database as a dataframe
elif input1 == "pathway_overview" and input2 == "alle":
result = REST.kegg_list("pathway").read()
return to_df(result)
#Print alle biosynteseveje
elif input1 == "print_pathway":
if input2 == "alle":
result = REST.kegg_get("map01100", "image").read()
img = Image(result, width=1500, height=1000)
else:
result = REST.kegg_get(input2, "image").read()
img = Image(result)
return img
#Find the compund vanillin
elif input1 == "find_molekyle" and input2 != None:
result = REST.kegg_find("compound", input2).read() #cpd:C00755
return print(result)
elif input1 == "info_molekyle" and input2 != None: #cpd:C00755
# Get the entry information for vanillin
result = REST.kegg_get(input2).read()
return print(result)
# Display molekylær struktur for cpd:C00051 (vanillin)
elif input1 == "molekyle billede" and input2 != None:
result = REST.kegg_get(input2, "image").read() #"cpd:C00755"
return Image(result)
elif input1 == "Enzyme molekyle" and input2 != None:
result = REST.kegg_find("enzyme", input2).read()
return to_df(result)
#Enzym delen
from Bio.KEGG import Enzyme
request = REST.kegg_get(input1)
records = Enzyme.parse(request)
record = list(records)[0]
if input2 == "reaction":
return record.reaction
elif input2 == "classname":
return record.classname
elif input2 == "genes":
genes = list()
for g in record.genes:
gene_id, gene_symbol = g
genes.append(gene_id)
return genes
else:
print("Du har indskrevet nogget der ikke er gældende. Prøv igen")
