用python写的处理PDB的文件

def Deal_pdb(fileA,fileB,index,parameterA,replaceParaA,line_paraA,parameterB,replaceParaB,line_paraB):
 from_file=open(fileA,"r")
 to_file=open(fileB,"a")
 lines=from_file.readlines()
 from_file.close()
 for line in lines:
  line_data=line.split()
  print len(lines)
  if line_data[index]==parameterA and int(line_data[1])>line_paraA:
   line_data[index]=replaceParaA
   line=line.replace(" %s "%parameterA," %s "%replaceParaA)
   to_file.write(line)
  elif line_data[index]==parameterB and int(line_data[1])>line_paraB:
   line_data[index]=replaceParaB
   line=line.replace(" %s "%parameterB," %s "%replaceParaB)
   to_file.write(line)
  else:
   to_file.write(line)
 to_file.close()

 

 

if  __name__=="__main__":
 Deal_pdb("3BEQ-tet.pdb","1.pdb",4,"A","C",2983,"B","D",9856)