desc = GlobalDescriptor([database['Sequence'][i]])
desc.calculate_charge(ph=7, amide=True)
charge_array.append(desc.descriptor[0][0])
except:
charge_array.append('')
database['charge'] = charge_array
print("Estimate charge_density")
#calculate charge density for each sequence
charge_density_array = []
for i in range(len(database)):
try:
desc = GlobalDescriptor([database['Sequence'][i]])
desc.charge_density(ph=7, amide=True)
charge_density_array.append(desc.descriptor[0][0])
except:
charge_density_array.append('')
database['charge_density'] = charge_density_array
print("Estimate isoelectric")
#estimate isoelectric point
isoelectric_array = []
for i in range(len(database)):
try:
desc = GlobalDescriptor([database['Sequence'][i]])
desc.isoelectric_point()
isoelectric_array.append(desc.descriptor[0][0])
gravy_long = []
for gp in peptides:
#
eisenbergs[gp] = get_peptide_values(peptides[gp], 'eisenberg')
for val in eisenbergs[gp]:
eisenbergs_long.append([gp, val])
#
properties = GlobalDescriptor(peptides[gp])
properties.calculate_charge(ph=7.4, amide=True)
charges[gp] = [x[0] for x in properties.descriptor]
for val in charges[gp]:
charges_long.append([gp, val])
#
properties = GlobalDescriptor(peptides[gp])
properties.charge_density(ph=7.4, amide=True)
charge_densities[gp] = [x[0] for x in properties.descriptor]
for val in charge_densities[gp]:
charge_densities_long.append([gp, val])
#
polarities[gp] = get_peptide_values(peptides[gp], 'polarity')
for val in polarities[gp]:
polarities_long.append([gp, val])
#
gravy[gp] = get_peptide_values(peptides[gp], 'gravy')
for val in gravy[gp]:
gravy_long.append([gp, val])
eisenbergs_long = fxn.list_to_df(eisenbergs_long, ['Group', 'Hydrophobicity'],
False)
charges_long = fxn.list_to_df(charges_long, ['Group', 'Charge'], False)
def exec(peptide, time_node):
file = open("../src/public/jobs/service1/service1.fasta", "w")
file.write(peptide)
file.close()
fasta = SeqIO.parse("../src/public/jobs/service1/service1.fasta", "fasta")
if(any(fasta) == False): #False when `fasta` is empty
return "error"
cantidad = 0
for record in SeqIO.parse("../src/public/jobs/service1/service1.fasta", "fasta"):
cantidad = cantidad+1
if (cantidad == 1):
properties = {}
for record in SeqIO.parse("../src/public/jobs/service1/service1.fasta", "fasta"):
properties[str(record.id)] = {}
#save properties
properties[str(record.id)]["length"] = len(record.seq)
#formula
try:
desc = GlobalDescriptor(str(record.seq))
desc.formula(amide=True)
properties[str(record.id)]["formula"] = desc.descriptor[0][0]
except:
properties[str(record.id)]["formula"] = "-"
#molecular weigth
try:
desc = GlobalDescriptor(str(record.seq))
desc.calculate_MW(amide=True)
properties[str(record.id)]["molecular_weigth"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["molecular_weigth"] = "-"
#boman_index
try:
desc = GlobalDescriptor(str(record.seq))
desc.boman_index()
properties[str(record.id)]["boman_index"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["boman_index"] = "-"
#charge
try:
desc = GlobalDescriptor(str(record.seq))
desc.calculate_charge(ph=7, amide=True)
properties[str(record.id)]["charge"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["charge"] = "-"
#charge density
try:
desc = GlobalDescriptor(str(record.seq))
desc.charge_density(ph=7, amide=True)
properties[str(record.id)]["charge_density"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["charge_density"] = "-"
#estimate isoelectric point
try:
desc = GlobalDescriptor(str(record.seq))
desc.isoelectric_point()
properties[str(record.id)]["isoelectric_point"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["isoelectric_point"] = "-"
#estimate inestability index
try:
desc = GlobalDescriptor(str(record.seq))
desc.instability_index()
properties[str(record.id)]["instability_index"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["instability_index"] = "-"
#estimate aromaticity
try:
desc = GlobalDescriptor(str(record.seq))
desc.aromaticity()
properties[str(record.id)]["aromaticity"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["aromaticity"] = "-"
#estimate aliphatic_index
try:
desc = GlobalDescriptor(str(record.seq))
desc.aliphatic_index()
properties[str(record.id)]["aliphatic_index"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["aliphatic_index"] = "-"
#estimate hydrophobic_ratio
try:
desc = GlobalDescriptor(str(record.seq))
desc.hydrophobic_ratio()
properties[str(record.id)]["hydrophobic_ratio"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["hydrophobic_ratio"] = "-"
#profile hydrophobicity
try:
desc = PeptideDescriptor(str(record.seq), scalename='Eisenberg')
desc.calculate_profile(prof_type='H')
properties[str(record.id)]["hydrophobicity_profile"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["hydrophobicity_profile"] = "-"
#profile hydrophobic
try:
desc = PeptideDescriptor(str(record.seq), scalename='Eisenberg')
desc.calculate_profile(prof_type='uH')
properties[str(record.id)]["hydrophobic_profile"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["hydrophobic_profile"] = "-"
#moment
try:
desc = PeptideDescriptor(str(record.seq), scalename='Eisenberg')
desc.calculate_moment()
properties[str(record.id)]["calculate_moment"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["calculate_moment"] = "-"
try:
os.mkdir("../src/public/jobs/service1/"+time_node)
except:
print("Error")
#generate plot profile
plot_profile(str(record.seq), scalename='eisenberg', filename= "../src/public/jobs/service1/"+time_node+"/profile.png")
#generate helical wheel
helical_wheel(str(record.seq), colorcoding='charge', lineweights=False, filename= "../src/public/jobs/service1/"+time_node+"/helical.png")
return(properties)
if (cantidad > 1):
properties = {}
for record in SeqIO.parse("../src/public/jobs/service1/service1.fasta", "fasta"):
properties[str(record.id)] = {}
properties[str(record.id)]["length"] = len(record.seq)
#formula
try:
desc = GlobalDescriptor(str(record.seq))
desc.formula(amide=True)
properties[str(record.id)]["formula"] = desc.descriptor[0][0]
except:
properties[str(record.id)]["formula"] = "-"
#molecular weigth
try:
desc = GlobalDescriptor(str(record.seq))
desc.calculate_MW(amide=True)
properties[str(record.id)]["molecular_weigth"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["molecular_weigth"] = "-"
#boman_index
try:
desc = GlobalDescriptor(str(record.seq))
desc.boman_index()
properties[str(record.id)]["boman_index"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["boman_index"] = "-"
#charge
try:
desc = GlobalDescriptor(str(record.seq))
desc.calculate_charge(ph=7, amide=True)
properties[str(record.id)]["charge"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["charge"] = "-"
#charge density
try:
desc = GlobalDescriptor(str(record.seq))
desc.charge_density(ph=7, amide=True)
properties[str(record.id)]["charge_density"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["charge_density"] = "-"
#estimate isoelectric point
try:
desc = GlobalDescriptor(str(record.seq))
desc.isoelectric_point()
properties[str(record.id)]["isoelectric_point"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["isoelectric_point"] = "-"
#estimate inestability index
try:
desc = GlobalDescriptor(str(record.seq))
desc.instability_index()
properties[str(record.id)]["instability_index"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["instability_index"] = "-"
#estimate aromaticity
try:
desc = GlobalDescriptor(str(record.seq))
desc.aromaticity()
properties[str(record.id)]["aromaticity"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["aromaticity"] = "-"
#estimate aliphatic_index
try:
desc = GlobalDescriptor(str(record.seq))
desc.aliphatic_index()
properties[str(record.id)]["aliphatic_index"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["aliphatic_index"] = "-"
#estimate hydrophobic_ratio
try:
desc = GlobalDescriptor(str(record.seq))
desc.hydrophobic_ratio()
properties[str(record.id)]["hydrophobic_ratio"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["hydrophobic_ratio"] = "-"
#profile hydrophobicity
try:
desc = PeptideDescriptor(str(record.seq), scalename='Eisenberg')
desc.calculate_profile(prof_type='H')
properties[str(record.id)]["hydrophobicity_profile"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["hydrophobicity_profile"] = "-"
#profile hydrophobic
try:
desc = PeptideDescriptor(str(record.seq), scalename='Eisenberg')
desc.calculate_profile(prof_type='uH')
properties[str(record.id)]["hydrophobic_profile"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["hydrophobic_profile"] = "-"
#moment
try:
desc = PeptideDescriptor(str(record.seq), scalename='Eisenberg')
desc.calculate_moment()
properties[str(record.id)]["calculate_moment"] = float("%.4f" % desc.descriptor[0][0])
except:
properties[str(record.id)]["calculate_moment"] = "-"
return(properties)