Question

Use python 2.7 to write this code: This code takes in a bunch of protein sequences as input listed below and tracks the highest average TM value (number under each protein) and If a protein has the highest TM value of all proteins found so far, track it. At the very end, the program should print out that as well:

For example, at the very end when run, it should print something like this:

Protein with the highest TM value: Protein 18: MFFFFKK...|

Average TM value:   0.6

test code to run the code:

AverageTMValues=

protein 1:'MKLVVRPWAGCWWSTLGPRGSLSPLGICPLLMLLWATLR'

aveTMValue:0.6

protein2:'MARKCSVPLVMAWLTWTTSRAPLPH',

aveTMValue: 0.5

protein3:'MPWPTSITXXXXXXSWSPEWLSSGLRSILGWEQPRVSHKGHSHEWHRRP'

aveTMValue:0.1

protein4: 'MAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG'

aveTMValue:0.5

HighestAverageTMvalue=highestAveTMValue(AverageTMValues)

The code should output/print out: (it should also print out the protein corresponding with the highest TM value)

protein 1:'MKLVVRPWAGCWWSTLGPRGSLSPLGICPLLMLLWATLR'

aveTMValue:0.6

clearly show work and run it. THANK YOU a ton!!

Answer 1

Screenshot

#Function to find highest average TM value from the passed list e C:\WINDOWS\system32\cmd.exe def highestAve TMValue (Average TMValues): protein 1: 'MKLVVRPWAGCWWSTLGPRGSLSPLGICPLLMLLWATLR' #index start aveTMValue:0.6 i=1 Press any key to continue ... #result index index=i #set highest TM avg=float(Average TMValues[i].split(':')[1]) #Loopuntil end of the list while(i<len (Average TMValues )): #Compare average with eac protein averages in the list #Accordingly change index and average if(avg<float(AverageTMValues[i].split(':')[1])): avg=float(AverageTMValues[i].split(':')[1]) index=i i=i+2 #Return protein details return Average TMValues[index-1]+"\n"+AverageTMValues[index] #Test list Average TMValues=["protein 1: 'MKLWRPWAGCWWSTLGPRGSLSPLGICPLLMLLWA "protein2: MARKCSVPLVMAWLTWTTSRAPLPH'","aveTMValue: 0.5", "protein3: 'MPWPTSITXXXXXXSWSPEWLSSGLRSILGWEQPRVSHKGHSHEWHRRP'"," "protein4: 'MAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG'","a #Call method HighestAverageTMvalue=highestAve TMValue (AverageTMValues)

Program

#Function to find highest average TM value from the passed list
def highestAveTMValue(AverageTMValues):
    #index start
    i=1
    #result index
    index=i
    #set highest TM
    avg=float(AverageTMValues[i].split(':')[1])
    #Loopuntil end of the list
    while(i<len(AverageTMValues)):
        #Compare average with eac protein averages in the list
        #Accordingly change index and average
        if(avg<float(AverageTMValues[i].split(':')[1])):
            avg=float(AverageTMValues[i].split(':')[1])
            index=i
        i=i+2
    #Return protein details
    return AverageTMValues[index-1]+"\n"+AverageTMValues[index]
#Test list
AverageTMValues=["protein 1:'MKLVVRPWAGCWWSTLGPRGSLSPLGICPLLMLLWATLR'","aveTMValue:0.6",
"protein2:'MARKCSVPLVMAWLTWTTSRAPLPH'","aveTMValue:0.5",
"protein3:'MPWPTSITXXXXXXSWSPEWLSSGLRSILGWEQPRVSHKGHSHEWHRRP'","aveTMValue:0.1",
"protein4:'MAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG'","aveTMValue:0.5"]
#Call method
HighestAverageTMvalue=highestAveTMValue(AverageTMValues)
#Display result
print(HighestAverageTMvalue)

Output

protein 1:'MKLVVRPWAGCWWSTLGPRGSLSPLGICPLLMLLWATLR'
aveTMValue:0.6