read 
 mol pqr comp17.pqr #System 
 mol pqr Protein17.pqr #Protein 
 mol pqr BEC17.pqr #BEC 
end 

ELEC name comp17 
 mg-auto 
 dime   161 161 193 #System 
 cglen 84.0990 106.8790 110.3130 #System
fglen 49.470 62.870 64.890 #System 
 cgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #System 
 fgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #System 
 npbe 
 bcfl sdh
ion  charge 1 conc 0.1500 radius 2.000 
 ion charge -1 conc 0.1500 radius 2.000 
 pdie 2 #System 
 sdie 80 #System
chgm spl2 
 mol 1 
 srfm smol 
 srad 1.4 #System 
 swin 0.3 #System 
 temp 300 #System 
 sdens 10.0 #System
calcenergy total #comps #total 
 calcforce no 
end 

ELEC name Protein17 
 mg-auto 
 dime   161 161 193 #Protein 
 cglen 84.0990 106.8790 110.3130 #Protein
fglen 49.470 62.870 64.890 #Protein 
 cgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #Protein 
 fgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #Protein 
 npbe 
 bcfl sdh
ion  charge 1 conc 0.1500 radius 2.000 
 ion charge -1 conc 0.1500 radius 2.000 
 pdie 2 #Protein 
 sdie 80 #Protein 
chgm spl2 
 mol 2 
 srfm smol 
 srad 1.4 #Protein 
 swin 0.3 #Protein 
 temp 300 #Protein 
 sdens 10.0 #Protein
calcenergy total #comps #total 
 calcforce no 
end 

ELEC name BEC17 
 mg-auto 
 dime   161 161 193 #BEC 
 cglen 84.0990 106.8790 110.3130 #BEC 
fglen 49.470 62.870 64.890 #BEC 
 cgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #BEC 
 fgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #BEC 
 npbe 
 bcfl sdh 
ion charge 1 conc 0.1500 radius 2.000 
 ion charge -1 conc 0.1500 radius 2.000 
 pdie 2 #BEC 
 sdie 80 #BEC
chgm spl2 
 mol 3 
 srfm smol 
 srad 1.4 #BEC 
 swin 0.3 #BEC 
 temp 300 #BEC 
 sdens 10.0 #BEC
calcenergy total #comps #total 
 calcforce  no 
end 

ELEC name vaccomp17 
 mg-auto 
 dime   161 161 193 #System 
 cglen 84.0990 106.8790 110.3130 #System 
 fglen 49.470 62.870 64.890 #System
cgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #System 
 fgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #System 
 npbe 
 bcfl sdh
ion  charge 1 conc 0.00 radius  2.000 
 ion charge -1 conc 0.00 radius 2.000 
 pdie 2 #System 
  sdie 1 #System
chgm spl2 
 mol 1 
 srfm smol 
 srad 1.4 #System 
 swin 0.3 
 temp 300 #System 
 sdens 10.0 #System
calcenergy total #comps #total 
 calcforce no 
end 

ELEC name vacProtein17 
 mg-auto 
 dime   161 161 193 #Protein 
 cglen 84.0990 106.8790 110.3130 #Protein 
 fglen 49.470 62.870 64.890 #Protein
cgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #Protein 
 fgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #Protein 
 npbe 
 bcfl sdh 
ion  charge 1 conc 0.00 radius 2.000 
 ion charge -1 conc 0.00 radius 2.000 
 pdie 2 #Protein 
 sdie 1 #Protein
chgm spl2 
 mol 2 
 srfm smol 
 srad 1.4 #Protein 
 swin 0.3 #Protein 
 temp 300 #Protein 
 sdens 10.0 #Protein
calcenergy total #comps #total 
 calcforce no 
end 

ELEC name vacBEC17 
 mg-auto 
 dime   161 161 193 #BEC 
 cglen 84.0990 106.8790 110.3130 #BEC 
 fglen 49.470 62.870 64.890 #BEC
cgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #BEC 
 fgcent 62.85500000000000000000 62.82500000000000000000 29.64500000000000000000 #BEC 
 npbe 
 bcfl sdh
ion  charge 1 conc 0.00 radius 2.000 
 ion charge -1 conc 0.00 radius 2.000 
 pdie 2 #BEC 
 sdie 1 #BEC
chgm spl2 
 mol 3 
 srfm smol 
 srad 1.4 #BEC 
 swin 0.3 #BEC 
 temp 300 #BEC 
 sdens 10.0 #BEC
calcenergy total #comps #total 
 calcforce no 
end 

print elecEnergy comp17 - Protein17 - BEC17 - vaccomp17 + vacProtein17 + vacBEC17 end 
quit
