Getting same total SASA results with different given probes

[ANUJGAHLAWAT]

CODE:

freesasa.Parameters({'probe-radius' : Probe}) freesasa.setVerbosity(freesasa.nowarnings) structure = freesasa.Structure("G:/Anuj/delaulney trangulation/check different PDB/"+File+".pdb",{'hetatm' : True, 'hydrogen' : True}) result = freesasa.calc(structure) print("pdb sasa :",result.totalArea())

PROBLEM: I am getting the same results i.e. 17478.4794656937 for 3o7p.pdb using the result.totalArea() function even if I am giving different probes (1.4 A and 3.0 A) to the freesasa.Parameters function. Generally, SASA for the PDB file should be lesser with a higher radius probe than that of a small probe.

Please provide your useful comments Thank you

[mittinatten]

Hi, thanks for reaching out. That sounds like a bug. I'll have a look.

[mittinatten]

I checked this out and found the problem, you need to pass the parameters to the calc function. (in your code you create a parameters object and then discard it, so it won't be used by the calculation). Here's code that works for me:

structure = freesasa.Structure("3o7p.pdb")
result = freesasa.calc(structure, freesasa.Parameters({'probe-radius': 3.0 }))

[ANUJGAHLAWAT]

Thank you for your response 😊

On Thu, 28 Apr 2022, 9:52 pm Simon Mitternacht, @.***> wrote:

Closed #77 https://github.com/mittinatten/freesasa/issues/77.

— Reply to this email directly, view it on GitHub https://github.com/mittinatten/freesasa/issues/77#event-6516563012, or unsubscribe https://github.com/notifications/unsubscribe-auth/AIAFUKCAVG2JFH7THMUOOMTVHK3MPANCNFSM5UCBMAHQ . You are receiving this because you authored the thread.Message ID: @.***>